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ABSTRACT 

A  quantitative  analysis  was  carried  out  on  the  performance  of 
turboprop  aircraft  within  a  microburst  windshear.  The  objective  of 
the  analysis  was  to  provide  specific  flight  procedures  for  optimal 
navigation  through  the  windshear.  The  microburst  windshear  model 
used  in  the  analysis  embodied  the  severe  characteristics  of  the 
microburst  encountered  by  Delta  Flight  191  during  an  approach  to 
landing  at  Dallas/Ft.  Worth,  2  August,  1985.  Different  escape 
strategies  were  tested  using  the  flight  performance  characteristics  of 
the  U.S.  Navy's  P-3  "Orion"  and  T-44  "Pegasus"  aircraft.  The  three 
flight  phases  investigated  were  approach  to  landing,  takeoff,  and  the 
low  altitude  ASW  mission.  Results  from  the  analysis  were  coupled 
with  the  pilot's  view-point  from  which  conclusions  were  drawn.  The 
results  of  the  analysis  support  a  constant-pitch-angle  escape 
procedure.  The  same  procedural  steps  can  be  used  for  both  aircraft 
in  any  configuration  or  situation  with  the  difference  being  the  degree 
of  pitch  to  employ.  The  conclusions  are  in  a  format  for  integrating 
specific  microburst  escape  procedures  within  the  NATOPS  programs 
for  the  P-3  and  T-44. 
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I.    INTRODUCTION 

The  microburst  windshear  is  documented  as  a  serious  hazard  to 
aviation.  The  microburst  has  been  a  contributing  factor  in  several 
airline  accidents  in  the  past  two  decades.  The  seriousness  of  these 
events  has  lead  to  a  significant  amount  of  research  on  the  effect  of 
microbursts  upon  airline  transport  aircraft.  However,  little  has  been 
documented  on  the  effect  of  a  microburst  on  light-to-medium  weight 
turbopropeller  aircraft.  The  U.S.  Navy  operates  a  fleet  of  P-3  "Orion" 
turboprop  aircraft  in  all-weather,  long-range  patrol  missions.  The 
P-3  aviation  community  requested  a  study  of  the  effects  of  a 
microburst  upon  the  aircraft's  dynamic  flight  performance  and  the 
appropriate  escape  procedure  if  one  were  encountered.  Training  of 
the  P-3  pilot  starts  in  the  T-44,  a  Beechcraft  King  Air  twin  turboprop. 
Therefore,  the  ultimate  goal  of  this  thesis  was  to  derive  viable 
microburst  windshear  escape  procedures  which  can  be  incorporated 
into  the  Naval  Air  Training  and  Operating  Procedures  Standardization 
(NATOPS)  Program  for  both  the  P-3  and  the  T-44  aircraft. 

A.  MICROBURSTS  AND  AVIATION 
1.  Microburst  Definition 

The     microburst     is     a     form     of     downburst     windshear. 
Downbursts  are  associated  with  strong  convectlve  activity  (namely 
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thunderstorms).    Downbursts  are  formed  by  a  cold  air  mass  dropping 

from   within    a    thunderstorm    producing   strong   downdrafts.      The 

downdraft  turns  outward  as  it  approaches  the  earth.     Fujita  is  a 

prominent    pioneer    in    identifying    and    classifying    this    type    of 

windshear  [Ref.  1].    He  defines  the  microburst  as  [Ref.  2:p.  8]: 

A  small  downburst  with  its  outburst,  damaging  winds 
extending  only  4  km  (2.5  miles)  or  less.  In  spite  of  its  small 
horizontal  scale,  an  intense  microburst  could  induce  damaging 
winds  as  high  as  75  m/sec  (168  mph). 

Conclusions  drawn  from  the  Northern  Illinois  Meterological  Research 

on  Downbursts  (NIMROD)  and  Joint  Airport  Weather  Studies  (JAWS) 

[Ref.  2]  programs  show  that  the  microburst  structure  consists  of  a 

down    flow    with    a    structure    of    one    or    more    vortex    rings    as 

illustrated    in    Figure  1.       The    microburst    may    be    wet    or    dry; 

precipitation  may  or  may  not  be  present  in  the  downdraft.     The 

design    of    a    microburst    structure   was   not    understood,    or    even 

accepted  as  a  viable  phenomenon  prior  to  1982.     Through  the  1980s, 

meterological  researchers  and  airline  accident  investigators  came  to 

the    conclusion    that    the    microburst    weather    phenomenon    is    a 

significant  hazard  to  aviation  and  has  been  a  contributing  factor  in 

many  airline  accidents.    The  accumulation  of  the  NIMROD  and  JAWS 

project  results  and  Inflight  data  conclusively  shows  that  ring  vortex 

flow  is  prevalent  in  microburst  windshear  [Ref.  3]. 


Microburst  windshear  is  uniquely  different  from  other  types  of 
windshear.  The  other  prevalent  type  of  windshear  hazard  to 
aviation  is  the  thunderstorm  gustfront  [Ref.  4].  A  gustfront 
windshear  is  associated  with  massive  settling  of  rain-cooled  air 
resulting  from  thunderstorm  activity.  A  sudden  wind-shift  with 
gusty  conditions  that  precedes  a  thunderstorm  is  characteristic  of  a 
gustfront.  The  main  difference  between  a  gustfront  and  a 
microburst  is  the  size  of  the  area  effected.    A  gustfront  is  widespread 


Figure     1. 
microburst. 


Wind     streamlines     characteristic     of     a 


across  the  face  of  one  or  more  thunderstorms.  A  microburst  is  small 
and  isolated  within  convective  activity.  The  local  magnitude  of  the 
horizontal  and  vertical  windshear  is  what  matters  to  an  aircraft 
penetrating  convective  activity.  Although  smaller,  a  microburst  can 
be  much  more  of  a  detriment  to  an  aircraft's  flight  performance. 


2.  The  Effect  Of  Microbursts  On  Aviation 

The  aviation  community  has  always  been  aware  of  the 
potential  danger  of  windshear.  Classification  of  different  types  was 
not  addressed  until  1976,  when  The  Federal  Aviation  Administration 
(FAA)  published  an  advisory  circular  discussing  the  effects  of 
windshear  on  flight  performance.  It  was  updated  in  1979  [Ref .  5] 
where  the  term  "downburst"  was  used  to  describe  a  strong 
downdraft  associated  with  the  center  of  a  thunderstorm.  There  was 
considerable  debate  on  the  hazard  of  the  downburst.  Many  argued 
that  the  downdraft  rapidly  weakened  as  it  approached  the  ground. 
The  research  by  Fujita  provided  overwhelming  evidence  to  the 
contrary,  thereby  defining  the  microburst.  By  the  mid-eighties,  the 
term  "microburst"  was  readily  being  used  by  the  aviation 
community.  However,  the  potential  effect  of  a  microburst  upon 
aircraft  performance  was  not  readily  understood  by  aircrew. 

The  danger  of  a  microburst  became  very  apparent  following 
the  crash  of  Delta  Flight  191  (DAL  191)  at  Dallas/Ft.  Worth,  August  2, 
1985.  The  crash  ensued  following  the  penetration  of  a  thunderstorm 
during  the  approach  to  landing.  The  aircraft  involved  was  a 
Lockheed  L-1011,  a  three  engine  heavy  airline  transport  aircraft.  The 
National  Transportation  Safety  Board  (NTSB)  investigating  the  DAL  191 
accident  published  in  their  report  the  following    [Ref.  6:p.  81]: 


The  National  Transportation  Safety  Board  determines  the 
probable  causes  of  the  accident  were  the  flight  crew's  decision 
to  initiate  and  continue  the  approach  into  a  cumulonimbus 
cloud  which  they  observed  to  contain  visible  lightning;  the  lack 
of  specific  guidelines,  procedures,  and  training  for  avoiding  and 
escaping  low-altitude  windshear;  and  the  lack  of  definitive,  real- 
time windshear  hazard  information.  This  resulted  in  the 
aircraft's  encounter  at  low  altitude  with  a  microburst-induced, 
severe  windshear  from  a  rapidly  developing  thunderstorm 
located  on  the  final  approach  course. 

Several    pertinent    questions    arise    from    this    conclusion.       Airlines 

routinely  fly  through  "bad"  weather  every  day.     How  do  a  flight 

crew  and  Air  Traffic  Control  (ATC)  authorities  recognize  microburst 

type  weather?      Can  any  aircraft  navigate  successfully   through  a 

microburst    or    is    it    beyond    the    performance    capabilities    of    the 

average  airliner?     Is  there  a  significant  departure  between  normal 

piloting    techniques    and    microburst    escape    procedures?    A    large 

amount    of    research    has    been    done    to    answer    these    questions. 

However,  the  focus  has  been  directed  predominantly  toward  airline 

transport  aircraft. 

The   advent   of   Digital   Flight  Data   Recorders  (DFDR)   provided 

conclusive  evidence  of  the  structure  and  the  potential  effects  of  a 

microburst  windshear  on  aircraft.      The  DFDR   records  the  last  30 

minutes  of  voice  transcripts,  as  well  as  several  aircraft    performance 

parameters.     Aircraft  body  axis  angles  and  accelerations,  airspeed, 

altitude,  and  time  are  a  few  of  the  recorded  parameters.     These 

records  allow  the  reconstruction  of  the  wind  field  vectors  encountered 


by  the  aircraft  [Ret.  7].  DAL  191  was  equipped  with  a  DFDR. 
American  Airline  Flight  539  (A A  539)  followed  DAL  191  and  penetrated 
the  same  convective  weather  110  seconds  later  [Ref.  6]  as  depicted  in 
Figure  2.  This  second  aircraft  was  a  McDonnell/Douglas  MD-80,  also 
equipped  with  a  DFDR.  The  availability  of  digital  flight  data  allowed 
a  comprehensive  database  to  be  built  on  the  wind  field  generated  by 
that  microburst. 


Figure  2.     Flight  paths  of  the  two  airlines  through  the 
Dallas/Ft.  Worth  microburst. 


The  center  of  the  thunderstorm  activity  was  located 
approximately  10,300  feet  from  the  approach  end  of  Runway  17L  at 
Dallas/Ft.  Worth  Airport.  The  DFDR  data  show  that  DAL  191  initially 
encountered  a  headwind  followed  by  a  very  abrupt  tailwind  [Ref.  7]. 


retrospect,     several     previous     airline    accidents     can    be 


A  significant  downward  airflow  also  accompanied  the  horizontal 
windshear  which  resulted  in  a  loss  of  airspeed  accompanied  by  a 
sharp  nose  down  pitch.  A  very  high  sink  rate  resulted.  The  aircraft 
exited  the  bottom  of  the  thunderstorm  approximately  100  feet  above 
ground  level  (AGL).  Ground  impact  ensued  with  the  aircraft  in  a 
slight  nose  up  attitude  and  at  a  high  airspeed  [Ref  5]. 

In 

attributed  to  windshear  of  this  nature.  Fujita  had  analyzed  several 
takeoff  and  landing  airline  accidents  involving  windshear  prior  to  the 
DAL  191  crash  [Ref.  2].  However,  the  Dallas/Ft.  Worth  accident 
showed  the  need  to  address  the  microburst  hazard  and  gain  a  more 
thorough  understanding  of  how  to  alleviate  the  danger. 
3.  Reducing  The  Hazard  To  Aviation 

One  approach  to  reducing  the  hazard  of  a  microburst  on 
aviation  is  recognition  and  avoidance.  This  solution  is  very  difficult 
given  available  aviation  weather  avoidance  equipment.  The  primary 
windshear  equipment  available  today  is  the  aircraft  weather  radar 
and  ground  based  Low  Level  Windshear  Alerting  System  (LLWAS). 
On-board  weather  radar  provides  a  picture  of  the  amount  of 
convective  activity  ahead,  but  cannot  recognize  microburst  activity. 
Ground  based   LLWAS  provides  warning  when  wind   speed  sensors 


located  about  the  airfield  measure  significantly  different  wind 
vectors.  However,  the  microburst  must  occur  at  the  location  of  the 
LLWAS  system.  An  LLWAS  system  was  installed  at  the 
Dallas/Ft.  Worth  airport  at  the  time  of  the  DAL  191  accident.  The 
microburst  occurred  on  the  final  landing  approach  corridor  two  miles 
from  the  center  of  the  airport.  Because  of  the  location  of  the  LLWAS 
sensors  around  the  airport  perimeter,  the  system  did  not  give  any 
windshear  warning  prior  to  the  accident.  An  LLWAS  alert  was 
sounded  as  the  thunderstorm  passed  overhead  the  airfield  several 
seconds  after  DAL  191  crashed.  Neither  weather  radar  nor  LLWAS 
provide  real-time  information  on  microburst  windshear. 

Very  recent  technology  has  produced  some  viable,  and 
expensive,  means  to  recognize  and  avoid  related  microburst  weather 
activity.  Ground  based  Doppler  Radar  has  been  successfully  proven 
at  Denver's  Stapleton  airport  [Ref.  8].  On-board  forward  looking 
devices  such  as  infrared  sensors,  lasers,  and  doppler  radar  are  being 
developed  to  provide  several  seconds  of  warning  prior  to  a 
microburst  encounter.  But  until  avoidance  systems  are  made  readily 
available,  survival  of  a  microburst  encounter  will  be  dependent  on 
the  use  of  a  successful  escape  procedure  by  the  flight  crew. 
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The  second  approach  in  reducing  the  microburst  hazard  is 
developing  viable  escape  procedures.  Research  has  shown  that  quick 
recognition  of  a  microburst  encounter  is  paramount.  The  reaction 
time  of  a  flight  crew  coupled  with  the  type  of  escape  maneuver  will 
govern  the  success  of  escape  [Ref.  9].  The  escape  procedure  the  flight 
crew  executes  depends  on  the  flight  phase  and  configuration  of  the 
aircraft.  Much  research  has  been  performed  concerning  aircraft 
performance  and  escape  procedures  with  the  focus  on  airline 
transport  aircraft  operations  to  be  described  below. 

B.  MICROBURST  ESCAPE  PROCEDURES 

A  large  amount  of  research  has  been  done  on  defining  the  optimal 
escape  maneuvers  for  airline  transport  aircraft.  Much  of  that 
research  is  based  upon  control  theory.  Comparisons  of  different  flight 
strategies  for  microburst  encounters  have  been  addressed  [Ref.  10, 
Ref.  11].  The  total  energy  concept  is  readily  used  to  compare  the 
outcomes  of  different  maneuvers  [Ref.  12,  Ref.  13].  One  very 
important  aspect  in  determining  the  optimal  escape  procedure  is  the 
ability  of  a  flight  crew  to  execute  the  appropriate  maneuver  utilizing 
the  available  flight  performance  information  available. 


The  FAA  has  generated  an  exhaustive  Windshear  Training  Aid 
[Ref.  14]  aimed  at  modern  day  transport  aircraft.  It  addresses  the 
crew  training  requirements  and  suggests  viable  microburst  escape 
strategies.  Analytical  research,  as  well  as  tried  and  tested 
procedures,  have  verified  the  suggestions  presented  by  the  FAA.  All 
the  major  airlines  have  incorporated  microburst  recognition  and 
escape  into  their  recurrent  training.  The  majority  of  the  civil  flight 
simulators  possess  some  degree  of  windshear  emulating  a  microburst. 
As  of  today,  the  aviation  community  recognizes  the  microburst  as  a 
aviation  hazard  and  has  taken  significant  steps  in  reducing  the 
vulnerability. 

The  different  escape  strategies  fall  within  three  general  categories: 

1)  Optimum  aircraft  performance  through  the  airmass. 

2)  Optimum  aircraft  performance  related  to  an  inertial  reference. 

3)  A  combination  of  the  two. 

Optimum  aircraft  performance  Is  defined  as  the  performance 
providing  the  largest  cushion  of  escape  from  a  microburst  encounter. 
The  parameters  measured  for  evaluating  an  escape  strategy's 
outcome,  are  airspeed  and  altitude  gained  (or  lost). 

The  first  category  includes  maneuvers  that  trade  altitude  for 
airspeed.    The  second  category  uses  measured  data  such  as  body  axis 
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acceleration  and  windshear  magnitude  in  calculating  the  optimal 
flight  path  [Ref.  12,  Ref.  10,  and  Ref.  11].  Obviously  the  second 
category  requires  significant  airborne  computing  capabilities  not 
normally  available.  The  FAA's  Windshear  Training  Aid  recommends 
an  escape  strategy  that  falls  within  the  third  category,  The 
recommendation  is  to  maintain  a  constant  pitch  attitude  with 
maximum  engine  thrust  applied.  In  general,  increasing  pitch 
attitude  toward  &  15°  pitch  angle  shows  optimal  performance  on  most 
of  the  large  airline  transport  aircraft.  The  constant  pitch  attitude  is 
maintained  with  disregard  to  airspeed. 

C.  EFFECTS  OF  A  MICROBURST  ON  FLIGHT  INSTRUMENTS 

How  a  flight  crew  perceives  and  navigates  through  a  microburst 
encounter  is  strongly  influenced  by  their  flight  instrument 
indications.    The  primary  flight  Instruments  available  to  pilots  are: 

1)  Attitude  Indicator  (AI)  -  a  gyroscopic/inertial  stabilized  aircraft 
attitude  indicator.  Modern  design  allows  accurate  attitude 
information  in  the  highest  levels  of  turbulence. 

2)  Airspeed  Indicator  -  a  pltot/statlc  system  instrument  which 
displays  airspeed  in  knots.  It  is  sensitive  to  airmass  pressure 
changes  and  position  error. 

3)  Pressure  Altimeter  -  a  static  system  instrument  which  relates 
altitude,  in  feet,  to  static  pressure.  It  is  sensitive  to  airmass 
pressure  changes. 
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4)  Vertical  Speed  Indicator  (VSI)  -  a  static  system  instrument 
which  measures  the  change  in  static  pressure  in  feet  per 
minute.  It  is  unreliable  in  turbulence  due  to  mechanical  and 
pressure  lag. 

All  aircraft  have  these  primary  instruments.    Inertial  navigation  and 

true    airspeed    computers    do    provide    other    flight    performance 

references  usually  through  flight  directors.    However,  the  majority  of 

the  flight  stations  rely  primarily  on  the  basic  flight  instruments. 

Microbursts  are  always  associated  with  some  degree  of  turbulence. 
Also,  pressure  differentials  can  be  expected  within  convective 
activity.  Ground  and  airborne  data  show  that  the  pressure  within 
a  microburst  varies  about  ±  3  millibar  [Ref.  2,  Ref.  3].  This  pressure 
change  equates  to  about  ±  80  feet  in  altimeter  variation. 

The  attitude  indicator  will  be  the  flight  Instrument  least  affected 
by  microburst  atmospherics.  This  instrument  will  always  give 
reliable  aircraft  attitude  information.  Due  to  the  minor  static 
pressure  changes,  the  airspeed  indicator  and  altimeter  can  still  be 
relied  upon  as  performance  Instruments.  The  airspeed  indicator  will 
be  one  of  the  primary  instruments  to  indicate  a  windshear 
penetration.  The  VSI  would  probably  be  erratic  and  unreliable  due 
to  turbulence  and  instrument  lag. 
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There  exits  a  means  of  measuring  angle-of-attack  (AOA)  in  most 
larger  aircraft.  Transport  aircraft  utilize  AOA  as  a  stall  warning 
device  (stick  shaker).  Navy  aircraft  incorporate  a  means  to  read  AOA 
in  units.  AOA  measuring  devices  may  have  a  high  damping  factor  in 
their  measurement,  resulting  in  substantial  lag  during  turbulent 
conditions.  Therefore,  AOA  indication  may  not  be  a  reliable 
performance  indictor  during  a  microburst  encounter. 

D.  PROBLEM  STATEMENT 

The  P-3  has  no  means  to  recognize  and  then  avoid  microburst 
windshear.  It  must  rely  solely  on  escape.  The  published  guidelines 
for  microburst  escape  are  based  upon  airline  transport  aircraft  which 
significantly  differ  from  the  P-3.  In  this  study,  the  results  of 
different  escape  procedures  based  upon  the  available  flight 
Instruments  are  compared  using  the  parameters  of  altitude,  airspeed, 
and  specific  energy.  Characterizing  the  effect  of  a  microburst 
windshear  upon  the  aircraft's  performance  followed  by  the 
implementation  of  the  appropriate  escape  procedure  form  the 
objective  for  the  research.  The  final  conclusions  must  provide  viable 
escape  procedures  that  are  based  upon  the  primary  flight 
instruments. 
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The  microburst  windshear  model  selected  must  conform  to  a 
known,  measured  phenomenon.  The  windshear  experienced  by  DAL 
191  provided  a  suitable  database  for  emulation. 

A  series  of  equations  of  motion  for  the  aircraft  were  developed 
and  designed  to  be  controlled  through  pitch  angle  and  thrust  inputs. 
These  two  parameters  emulate  the  pilot's  available  flight  control 
inputs. 

Three  flight  phases  critical  to  a  microburst  encounter  were 
considered.  Approach  to  landing,  takeoff,  and  on-statlon  (P-3  only) 
phases  require  operation  low  to  the  ground  in  all-weather  conditions. 
Viable  escape  procedures  using  attitude,  airspeed,  altitude,  and  angle 
of  attack  were  analyzed  for  all  three  flight  phases. 

Three  questions  were  considered  for  each  flight  phase: 

1.  What  is  the  optimum  microburst  escape  procedure  given 
available  flight  Information  to  the  pilot? 

2.  Does  the  optimum  escape  procedure  change  with  gross  weight 
or  available  thrust? 

3.  What  flight  instrument  indications,  flight  control  feeling,  and 
dynamic  response  would  be  expected  during  the  optimum  escape 
maneuver? 

Viable  microburst  escape  procedures  were  derived  for  the  P-3  and 

T-44  aircraft. 
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II.    MATH  MODELS  AND  EXPERIMENTAL  DESIGN 

This  analysis  is  based  upon  mathematical  representations  of 
microburst  windshear  and  aircraft  performance.  Microburst 
windshear,  inertial  reference,  aircraft  performance,  and  critical  flight 
phases  form  the  major  segments  of  the  total  algorithm  suite.  Most 
of  the  mathematical  simulation  performance  was  carried  out  on  a 
Macintosh  SE®  computer  utilizing  MatLab®  software. 
A.  MICROBURST  WINDSHEAR  MODEL 

A  double  ring  vortex  was  chosen  to  emulate  microburst 
windshear.  Fujita  [Ref.  2:p.  14]  demonstrated  that  the  microburst 
was  a  settling  airmass  with  vortices  generated  near  the  earth.  This 
phenomenon  produces  varying  degrees  of  three-dimensional 
windshear  when  close  to  the  earth's  surface.  Figure  1  showed  the 
general  concept  of  a  microburst. 

A  strong  microburst  windshear  was  experienced  at  Dallas/Ft. 
Worth  on  August  2,  1985.  A  Delta  Airlines  L-1011  (DAL  191)  crashed 
while  on  final  approach  during  a  microburst  windshear  encounter. 
An  American  Airlines  MD-80  (AA  539)  performed  a  missed  approach 
110  seconds  following  the  DAL  191  crash.  Both  aircraft  were  equipped 
with  DFDRs  which  register  a  time  history  of  the  aircraft's 
parameters.       Body    axis    accelerations,   velocities,    and    Euler-angle 
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values  were  used  to  calculate  the  flow-field  winds  generated  from 
the  microburst.  Wingrove  and  Bach  [Ref  7]  provided  the  analytical 
means  to  calculate  the  windshear.  From  these  data,  an  insight  was 
gained  of  the  magnitude  and  characteristics  of  the  windshear. 

A  double-vortex-ring  model  was  chosen  for  the  windshear  model. 
Schultz  [Ref.  15]  devised  a  double-vortex-ring  mathematical  model 
which  closely  approximated  the  windshear  experienced  by  the  AA  539 
flight.  The  aircraft  flew  through  the  microburst  windshear  at  2500 
feet  AGL.  The  original  algorithm  was  developed  and  tested  against 
recorded  flight  data.  The  double-vortex-ring  algorithm  was  modified 
in  the  current  study  with  the  addition  of  source  flows  to  better  fit 
low  altitude  windshear.  This  modified  windshear  model  was  fit  to 
the  DAL  191  flight  data  by  the  application  of  a  parameter  sensitivity 
scheme.  Ring  location  and  vortex  Intensity  were  found  to  differ 
from  the  A  A  539  model.  However,  a  multiple  microburst  structure 
can  be  expected  during  strong  convective  activity  [Ref.  2:p.  35]. 
MatLab®  programing  was  utilized  in  this  study  and  is  listed  In 
Appendix  A. 

1.  Vortex  Ring  Model  Development 

The  development  of  this  vortex  ring  model  comes  directly  from 
the  work  of  Schultz  [Ref.  15].  Source  flow  was  integrated  into  this 
algorithm  to  more  closely  match  the  DAL  191  windshear.  The  wind 
component   in  the   y  direction  was  not   considered   in  the  current 
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model.  The  y  component  is  not  easily  modeled  by  symmetric  vortex 
rings  [Ref.  15]  and  has  no  direct  impact  on  the  longitudinal  dynamics 
of  an  aircraft.  Figure  3  illustrates  the  geometric  aspects  of  the 
vortex  ring  model.  Imaging  of  the  vortex  rings  leads  to  in  viscid 
ground  effect  upon  the  wind  field.    The  wind  field  becomes  horizontal 


Figure  3.    Vortex  ring  mathematical  representation, 
next  to  the  ground  plane  as  seen  In  nature.    Boundary  layer  effects 
were  neglected  in  this  model.    Bowles  and  Oseguera  have  performed 
research   in   boundary   layer   effects  on   windshear   [Ref.  16].      From 
Figure  3,  the  following  relationships  can  be  seen  to  hold: 
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where  G,  P,  and  I  are  ground,  primary,  and  image  reference 
respectively.  The  velocity  components  for  each  vortex  were  obtained 
by  differentiating  the  vortex  ring  stream  function  <J>: 
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where  r  is  the  radial  distance  from  the  z  axis: 
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The  stream  function  for  a  vortex  ring  filament  was  obtained 
from  the  evaluation  of  an  elliptical  Integral  [Ref.  15]: 


0.788X 
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where  n  and  r2  represent  the  closest  and  farthest  distances  to  the 
point   of  interest   from  the  ring's  filament;  k  is   the  ring  filament 
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strength;  and  A  is  a  scaling  term.    Their  algebraic  relationships  are  as 
follows: 


r  -Vz^Cr-li 
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where  R  is  the  radius  of  the  filament. 

Equation  6  was  differentiated  to  obtain  expressions  for  3r   and 

3z: 
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Equations  10  and  11  were  substituted  into  Equation  4  to  get  the 
individual  velocity  components  produced  by  each  ring  at  a  given 
point  in  space  (with  respect  to  the  ground  reference). 

A  velocity  damping  factor  (Q  was  calculated  for  each  ring  to 
prevent  erroneously  high  values  of  velocity  near  the  filament  cores. 
The  value  ranged  from  zero  to  one;  it  is  zero  when  the  point  of 
interest  is  at  the  viscous  core  and  approaches  one  at  increased 
distances.     The  algebraic  relationship  for  C  is  as  follows  [Ref.  15]: 


C .  ■  (1  -  exp 
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(19) 
where  ru  is  the  closest  distance  from  the  point  of  interest  to  the  /th 

ring's  viscous  core.  Schultz  obtained  the  weighting  factors,  I  and  p, 
by  a  visual  comparison  of  the  vortex  ring  strength  and  through  a 
parameter  estimation  scheme  [Ref.  15].  A  total  damping  factor  (Z)  is 
obtained  for  the  point  of  Interest  from  the  product  of  the  damping 
factors  for  each  of  the  four  vortex  rings: 
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i-1  (20) 

The  low  altitude  windshear  model  required  the  addition  of  two 

point  source  flows  located  10,000  feet  above  and  below  the  ground 

plane  at  the  center  of  the  microburst.    The  radial  flow  is  represented 

by  the  following  expression: 


r    2,Rj  Q1) 

where  A  is  the  source  strength  and  Rj  is  the  radial  distance  from  the 

point  source  to  the  point  of  interest.     The  horizontal  and  vertical 
wind  vectors  contributions  are  estimated  by: 
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Combining  the  upper  and  lower  (mirror  image)  sources  produce: 
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2  1 

In  summary,  wind  velocity  for  a  point  of  interest  In  space 

referenced  to  earth  was  calculated  in  the  following  manner: 

1)  X  (horizontal)  and  z  (altitude)  positions,  measured  in  feet, 
are  inputs  for  Equation  4.  Calculations  are  made  for  the 
four    vortex    rings    (large  primary,    small    primary,    large 
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image,  and  small  image).  A  horizontal  wind  velocity  (vx*) 
and  a  vertical  wind  velocity  (vzj)  are  calculated  for  each 
ring. 

2)  A   damping  factor   ty)  was  calculated   for  each  ring   using 
Equation  19. 

3)  The    induced    velocities    for    each    ring    were    added    and 
multiplied  by  the  total  velocity  damping  factor  (Z): 
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where  Wx  and  Wh  are  windshear  values  (ft/s). 

4)  For  the  low  altitude  model,  the  source  flow  velocities  were 
added  to  the  vortex  velocities  to  give  the  final  horizontal 
and  vertical  windshear  components. 

2.  Windshear  Model  Fit 

The  unmodified  double  vortex  windshear  model  was  applied  to 
AA  539  flight  data  by  Schultz,  using  a  parameter  minimization 
scheme  to  fit  the  double-vortex-ring  model  to  the  actual  DFDR  data. 
Table  1  lists  the  parameters  varied  and  the  final  results.  Schultz 
concluded  that  the  difference  between  the  model  results  and  recorded 
wind  data  was  an  rms  (root  mean  square)  of  16  ft/s  for  the  entire 
system.  Figure  4  shows  the  comparison  of  recorded  winds  to  the 
model  prediction. 
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The  AA  539  windshear  model  was  fit  in  the  current  study  for 
model  verification  only.       The    windshear  that    DAL  191  encountered 
was  the  windshear  to  emulate.    This    windshear  model    was    used  to 
TABLE  1.    2500  FEET  ALTITUDE  WINDSHEAR  MODEL  PARAMETERS 


parameter 

large  ring 

small  ring 

ring  radius,  ft 

8503.3 

1701.7 

core  radius,  ft 

20004.1 

323.9 

ring  circulation,  ft'2/s 

431968.8 

57204.9 

x  position,  ft 

0.0 

50.0 

V  position,  ft 

3350.4 

830.9 

z  position,  ft 

3400.6 

2333.6 

analyze    aircraft    performance    during    takeoff    and    landing    flight 

segments.    Therefore,  the  double-vortex-ring  model  was  modified  by 

th  parameters  listed  in  Table  2.    The  data  from  DAL  191  showed  that 

TABLE  2.    LOW  ALTITUDE  WINDSHEAR  PARAMETERS 


parameter 

large  ring 

small  ring 

ring  radius,  ft 

7000 

1300 

core  radius,  ft 

2004.1 

323.9 

ring  circulation,  ft*2/s 

431968.8 

131571.3 

x  position,  ft 

2500.0 

300.0 

v  position,  ft 

-300.0 

1.0 

z  position,  ft 

3400.6 

800.0 

x  dir.  source  strength 

1355396  fT2/s 

z  dir.  source  strength 

3049200  fT2/s 

the  two  rings  were  not  concentric  about  the  center  [Ref.  3].  The 
small  vortex  ring  was  lower  to  the  ground  and  displaced  in  the  x 
direction  from  the  large  ring.  Point  source  flows  located  10,000  feet 
above  and  below  the  surface  plane  were  added  to  the  double-vortex- 
ring  model  to  increase  the    outflow.    A  simple   parameter    sensitivity 
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Figure    4.       Vortex   ring    model   comparison    to   recorded 
flight  data  [Ref.  3]. 
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scheme  was  then  used  to  estimate  the  parameters  that  best  fit  the 
modified  windshear  model  to  the  recorded  DAL  191  flight  data. 

The  sensitivity  scheme  (listed  in  Appendix  B)  used  an  iterative 
process  of  varying  each  parameter  individually,  then  calculating  the 
x  direction,  z  direction,  and  total  system  rms  difference  between 
model  and  recorded  winds.  The  initial  parameter  matrix  was  built 
up  as  an  8x3  matrix.  The  center  column  was  composed  of  the  "best 
guess"  values.  Columns  one  and  three  were  values  obtained  from  an 
interval  surrounding  the  initial  guess  value.  First  guess  parameters, 
including  the  corresponding  intervals,  were  obtained  by  graphical 
comparison  of  model  and  recorded  data.  The  parameter  matrix  was 
refined  after  each  successive  run  of  the  sensitivity  scheme.  Final 
parameter  resolution  was  less  than  3S  change  for  the  total  system 
rms.  Table  2  lists  the  best-fit  parameters  to  the  DAL  191  flight 
windshear  measurements.  Figure  5  is  a  graphical  representation  of 
the  current  model's  vortex  rings  and  source  flow  point  relative 
locations.  The  wind  vectors  produced  from  the  model  are  shown  in 
Figure  6.  Figure  7  compares  the  model  windshear  with  the 
measured  wind  velocities  from  DAL  191. 

Total  system  rms  was  calculated  as  16.5  ft/s  for  the  microburst 
windshear  model.  The  final  model  shows  a  close  approximation  of  an 
actual  microburst,  serving  as  a  good  windshear  model  for  testing 
aircraft  performance  during  low-level  encounters. 
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Figure  5.    Geometric  relationship  of  the  vortex  rings  and 
point  source  within  the  microburst  model. 
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Figure  6.   Vector  plots  of  the  wind-fields  generated  by 
the  microburst  windshear  model. 
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Figure  7.    DAL  191  recorded  winds  compared  to 
mathematical  model. 
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B.  INERTIAL  REFERENCE  FRAME 

The  concept  of  aircraft  specific  energy  was  used  to  compare  the 
outcomes  of   different  windshear   escape  strategies.      The  aircraft's 
total  energy  is  defined  as  the  sum  of  the  airmass-relative  kinetic 
energy  and  the  inertial  potential  energy.    A  vectorial  relationship  of 
an  aircraft's  motion  through  the  airmass  is  presented  in  Figure  8.    A 
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Figure  8.    Inertial  reference  coordinate  system, 
set   of   coupled,   non-linear   differential   equations   can   be   developed 
from  this  vectorial  relationship.     These  differential  equations  can  be 
solved  numerically  where  specific  energy  can  be  determined  at  any 
point  in  Inertial  space. 
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1.  Total  Energy  Concept 

The  aircraft's  altitude  and  airspeed  at  any  given  point  can  be 
used  to  determined  total  energy  from  the  following  relationship: 

E  =  &mV2  +  mgh  (26) 

where  m  is  mass,  V  is  airspeed,  and  h  is  the  altitude.     The  specific 
energy  (sometimes  called  energy  height),  Es,  Is  defined  by: 

Es  ■  -=—■  +  h  ,    . 

2g  (27) 

The   time   rate   of   change   of    specific   energy    is   equal   to   aircraft 

acceleration   plus  rate  of   climb.      Differentiating  Equation  (27)  with 

respect  to  time  gives  specific  power: 


,v(v)  +  h 


(28) 

Values  for  V,  V,  h,  and  h  can  be  obtained  at  any  instant  of  time 
from  the  relevant  equations  of  motion. 
2.  Equations  of  Motion 

The  relationship  of  an  aircraft's  motion  through  a  moving 
airmass  relative  to  earth  leads  to  the  development  of  six  coupled, 
non -linear  equations  [Ref.  12].  Using  the  coordinate  system  in 
Figure  8,  the  following  equations  of  motion  evolve: 

x-Vcosya+Wx  (29) 

n*  Vsin?a+  Wh  (29) 
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^■fX"f'^f  (30) 

V  =  p^cosa-  ^-  gsin?a-  W  x  cosy&-  W  h  sin?a  (32) 

T     eir,  L         S  ^c  Wx     .  W  h  _ 

V?nVSm*  +  W-VC0S*a*— sin*a-  — cos*a  (33) 

3Wx  3Wh 

Note  that     3t     and     3t     were  set  equal  to  zero  due  to  steady  state 

conditions  being  assumed  for  the  windshear  model.    This  assumption 

is  valid  due   to  the  minimal  time   the  airplane  is  exposed   to  the 

microburst  (30  to  60  sec).     Note  that  the  equations  are  for  point 

mass  analysis.    No  dynamic  cross-coupling  is  considered. 

The  equations  of  motion  were  solved  using  a  predictor-corrector 

numerical     scheme.     Euler     First-Forward,     Euler     Half-Step     and 

Richardson  Extrapolation   [Ref.   17]  were  combined   in  the  following 

manner: 

yntl-yn*4tyn  (34) 

ynrf/2     y«     2   yn  (35) 

y^.yO       +41y 

7n+l     'n+1/2     2    'ntl/2  (36) 

yjul'^nil-yml  (37) 
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Second-order  accuracy  is  expected  from  this  numerical  scheme.  By 
providing  input  values  for  aircraft  constants  and  initial  values  for 
the  inertial  reference,  the  equations  of  motion  were  solved  to  obtain 
the  aircraft  performance. 

C.  AIRCRAFT  PERFORMANCE  MODEL 

A  MatLab®  program  was  written  to  take  aircraft  performance 
parameters  and  use  the  set  of  motion  equations  to  measure  aircraft 
response  in  the  windshear.  Initial  values  were  determined  from 
aircraft  performance  parameters  and  the  initial  position.  Iteration  of 
the  set  of  motion  equations  was  based  upon  the  time  interval 
chosen.         Values   for   x,   h,  V,    v,   a,   0,   Yi>  Ya>   Es,   and    Es  were 

tabulated  after  each  iteration.  The  control  of  the  model  was  taken 
from  a  pilot's  point  of  view.  The  controlling  variables  were  8  (deck 
angle  or  pitch  angle)  and  T  (thrust).  Both  variables  can  be  changed 
between  iterations. 

Certain  aircraft  performance  parameters  were  provided  for 
calculating  lift,  drag,  and  AOA.    The  parameters  were  Clo>  C]_c«,  Cd0>  S, 

K,  W,  maximum  thrust  available,  and  maximum  AOA.     Cl  and  Cd 

were  calculated  as  follows: 

Cl"CLo+Cl«*  (38) 

VSo*1^  (39) 
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Note  that  the  lift-curve  slope  was  referenced  from  the  longitudinal 
axis  of  the  aircraft  (see  Figure  8).  Lift  and  drag  were  calculated 
using  the  familiar  relationships: 

LbQSCL  (40) 

D-QSC  ,     v 

D  (41) 

where  the  dynamic  pressure,  Q,  is  based  upon  equivalent  airspeed: 

Q"2poV  (42) 

Initial  values  for  airspeed,  thrust,  and  theta  were  chosen  to 
match  the  aircraft  requirements  for  the  phase  of  flight  under 
scrutiny.  Initial  values  for  x  and  h  were  chosen  referenced  to  the 
microburst  windshear  center.  The  iteration  of  the  model  was  based 
upon  the  time  step,  At. 

Three  aircraft  were  analyzed  for  their  response  in  a  microburst 
windshear.  They  were  the  U.S.  Navy  P-3  (Lockheed  L-188),  the  U.S. 
Navy  T-44  (Beechcraft  King  Air  H-90),  and  a  generic  3-engine  "heavy" 
airline  transport.  The  P-3  is  a  four-engine  turboprop  with  gross 
weights  in  the  medium  range  (75,000  to  135,000  lbs.).  The  T-44  is 
powered  by  two  turboprop  engines  and  falls  within  the  category  of  a 
"light-twin"  transport.  The  3-engine  heavy  airline  transport  is 
powered  by  turbofans  and  as  the  name  implies,  is  capable  of  high 
gross  weights.      Tables    3,    4,    and    5    delineate     the     performance 
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parameters  of  the  three  aircraft.  The  lift  curve  slope  for  the  P-3  is 
greater  than  the  theoretical  prediction.  This  result  is  due  to  the 
influence  of  the  propeller  induced  flow-field.  (Power-off  Clc=5.7  rad-i.) 


TABLE  3: 

P-3  PERFORMANCE  PARAMETERS  [ 

Ref.  18]. 

Landing 

Performance 

Approach 

Takeoff 

On-station 

Parameters 

Configuration 

Configuration 

Configuration 

W  (lbs) 

89,500,  114,000 

89,500;  120,000; 
135,000 

120,000 

S  (ftA2) 

1300 

1300 

1300 

K 

0.05041 

0.05041 

0.05041 

Cdo      4  engine 

0.0567 

0.0551 

0.0213 

3  engine 

0.0630 

- 

- 

CLa  (rad-i) 

6.38 

6.38 

6.25 

Clo 

0.800 

0.800 

0.350 

a  max  (rad) 

0.244 

0.244 

0.209 

TmaxQbf)  4  eng. 

33400 

33400 

33400 

3  eng. 

25050 

- 

- 

q  (rad/s) 

0.0873 

0.0873 

0.0873 

Vref  (ft/s) 

236;  262 

- 

354 

VI  (ft/s) 

- 

204;  214;  229 

- 

V2  (ft/s) 

- 

220;  227;  239 

- 

As  mentioned,  the  AOA  used  in  the  calculations  was  based  upon  the 
longitudinal  axis  of  the  aircraft.  Maximum  thrust  available  and 
gross  weight  were  adjusted  to  meet  the  scenario  requirements. 


34 


TABLE  4. 

3-ENGINE  HEAVY  AIRLINE  TRANSPORT 

PERFORMANCE  PARAMETERS  [Ref. 

10]. 

Landing 

Performance 

Approach 

Takeoff 

Parameters 

Configuration 

Configuration 

W  (lbs) 

362,000 

462,000 

S  (ft~2) 

4578 

4578 

K 

0.059 

0.059 

Cdo 

0.108 

0.098 

CLa  (rad-l) 

4.96 

4.96 

Clq 

0.532 

0.532 

«  max  (rad) 

0.317 

0.314 

TmaxObf) 

126,000 

126,000 

q  (rad/s) 

0.0873 

0.0873 

Vref  (ft/s) 

227 

- 

VI  (ft/s) 

- 

238 

V2  (ft/s) 

- 

255 

TABLE  5:    T-44  PERFORMANCE  PARAMETERS  [Ref.  19]. 


Performance 
Parameters 
W  (lbs) 
S  (ftA2) 
K 

Cdo 

CLa  (rad-l) 

Clo 

a  max  frad) 

TmaxObf) 

q  (rad/s) 

Vref  (ft/s) 

VI  (ft/s) 

V2  (ft/s) 


Landing 

Approach 

Takeoff 

Configuration 

Configuration 

8,280 

8,280 

210 

210 

0.0503 

0.040 

0.120 

0.100 

6.24 

6.24 

0.587 

0.0523 

0.244 

0.227 

3,023 

3,023 

0.0873 

0.0873 

203 

- 

- 

152 

- 

202 
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For    each    Iteration    of    At,    the    following    aircraft,   inertial,    and 
windshear  model  values  were  calculated  and  recorded: 

1)  x  -  Distance  from  microburst  center  in  feet. 

2)  h  -  Absolute  altitude  in  feet. 

3)  V  -  Equivalent  airspeed  in  knots. 

4)  e  (theta)  -  Pitch  angle  in  degrees. 

5)  a  (alpha)  -  Angle-of-attack  in  degrees. 

6)  *ya  -  Flight  path  in  degrees  referenced  to  the  airmass. 

7)  Yj  -  Flight  path  in  degrees  referenced  to  the  earth. 

8)  Vg  -  Ground  speed  in  knots. 

9)  ROC  -  Rate  of  climb  in  ft/sec.  referenced  to  the  earth. 

10)  Thrust  -  instantaneous  thrust  in  lbf . 

11)  Es  -  Specific  energy  in  feet. 

12)  Es  (Esdot)  -  Time  rate  change  of  specific  energy  in  ft/sec. 

13)  Wx  -  Horizontal  wind  speed  in  ft/sec. 

14)  Wh  -  Vertical  wind  speed  in  ft/sec. 

From  the  collection  of  data,  plots  were  made  of  altitude,  theta,  alpha, 
airspeed,  specific  energy,  and  time  with  comparison  to  distance  from 
microburst  center. 
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Four  basic  escape  maneuvers  were  examined  for  critical  phase  of 
flight.  Constant  airspeed,  constant  altitude,  constant  pitch  angle  (0), 
and  constant  angle-of-attack  (a)  escape  maneuvers  were  evaluated 
quantitatively  (engineering  standpoint)  and  qualitatively  (piloting 
standpoint). 

D.  CRITICAL  FLIGHT  PHASES  MODELED 

There  are  three  phases  of  flight  in  which  a  microburst  encounter 
becomes  critical  for  airplane  survival.  Approach  to  landing,  initial 
takeoff  and  climb  out,  and  low-altitude  maneuvering  are  the  three 
phases  of  flight  considered  in  this  analysis.  All  three  phases  occur  at 
an  altitude  of  less  than  1000  feet  AGL.  The  low-altitude 
maneuvering  flight  phase  is  applicable  to  the  P-3  community.  The 
P-3  routinely  operates  200  feet  above  the  water  in  all  weather 
conditions  during  Anti-Submarine  Warfare  (ASW)  missions. 

1.  Landing  Approach 

The  landing  approach  scenario  is  based  upon  a  3  degree 
glideslope  descent  to  landing.  In  this  scenario,  the  microburst  center 
is  placed  10,300  feet  from  the  end  of  runway.  The  simulation  starts 
500   feet   (x=-500   feet)  before   the  microburst   center   on   glideslope 
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(h=566  feet).  This  situation  closely  represents  the  scenario  of  DAL  191. 
The  aircraft  is  exposed  to  the  windshear  immediately  at  start,  t=0. 
The  simulation  is  run  for  a  programmed  length  of  time  or  until 
ground  impact  occurs.  Five  aircraft/weight  combinations  were 
analyzed  for  the  approach  to  landing  scenario.  They  were  a  P-3  at 
89,500  pounds  gross  weight,  a  P-3  at  114,000  pounds  gross  weight,  a  3- 
engine  heavy  airline  transport,  a  T-44,  and  a  P-3  at  89,500  pounds 
gross  weight  with  one  engine  shut  down. 

Starting  aircraft  parameters  were  such  that  the  given  pitch 
angle  (theta)  and  thrust  would  maintain  a  3°  glideslope  at  target 
approach  airspeed  (Vref).  Pitch  angle  and  thrust  were  maintained 
until  a  loss  of  airspeed  equated  to  Vref  minus  20  knots.  At  such 
time,  the  aircraft  performed  one  of  the  following  programmed  escape 
procedures: 

1)  constant  airspeed  -  Maximum  thrust  was  applied  and  a  pitch 
angle  set  to  0°.  This  pitch  angle  was  maintained  until  the 
airspeed  equaled  Vref.  Pitch  angle  was  then  adjusted  to 
maintain  the  airspeed  at  Vref  +  5  knots. 

2)  constant  altitude  -  Maximum  thrust  was  applied  and  a  pitch 
angle  was  set  to  obtain  a  positive  rate  of  climb.  This  pitch 
angle  was  maintained  until  the  target  altitude  (altitude  at 
which  the  maneuver  began)  was  established.  Pitch  angle  was 
then  adjusted  to  maintain  the  target  altitude  by  +  20  feet. 
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3)  constant  thcta  -  Maximum  thrust  was  applied  and  a  pitch 
angle  (theta)  was  set  and  maintained.  Specifically,  theta 
values  of  5°,  8°,  10°,  and  15°  were  used. 

4)  constant  alpha  -  Maximum  thrust  was  applied  and  a  pitch 
angle  was  adjusted  to  get  a  given  angle-of-attack  (alpha).  AOA 
values  of  12,  15,  and  20  units  were  used  for  the  P-3  model. 

The  above  escape  maneuvers  were  constrained  by  certain  limits. 

Pitch  rate  (q)  was  set  at  5°/sec.    Thrust  application  rate  for  the  P-3 

model  was  set  at  0.5  maximum  thrust/sec.     Thrust  application  rate 

for  the  generic  3-engine  heavy  and  T-44  aircraft  was  set  at  0.2 

maximum  thrust/sec.    This  rate  accounts  for  engine  "spool-up"  time. 

Maximum  pitch  angle  was  limited  not  to  exceed  maximum  AOA. 

Appendix  C  contains  the  MatLab®  program  used  for  the  approach  to 

landing  scenario. 

2.  Takeoff 

The    takeoff    scenarios    primarily    explored    the    effects    of 

microburst  windshear  on  the  takeoff  performance  with  penetration 

at  liftoff.    However,  some  analysis  was  performed  on  the  microburst 

center  distance  from  liftoff  point.     Generally,  the  simulation  began 

with  the  aircraft  lifting  off  1200  feet  from  the  microburst  center 

(x=-1200  feet)  at  the  appropriate  liftoff  speed  (Vi).    Initial  pitch  angle, 
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theta,  was  such  to  achieve  takeoff  safety  speed  (V2)  at  50  feet,  no 
windshear.  Maximum  thrust  was  used  for  all  cases.  Execution  of 
the  particular  escape  maneuver  was  begun  when  the  rate  of  climb 
(ROC)  is  less  than  or  equal  to  zero  or  the  airspeed  fell  below  VY  Five 
aircraft/weight  combinations  were  looked  at.  They  were  a  P-3  at 
90,000  pounds,  114,000  pounds,  and  135,000  pounds  gross  weight,  as 
well  as  a  3-engine  heavy  airline  transport  and  a  T-44. 

Four  escape  methods  were  considered  for  a  microburst 
encounter  at  takeoff.  When  the  ROC  was  less  than  or  equal  to  zero, 
or  the  airspeed  fell  below  V2,  one  of  the  following  programs  was 
executed: 

1)  constant  airspeed  -  If  airspeed  was  less  than  V2  at  initiation, 
theta  was  reduced  to  Oo.  if  airspeed  was  greater  than  V2  at 
initiation,  theta  was  increased.  In  both  cases,  theta  was 
manipulated  to  maintain  airspeed  =  V2  ±  5  knots  once  V2  was 
achieved. 

2)  constant  altitude  -  Theta  was  varied  to  maintain  altitude  ±  20 
feet  about  the  target  altitude.  The  target  altitude  was  the 
altitude  at  which  the  escape  maneuver  began. 

3)  constant  theta  -  Theta  was  held  at  a  programmed  constant 
value  throughout  the  maneuver.  Theta  values  used  were  5°, 
10°,  and  15°. 

4)  constant  alpha  -  Theta  was  varied  to  maintain  a  constant 
A0A  value. 


40 


For  the  above  maneuvers,  a  pitch  rate  of  5<>/sec  was  used  and 
maximum  thrust  was  maintained.  Also,  theta  was  reduced  any 
time  critical  AOA  was  exceeded. 

Two  variants  of  the  takeoff  scenario  were  analyzed  for 
performance  sensitivity.  Both  used  the  120,000  pound  P-3.  For  one, 
Vi0f  was  increased  18  knots  to  145  knots.  This  equated  to  an  increase 
of  the  rotate  speed  (Vr)  to  145  knots.  The  second  variant  moved  the 
microburst  center  from  1200  feet  to  2000  feet  from  the  point  of 
liftoff.  This  resulted  in  the  aircraft  gaining  airspeed  and  altitude 
before  encountering  the  severe  horizontal  windshear.  The  second 
scenario  closely  simulated  a  microburst  encounter  after  takeoff. 
MatLab®  programming  for  the  takeoff  scenario  is  listed  in 
Appendix  D. 

3.  P-3  On-Station 

The  on-station  microburst  encounter  scenario  used  a  P-3  at 
120,000  pounds  gross  weight  in  a  4-engine  loiter  configuration.  No 
wing  flaps  or  landing  gear  were  extended.  Only  aircraft  reaction  and 
performance  were  analyzed.    No  specific  escape  maneuver  was  used. 
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This   scenario   tried   to    emulate   an   autopilot   maintaining   altitude 
through  pitch  authority  and  the  flight  crew  controlling  thrust. 

The  microburst  encounter  began  with  the  P-3  at  200  feet  AGL 
and  5000  feet  from  the  microburst  center  (x=-5000  feet).  The  initial 
airspeed  was  210  knots  (prescribed  loiter  airspeed).  Theta  was  varied 
throughout  the  encounter  to  maintain  altitude.  Thrust  was 
maintained  at  loiter  power  until  40  knots  of  airspeed  was  lost.  At 
this  point,  maximum  thrust  was  applied  while  still  maintaining 
altitude.  The  scenario  was  ended  at  5000  feet  on  the  opposite  side 
of  the  microburst  center. 
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HI.    RESULTS  AND  DISCUSSION 

The  microburst  and  inertial  reference  math  models  were 
validated  by  comparing  actual  flight  data  to  model  data.  The 
inertial  model  was  further  tested  for  the  stability  of  the  coupled 
differential  equation  scheme.  Any  weakness  noted  was  considered 
when  the  results  were  analyzed.  Each  critical  phase  of  flight  was 
studied  for  aircraft  response  when  performing  a  particular  escape 
maneuver.  Mathematical  results  were  combined  with  other 
observations  to  build  a  foundation  for  conclusions. 

A.  WINDSHEAR  AND  INERTIAL  REFERENCE  MODEL  VALIDATION 

The  windshear  math  model  was  developed  from  fitting  a 
vortex/source  flowfield  to  recorded  flight  data.  The  wind  field 
recorded  by  the  AA  539  DFDR  (Digital  Flight  Data  Recorder)  showed  a 
definite  vortex  flow.  However,  the  flight  data  from  the  DAL  191  DFDR 
showed  a  different  vortex  flow  arrangement.  In  addition,  a  strong 
outflow  at  the  surface,  and  Increased  vertical  sink,  required  the 
addition  of  source  flows  to  the  windshear  model.  This  use  of  point 
source  flows  led  to  a  much  closer  fit  of  the  original  vortex  model  to 
the  wind  field  DAL  191  experienced.  The  source  flows  have  no  range 
damping  terms.     Therefore,  the  DAL  191  emulation  windshear  model 
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becomes  invalid  at  distances  greater  than  6000  feet  (x=6000  feet) 
from  the  microburst  center.  All  analysis  was  easily  done  within  this 
distance  limit. 

It  is  important  to  note  that  the  aircraft  performance  results  are 
insensitive  to  the  exact  modeling  of  a  particular  windshear.  The 
rapid  change  of  a  headwind  to  a  tailwlnd  is  the  governing  factor  in  a 
microburst.  Secondary  are  the  vertical  down  drafts  that  can  be 
experienced  above  100  feet  AGL  (above  ground  level).  The  windshear 
math  model  contains  both  characteristics  to  the  same  degree  as  the 
windshear  experienced  by  DAL  191. 

Validation  of  the  inertial  reference  model  was  scrutinized  for 
proper  aircraft  response  to  changing  conditions  and  the  effect  of  the 
time  step  size.  Stability  of  the  solutions  obtained  from  the 
differential  equation  numerical  scheme  was  the  greatest  concern. 
Results  generated  from  the  inertial  model  were  studied  for  the  light 
P-3  under  stable  and  turbulent  conditions.  Also,  a  comparison  of  the 
3-engine  transport  model  was  made  with  the  DAL  191  DFDR  data 
during  the  final  seconds  of  the  fateful  flight. 

The  first  approach  was  to  check  the  response  of  a  light  P-3 
initially  stabilized  on  a  3©  descent  path  with  no  windshear.  The  input 
came  from  a  subroutine  that  would  not  vary  pitch  angle  (theta)  or 
thrust.  No  excursions  were  noted  in  descent  path,  airspeed,  or 
alpha,  as  shown  in  Figure  9.    The  time  step  used  was  1  second. 
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The  second  validation  was  to  observe  the  stability  of  the  inertial 
model  solution  with  different  size  time  steps.  A  time  interval  of  0.5 
seconds  was  the  target  time  step  for  running  all  analyses.  Therefore, 
time  steps  of  1,  0.5,  and  0.1  second  were  investigated.  The  light  P-3 
on  approach  to  landing  was  again  used.     This   time,    the    windshear 


-500     0     500  10001500  2000  25003000  3500  4000  4500  5000  5500 
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Figure  9.    Light  P-3  on  a  stabilized  descent  path, 
model  was  incorporated  to  induce  changes  in  the  flight  environment. 
Aircraft   control   was   through   theta   and   thrust   inputs.      Control 
feedback  was  set  to  vary  theta  to  keep  within  ♦  100  feet  of  the 
descent  path  and  to  vary  thrust  to  maintain  Vrcf  within  +  5  knots. 

The  results  are  depicted   in  Figure   10.      Note  the  close  correlation 
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between  the  0.1  and  0.5  second  time  step  cases,  indicating  that  a  0.5 
second  time  step  provides  stability  and  adequate  accuracy  for  the 
desired  results. 


1000        1500        2000        2500        3000 
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Figure  10.     The  effect  of  time  step  size  on  the  inertial 
reference  calculations. 

The  final  validation  was  comparing  the  actual  aircraft  response  to 

model  predictions.    DAL  191  DFDR  data  were  used  to  compare  an  L-1011 

flight  path,  airspeed,  and  alpha  data  to  that  of  the  3-engine  heavy 

airline  transport  model.        Inputs  to  the  aircraft  control  were  the 

recorded  theta  and  thrust  of  DAL  191.    The  windshear  experienced  by 
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DAL  191  was  also  incorporated  for  this  comparison.  Figure  11 
graphically  displays  the  close  characterization  of  the  actual  flight 
path  by  the  model.    Airspeed  and  alpha  do  not  reflect  actual  aircraft 
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Figure  11.  Comparison  of  aircraft  response  in  windshear 
between  the  DAL  191  Flight  L-1011  and  the  3-engine 
transport  model. 

response  as  closely  as  the  P-3  model,  but  do  follow  the  same  general 

response.    The  lag  in  airspeed  and  alpha  is  probably  due  to  the  model 

Cl  and  Cd  equations  being  only  of  second  order. 
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B.  WEAKNESSES  AND  OTHER  INDEPENDENT  VARIABLES 

Some  weaknesses  in  the  math  models  and  unobserved 
independent  variables  exist  in  the  total  analysis.  The  most  notable 
are: 

1)  Cl  and  Cd  value  errors  at  angle  of  attack  near  stall. 

2)  Unobserved  effects  of  dynamic  pressure  and  AOA  upon  engine 
thrust. 

3)  Unobserved     effects     of     rain     and     turbulence     on     aircraft 
aerodynamic  performance. 

4)  Unobserved  effects  of  coupled  longitudinal  and  lateral  dynamic 
modes  excited  by  turbulence  and  pilot  induced  oscillations. 

None    of    the    weaknesses    or    unobserved    variables    is    believed    to 

significantly  impact  the  results  of  the  analysis.    The  overall  concept 

was  to  compare  outcomes  of  different  microburst  escape  strategies 

given  the  same  parameters. 

A  weakness  in  the  aircraft  model  is  the  ability  to  predict  the 

effects  of  flow  separation  at  low  airspeed  and  high  angle  of  attack. 

These   effects   would   be   an   Increase   in   drag   and   a    decrease   in 

generated  lift.    Table  look-up  or  higher  order  equations  could  be  used 

for  the  Cd  and  Cl  expressions.     The  model  does  limit  the  maximum 

theta  not  to  exceed  maximum  angle  of  attack  for  that  aircraft. 
Keeping  this  weakness  in  mind,  the  results  obtained  from  the  math 
models  are  valid  when  determining  the  "best"  escape  maneuver. 
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Engine  thrust  is  somewhat  effected  by  airspeed  and  AOA.  The 
turbojet  and  turbofan-type  engines  can  be  significantly  effected  while 
the  turboprop  is  effected  to  a  much  lesser  degree.  Only  nominal 
values  for  thrust  were  used  in  the  analysis. 

The  effects  of  turbulence  and  rain  on  aircraft  performance  were 
not  programmed  into  the  models.  Studies  by  Wingrove  and  Bach 
[Ref:  2]  analytically  determined  that  the  rain  had  negligible  effect 
upon  DAL  191.  The  NTSB  Report  (National  Transportation  Safety 
Board)  came  to  the  same  conclusion.  What  must  be  considered,  is 
the  effect  of  turbulence  upon  a  pilot's  ability  to  control  the  aircraft. 
The  effect  of  turbulence  on  the  ability  of  a  pilot  to  execute  a 
particular  escape  maneuver  was  kept  in  mind  during  the  writing  of 
the  final  conclusions. 

All  the  aircraft  models  are  limited  to  one  axis  of  freedom.  The 
effects  of  coupled  lateral  and  directional  modes  upon  the  longitudinal 
response  were  not  modeled.  From  past  analysis  and  personal 
aviation  experience,  this  limitation  should  not  significantly  impact 
aircraft  performance.  Again,  aircraft  control  may  become  difficult  if 
an  inherent  mode  is  excited  (such  as  the  phugoid  at  low  airspeed  and 
high  angle  of  attack).  Escape  maneuvers  that  may  produce  such 
dynamic  modes  were   noted. 
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The  effect  of  changing  dynamic  pressure  on  stick  forces  was 
evaluated.  This  is  an  important  consideration  from  the  pilot's 
perspective.  A  loss  of  20-30  knots  airspeed  can  induce  a  significant 
nose-down  pitch  force. 

C.  DATA  REDUCTION  AND  ANALYSIS 

Aircraft  performance  for  three  critical  flight  phases  was  evaluated 
when  exposed  to  a  DAL  191  type  microburst  windshear.  Approach  to 
landing  and  takeoff  flight  phases  were  scrutinized  for  the  best  escape 
maneuver  that  could  be  applied.  The  on-station  flight  phase  was 
examined  for  the  performance  of  a  P-3  while  encountering  the 
microburst  windshear  and  attempting  to  maintain  altitude. 
Calculated  data  were  recorded  and  converted  to  applicable  units  (eg., 
ft/s  to  knots).  The  data  were  then  presented  in  a  graphical  format 
for  analytical  and  subjective  comparison. 

For  all  cases,  tabulated  data  were  converted  to  graphical  form. 
Strip  graphs  for  key  dependent  variables  were  produced  for  each 
aircraft  performing  a  specific  escape  maneuver.  They  are  listed  in 
Appendices  F  and  G.  Combination  graphs  of  flight  path,  airspeed,  and 
specific  energy  compare  escape  maneuvers  for  each  aircraft  in  a 
given  microburst  encounter.  The  abscissa  axis  for  all  the  graphs 
represents   the   distance  from   the   microburst   center   in   feet.      A 
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three-step  approach  was  carried  out  to  deduce  the  "best"  escape 
manuever  for  a  particular  encounter  using  these  graphs. 

First,  the  inertial  reference  frame  of  flight  path  performance  was 
analyzed.  This  involved  using  the  combination  flight  path  graph.  A 
list  is  generated  containing  the  highest  to  lowest  altitude  obtained  at 
a  specific  point  from  the  mlcroburst  center.  For  the  approach  to 
landing  scenario,  only  the  maneuvers  that  resulted  in  flight  paths 
staying  above  the  descent  path  (landing  glideslope)  were  chosen. 
Altitude  values  were  compared  for  each  valid  escape  maneuver  at 
x=5000  feet.  For  the  takeoff  scenario,  only  the  escape  maneuvers 
that  did  not  lead  to  ground  impact  were  listed.  Here,  altitude  values 
were  compared  at  x=4000  feet. 

The  second  step  was  to  cross  reference  each  selected  maneuver's 
airspeed  from  the  combination  airspeed  graph.  This  supplied  an 
insight  of  the  airmass  reference  performance.  The  airspeed  values  at 
the  point  of  Interest,  and  throughout  the  microburst  encounter, 
allowed  a  subjective  analysis  of  the  validity  of  the  results.  Airspeeds 
that  fell  deep  within  the  power-on  stall  region  were  scrutinized  with 
the  strip  graphs  comparing  calculated  AOA  values.  Any  maneuver 
that  resulted  in  a  very  low  airspeed,  high  AOA  condition  was  rejected 
from  the  list. 

The  third  step  evolved  using  the  specific  energy  combination 
graph  for  a  final  resolution.     Specific  energy  incorporates  both  the 
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inertial  and  the  airmass  reference.    The  eligible  maneuvers  that  had 

the  highest  specific  energy  value  at  the  point  of  interest  analytically 

are   the    "best"    escape   maneuvers   for   that   aircraft   configuration, 

given  that  type  of  microburst  encounter. 

1.  Approach  to  Landing  Microburst  Encounter  Analysis 

Figure  12,  Figure  13,  and  Figure  14  depict  the  results  obtained 

for  different  escape  maneuvers  executed  by  a  P-3  at  89,500  pounds 

gross  weight.     Table  6  compares  the  results  observed  at  x=5000  feet. 

The   constant   airspeed   and   12   unit   AOA   escape   maneuvers  were 

rejected  initially  due  to  resulting  flight  trajectories  below  the  descent 

path  (3°  landing  approach  glideslope)  as  depicted  in  Figure  12.      The  20 

unit  AOA  escape  maneuver  was  subsequently  rejected  because  of  low 

airspeed,  high  AOA  observed  in  Figure  13.    Figure  14  shows  that  the 

5o  theta  escape  maneuver  results  in  the  highest  specific  energy  value 

of  the  remaining  list. 

TABLE  6.    RELATIVE  VALUE  RANKING  FOR  A  P-3  AT  89,5000LBS. 

maneuver  altitude                    airspeed             specific  energy 

20  unit  AOA  1                        rejected                    

const,  altitude  2                              5                              5 

15o  theta  344 

15  unit  AOA  4                                3                                3 

lOo  theta  5                              2                              2 

5o  theta  6                              1                              1 

Figure  15,  Figure  16,  and  Figure  17  depict  the  results  obtained 
for  different  escape  maneuvers  executed  by  a  P-3  at  114,000  pounds 
gross  weight.    Table  7  compares  the  results  observed  at  x=5000  feet. 


52 


< 
o 
< 


c 
o 


£  f  <-J     ■+• 

^   ■*■'     w 


o  to  10     O 


4' 


r 

* 
4- 

< 

+* 

Cl 

O 

W 

< 

-C 

** 

<v 

VI 

♦* 

C 

rr 

* 

♦' 

c 

<i> 

o 

10 

3 

TJ 

0> 

c 
o 

(N 

in 

■o 

o 

«— 

Figure  12.  Flight  path  data  comparison  of  a  P3  at 
89,500105  performing  different  escape  maneuvers  during 
an  approach  to  landing  microburst  encounter. 
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Figure  13.  Airspeed  data  comparison  of  a  P3  at  89,5001bs 
performing  different  escape  maneuvers  during  an 
approach  to  landing  microburst  encounter. 
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Figure  14.  Specific  energy  data  comparison  of  a  P3  at 
89, 500 lbs  performing  different  escape  maneuvers  during 
an  approach  to  landing  microburst  encounter. 


55 


The  constant  airspeed,  12  unit  AOA,  and  5°  theta  escape  maneuvers 

were   rejected   initially    due   resulting   flight   trajectories   below   the 

descent  path  as  depicted  in  Figure  15.     The  20  unit  AOA  and  constant 

altitude  escape  maneuvers  were  subsequently  rejected  because  of  low 

airspeeds,  high  AOAs  observed  in  Figure  16.    Figure  17  shows  that  the 

15  unit  AOA  escape  maneuver  results  in  the  highest  specific  energy 

value  of  the  remaining  list. 

TABLE  7.  RELATIVE  VALUE  RANKING  FOR  A  P-3  AT  114,000LBS. 

maneuver  altitude  airspeed  specific  energy 

20  unit  AOA  1  rejected 

const,  altitude  2  rejected                     

15o  theta  3                                3                               3 

10o  theta  4                               2                              2 

15  unit  AOA  5                               1                               1 

Figure  18,  Figure  19,  and  Figure  20  depict  the  results  obtained 

for  different  escape  maneuvers  executed  by  a  3-englne  heavy  airline 

transport.     Table  8  compares  the  results  observed  at  x=5000  feet. 

The  constant    airspeed,  5©    theta,    and    10o    theta    escape  maneuvers 

TABLE  8.    RELATIVE  VALUE  RANKING  FOR  A  3-ENGINE  AIRLINE 
TRANSPORT. 
maneuver  altitude  airspeed  specific  energy 

const,  altitude  1  2  2 

15o  theta  2  1  1 

were   rejected   initially   due   resulting   flight   trajectories   below   the 

descent  path  as  depicted  in  Figure  18.      Figure  20  shows  that  the  15o 

theta  escape  maneuver  results  in  the  highest  specific  energy  value. 
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Figure  15.  Flight  path  data  comparison  of  a  P3  at 
114,0001bs  performing  different  escape  maneuvers  during 
an  approach  to  landing  microburst  encounter. 
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Figure  16.  Airspeed  data  comparison  of  a  P3  at  114,0001bs 
performing  different  escape  maneuvers  during  an 
approach  to  landing  microburst  encounter. 
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Figure  17.  Specific  energy  data  comparison  of  a  P3  at 
114,0001bs  performing  different  escape  maneuvers  during 
an  approach  to  landing  microburst  encounter. 
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Figure  18.  Altitude  data  comparison  of  a  3-engine  heavy 
airline  transport  performing  different  escape  maneuvers 
during  an  approach  to  landing  microburst  encounter. 
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Figure  19.  Airspeed  data  comparison  of  a  3 -engine  heavy 
airline  transport  performing  different  escape  maneuvers 
during  an  approach  to  landing  microburst  encounter. 
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Figure  20.  Specific  energy  data  comparison  of  a  3-engine 
heavy  airline  transport  performing  different  escape 
maneuvers  during  an  approach  to  landing  microburst 
encounter. 
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Figure  21,  Figure  22,  and  Figure  23  depict  the  results  obtained 

for    different    escape    maneuvers    executed    by    a    T-44.       Table    9 

compares    the    results    observed    at    x=5000    feet.  The    constant 

airspeed,   10   unit   AOA,   and    16    unit   AOA   escape   maneuvers   were 

rejected  initially  due  to  resulting  flight  trajectories  below  the  descent 

path  as  depicted  In  Figure  21.    The  25  unit  AOA  and  constant  altitude 

escape  maneuvers  was  subsequently  rejected  because  of  low  airspeeds 

and  high  AOAs  observed  in  Figure  22.      Figure  23  shows  that  the  5° 

theta  escape  maneuver  results  in  the  highest  specific  energy  value  of 

the  remaining  list. 

TABLE  9.    RELATIVE  VALUE  RANKING  FOR  A  T-44. 

maneuver  altitude  airspeed  specific  energy 

const,  altitude  - 1  rejected  

15o  theta  2  3  3 

25  unit  AOA  3  rejected  

10o  theta  4  2  2 

5o  theta  5  1  1 

Figure  24,  Figure  25,  and  Figure  26  depict  the  results  obtained 

for  different  escape  maneuvers  performed  by  a  P-3  at  89,500  pounds 

gross  weight  operating  on  3  engines.     Table  10  compares  the  results 

observed  at  x=5000  feet.    The  constant  airspeed,  12  unit  AOA,  and  5o 

theta  escape  maneuvers  were  rejected  initially  due  to  resulting  flight 

trajectories  below  the  descent  path  as  depicted  in  Figure  24.    The  20 

unit  AOA  and  constant  altitude  escape  maneuvers  were  subsequently 

rejected    because    of    low    airspeeds    and    high    AOAs    observed    in 

Figure  25.    The  15o  theta  escape  maneuver  was  also  rejected  for  the 
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Figure     21.  Altitude     data     comparison     of     a     T-44 

performing     different     escape     maneuvers     during     an 
approach  to  landing  microburst  encounter. 
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Figure  22.  Airspeed  data  comparison  of  a  T-44 
performing  different  escape  maneuvers  during  an 
approach  to  landing  microburst  encounter. 
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Figure  23.  Specific  energy  data  comparison  of  a  T-44 
performing  different  escape  maneuvers  during  an 
approach  to  landing  microburst  encounter. 
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Figure  24.  Altitude  data  comparison  of  a  P-3  at  89,5001bs 
operating  on  3  engines  performing  different  escape 
maneuvers  during  an  approach  to  landing  microburst 
encounter. 
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Figure  25.  Airspeed  data  comparison  of  a  P-3  at 
89,5001bs  operating  on  3  engines  performing  different 
escape  maneuvers  during  an  approach  to  landing 
microburst  encounter. 
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Figure  26.  Specific  energy  data  comparison  of  a  P-3  at 
89,5001bs  operating  on  3  engines  performing  different 
escape  maneuvers  during  an  approach  to  landing 
microburst  encounter. 
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same    reasons.       Although    the    airspeed    is    only    within    the    flow 

separation  region  (steady  state  stall  buffet),  analysis  of  the  AOA  (see 

Figure  F34)  shows  that   critical   alpha  was  sustained   on  the  latter 

parts  of  the  maneuver.    Figure  26  shows  that  the  15  unit  AOA  escape 

maneuver    results    in    the    highest    specific    energy    value    of    the 

remaining  list. 

TABLE  10.    RELATIVE  VALUE  RANKING  FOR  A  P-3  AT  89,5000LBS 

OPERATING  ON  3  ENGINES. 

maneuver                 altitude                   airspeed  specific  energy 

const,  altitude                     1                          rejected  

15°  theta                             2                         rejected  

20  unit  AOA                      3                        rejected  

10o  theta                              4                               2  2 

15  unit  AOA                         5                                1  1 

2.  Takeoff  Mlcroburst  Encounter  Analysis 

Figure  27,  Figure  28,  and  Figure  29  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  P-3  at  90,0001bs  gross 
weight.  Altitude,  airspeed,  and  specific  energy  values  are  compared 
at  x=4000ft.  All  six  escape  maneuvers  resulted  in  ground  plane 
clearance  as  depicted  In  Figure  27.  The  15  unit  AOA  escape  maneuver 
was  subsequently  rejected  because  of  low  airspeed,  high  AOA  as 
observed  in  Figure  28.  Figure  29  shows  that  all  escape  maneuvers, 
except  the  15  unit  AOA  escape  maneuver,  resulted  in  a  grouped  single 
value.  This  shows  that  all  relevant  escape  maneuvers  provide  the 
same  final  specific  energy. 
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Figure  27.  Flight  path  data  comparison  of  a  P-3  at 
90,0 00 lbs  performing  different  escape  maneuvers  during 
a  takeoff  microburst  encounter. 
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Figure  28.  Airspeed  comparison  of  a  P-3  at 
9  0,000  lbs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 


72 


T? 

<!■ 

■1- 

1' 

x> 

I-! 

u 

V) 

•+*■ 

i- 

^- 

•r- 

<v 

«> 

•+s 

•+s 

Vi 

11 

c 

c 

o 

o 

o 

o 

*  r  9  9 

a.  -C 

-  o>  *:  ^ 

r-r.  a.  c  c 

"°  O  CM  IT) 

ID  —  —  — 


.1   <1   i 


O 
O 

o 
Tf 

o 
o 
id 
ro 

o 
o 
o 

o 
o 
id 

CM 

o 
o 
o 

CM 

o 
o 

ID 


o 
o 
o 


o 
o 

ID 


-  o 


o 
o 

ID 
I 

o 
o 
o 


o 
o 

ID 


o 
o 

CM 


O 

o 


o 
o 
o 


o 
o 


O 
O 
00 


o 
o 


O  I 

o 


Figure  29.  Specific  energy  comparison  of  a  P-3  at 
90,0001bs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  30,  Figure  31,  and  Figure  32  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  P-3  at  120,000  pounds 
gross  weight.  Value  comparison  is  made  at  x=4000  feet.  The 
constant  airspeed,  constant  altitude,  12  unit  AOA,  15  unit  AOA,  and  5° 
theta  escape  maneuvers  were  rejected  initially  due  ground  impact  as 
depicted  in  Figure  30.  The  15°  theta  escape  maneuver  was 
subsequently  rejected  because  of  low  airspeed  and  high  AOA  as 
observed  in  Figure  31.  The  only  viable,  however  marginal, 
performance  observed  was  for  the  10°  theta  escape  maneuver. 
Figure  31  shows  that  the  10°  theta  maneuver  resulted  in  an  airspeed 
drop  of  30   knots  below  V2   at  one  point  followed  by  a   sustained 

airspeed  of  20  knots  below  VY 

Included  with  the  120,000  pound  P-3  data  is  the  resulting  flight 
performance  when  Vr  was  increased.     A  large  gain  in  altitude  and 

specific  energy  over  other  performance  profiles  is  observed. 
Increasing  the  rotate  speed  by  18  knots  coupled  with  flying  a  constant 
10°  theta  provides  a  2  fold  increase  in  altitude  and  a  27?  increase  in 
specific  energy.  This  Increase  In  specific  energy  Is  important  when 
compared  to  the  spread  of  specific  energy  values  of  different  escape 
maneuvers. 
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Figure  30.  Flight  path  comparison  of  a  P-3  at 
12 0,000 lbs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure     31.  Airspeed     comparison     of     a     P-3     at 

120, 000 lbs    performing    different    escape    maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  32.  Specific  energy  comparison  of  a  P-3  at 
12 0,000 lbs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  33,  Figure  34,  and  Figure  35  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  P-3  at  135,0001bs 
gross  weight.  It  can  be  graphically  observed  that  a  P-3  loaded  to 
maximum  takeoff  weight  is  at  the  mercy  of  a  strong  microburst 
encountered  at  takeoff.  Although  the  15°  theta  maneuver  misses  the 
ground,  the  low  airspeed  and  high  AOA  experienced  disqualifies  it  as 
viable. 

Figure  36,  Figure  37,  and  Figure  38  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  T-44.  Altitude, 
airspeed,  and  specific  energy  values  are  compared  at  x=4000  feet.  All 
five  escape  maneuvers  resulted  in  ground  plane  clearance  as  depicted 
in  Figure  36.  All  five  escape  maneuvers  provided  adequate  margin  of 
airspeed  and  AOA  as  depicted  in  Figure  37.  Note  that  power-on  stall 
speed  for  the  T-44  Is  less  than  90  knots.  Figure  38  shows  that  all 
escape  maneuvers  resulted  in  a  grouped  single  value  of  specific 
energy.  As  seen  above,  all  escape  maneuvers  for  a  given  aircraft 
provide  the  same  final  specific  energy  during  a  takeoff  microburst 
encounter.    The  only  exception  is  when  Vr  is  increased. 

The  results  of  one  variant  to  the  takeoff  encounter  are 
depicted  in  Figure  39,  Figure  40,  and  Figure  41.  Here,  the  microburst 
center  is  moved  an  added  800  feet  from  the  liftoff  point.  This  results 
in  a  delayed  encounter  with  the  severest  part  of  the  windshear.  A 
P-3  at  120,00  pounds  gross  weight  was  used  to  perform  the  different 
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Figure  33.  Flight  path  comparison  of  a  P-3  at 
135,0001135  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  34.  Airspeed  comparison  of  a  P-3  at 
13 5,000 lbs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  35.  Specific  energy  comparison  of  a  P-3  at 
13 5,0 00 lbs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  36.  Flight  path  comparison  of  a  T-44 
performing  different  escape  maneuvers  during  a 
takeoff  microburst  encounter. 
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Figure      37.  Airspeed      comparison      of     a      T-44 

performing    different    escape    maneuvers    during    a 
takeoff  microburst  encounter. 
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Figure  38.  Specific  energy  comparison  of  a  T-44 
performing  different  escape  maneuvers  during  a 
takeoff  microburst  encounter. 
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Figure  39.  Flight  path  comparison  of  a  P-3  at 
120,0001bs  performing  different  escape  maneuvers 
during  an  after  takeoff  microburst  encounter. 
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Figure  40.  Airspeed  comparison  of  a  P-3  at 
1 20,00 Olbs  performing  different  escape  maneuvers 
during  an  after  takeoff  microburst  encounter. 
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Figure  41.  Specific  energy  comparison  of  a  P-3  at 
120,0001bs  performing  different  escape  maneuvers 
during  an  after  takeoff  microburst  encounter. 
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escape  maneuvers  including  the  increased  Vr  technique.      Table  11 

compares  the  results  observed  at  x=4000  feet.  The  12  unit  AOA 
escape  maneuver  was  rejected  initially  due  to  ground  impact  as 
depicted  in  Figure  39.  The  15  unit  AOA  and  15°  theta  escape 
maneuvers  were  subsequently  rejected  because  of  low  airspeeds  and 
high  AOAs  observed  in  Figure  40.  Figure  41  shows  that  the  close 
grouping  of  the  specific  energy  values  still  exists  for  the 
non-increased  Vr  maneuvers.    Again,  the  increased  Vr  resulted  in  a 

significant  improvement  in  altitude  and  specific  energy. 

TABLE  11.    RELATIVE  VALUE  RANKING  FOR  A  P-3  AT  120,000LBS 
WITH  AN  AFTER  TAKEOFF  ENCOUNTER, 

maneuver                 altitude                   airspeed  specific  energy 

increased  Vr                      1                              2  1 

15°  theta                             2                        rejected  

10o  theta                              3                               4  2 

const,  altitude                    4                              3  3 

15  unit  AOA                        5                        rejected  

const,  airspeed                  6                              1  4 

3.  On -Station  Encounter  Analysis 

The  analysis  showed  that  a  P-3  at  120,0001bs  gross  weight 
successfully  navigated  a  strong  microburst  during  loiter  operations. 
Figure  42  graphically  depicts  the  results.  Altitude  deviation  was  no 
greater  than  ±  20ft.  Theta  input  never  exceeded  10°.  AOA  remained 
below  critical  angle.  Additional  thrust  was  not  applied  until  40  knots 
of  airspeed  was  lost  in  the  encounter  (a  drop  from  210  to  170  knots). 
Power  was  added  at  x=2700  feet  and  a  time  lapse  of  23sec  from  the 
initial  point  of  the  encounter. 
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Figure  42.    P-3  at  120,0001bs  gross  weight  encountering  a 
microburst  windshear  while  on-station  loiter. 
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D.  OBSERVATIONS 

The  results  obtained  from  the  analytical  analysis  for  each 
encounter  are  integrated  with  other  relative  information  to  answer 
the  three  posed  questions: 

1)  What  is  the  optimum  microburst  escape  procedure  given  the 
type  of  encounter  and  the  available  flight  performance 
information? 

2)  Is  the  optimum  escape  procedure  effected  by  change  of  gross 
weight  or  available  thrust? 

3)  What  flight  instrument  indications,  flight  control  feeling,  and 
dynamic  response  would  be  expected  during  the  optimum 
escape  procedure? 

The  answers  to  these  questions  are  directed  to  the  P-3.     However, 

comparison  is  made  to  the  T-44  and  3-engine  heavy  airline  transport 

where  appropriate.  Other  important  issues  must  also  be  addressed 

for  comparison  of  results  to  other  wlndshear  studies.      These  Issues 

include  the  impact  of  wing  loading  (W/S),  thrust  to  weight  (T/W),  and 

early  liftoff  speeds  for  the  P-3. 

1.   Available  Aircraft  Flight  Instruments 

The  available  flight  information  and  type  of  presentation  are 

extremely  influential  on  the  choice  of  the  optimum  escape  procedure. 

Both  the  P-3  and  T-44  have  "conventional"  flight  instruments  [Ref:  8 

and  Ref:  9].     Conventional  implies  gyroscopic-stabilized  attitude  and 

heading  indication  in  combination  with  pitot-static  airspeed,  altitude, 

and  vertical  speed  indication  (VSI).    Angle  of  Attack  (AOA)  indication 
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is  installed  on  the  P-3  and  7-44.  However,  the  AOA  is  heavily 
damped  and  is  primarily  designed  for  steady-state  flight  conditions. 
Standard  radio  navigation  equipment  is  installed  on  both  aircraft 
encompassing  an  Integrated  Landing  System  (ILS),  VHF  Omni  Radio 
(VOR),  and  Tactical  Navigation  (TACAN).  The  P-3  in  addition  has 
Inertial  Navigation  Systems  (INS),  but  no  Information  is  presented  at 
the  flight  station.  Both  aircraft  incorporate  autopilots  and  flight 
directors.  The  P-3  autopilot  is  controlled  by  Control  Wheel  Steering 
(CWS)  input  only.  The  T-44  autopilot  is  controlled  through  thumb 
wheels  and  can  couple  to  the  flight  director.  Flight  directors  on  both 
aircraft  can  provide  navigational  steering.  However,  the  primary 
design  of  the  P-3  flight  director  is  tactical  while  the  primary  design 
of  the  T-44  flight  director  is  navigational.  In  summary,  both  the  P-3 
and  T-44  have  Attitude  Indicators  (AI),  damped  AOA  indication,  and 
pitot-static  instruments  for  reference  while  executing  an  escape 
maneuver. 

2.  Optimal  Escape  Procedure 

Flying  a  constant  theta  with  reference  to  the  AI  seems  to  be 
the  optimal  escape  procedure  during  a  microburst  encounter.  This 
conclusion  is  supported  analytically  and  qualitatively  for  the  P-3  and 
T-44.  Analytically,  superior  performance  was  obtained  in  a  few 
circumstances  from  constant  AOA  maneuvers  with  an  alpha  value 
between    approach    and    stall.       However,    qualitatively    the    AOA 
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indication  in  both  aircraft  does  not  lend  itself  as  a  viable  reference 
for  two  reasons.  First,  the  AOA  indication  may  be  significantly  in 
error  because  of  instrument  lag  in  turbulent  conditions.  Second,  P-3 
and  T-44  pilots  seldom  use  AOA  indication  during  the  approach  phase. 
The  transition  from  attitude/airspeed  to  attitude/ AOA  reference 
during  an  intense  situation  as  a  microburst  encounter,  is  more  than 
can  be  expected  from  a  pilotage  standpoint.  The  only  viable 
alternative  to  flying  constant  AOA  was  a  constant  theta  maneuver. 
Note  that  quantitatively  in  all  circumstances,  a  constant  theta 
maneuver  was  superior  to  or  closely  matched  with  any  other 
attempted  escape  maneuver. 

Theta  values  ranging  from  5o  to  15°  were  identified  as  optimal 
for  an  approach  to  landing  encounter.  The  theta  value  of  10°,  if  not 
optimal,  provided  suitable  recovery  for  the  P-3  In  the  configurations 
analyzed.  The  light,  the  heavy,  and  the  3-engine  P-3  successfully 
navigated  the  microburst  encounter  utilizing  the  10°  theta  escape. 
This  result  is  important  from  the  standpoint  that  10°  is  easy  to 
remember  and  easy  to  read  on  the  AI.  Also,  10°  theta  worked  well 
for  the  T-44  in  the  approach  configuration  (approach  flap  setting).  A 
theta  of  15°  was  optimum  for  the  T-44  in  the  takeoff  configuration 
(flaps  up).  This  is  important  because  future  P-3  pilots  receive  their 
first  training  In  T44s. 
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3.  The  Effect  Of  Specific  Energy  On  Survival  Probability 

Escaping  from  a  takeoff  microburst  encounter  has  a  lower 
probability  for  success  compared  to  the  approach  to  landing 
encounter.  The  type  of  escape  maneuver  selected  in  a  takeoff 
encounter  has  a  lesser  effect  on  the  outcome.  This  result  is  caused 
by  the  excess  airspeed  and  additional  altitude  available  during  the 
approach  phase.  (Note  that  the  P-3  community  flies  an  unusually 
fast  approach  speed  compared  to  their  civilian  counterparts.  The  P-3 
approach  Vrcf  airspeed  is  determined  from  1.35V6+5  knots  [Ref.  20], 

which  equates  to  1.4VS  for  a  gross  weight  of  114,000  pounds.     The 

industry  standard  for  Vrcf  is  1.3Vs  +  an  added  factor  as  dictated  by 

the  type  of  aircraft.)  Specific  energy  can  be  used  for  an  interesting 
comparison  between  an  approach  to  landing  and  a  takeoff  microburst 
encounter.  A  P-3  weighing  89,000  pounds  has  specific  energy  values 
of  1230  feet  and  1040  feet  at  x=500  feet  and  x=4000  feet  respectively 
during  an  approach  to  landing  encounter  (reference  Figure  14). 
During  a  takeoff  encounter,  a  P-3  at  90,000  pounds  has  a  specific 
energy  values  of  810  feet  and  980  feet  at  the  same  x  distances  from 
the  microburst  center  (reference  Figure  29).  Both  executed  a  10° 
theta  escape  maneuver  and  both  successfully  navigated  the 
microburst.  Note  that  the  approach  to  landing  phase  lost  energy 
while  the  takeoff  phase  gained  energy. 
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An  important  aspect  of  specific  energy  comparison  is  seen  in 
the  increased  rotate  speed  data.  Figure  41  graphically  depicts  a 
120,000-pound  P-3  starting  with  an  additional  specific  energy  of  200 
feet  compared  to  the  normal  rotate  speed  profiles.  The  exit  specific 
energy  was  175  feet  greater  than  the  comparable  10°  theta 
maneuver.  The  175  feet  were  translated  almost  completely  to 
altitude  (reference  Figure  30  and  Figure  31).  Data  support  that  an 
increased  rotate  speed  is  far  more  beneficial  than  any  particular 
escape  maneuver  during  a  takeoff  microburst  encounter. 
4.  Stick  Force  vs.  Off -Trim  Airspeed  For  The  P-3 

The  expected  change  in  flight  control  "feel"  for  the  P-3  is  light 
compared  to  large  aircraft.  Flight  tests  show  that  stick  force  versus 
change  from  trim  airspeed  results  in  a  shallow  gradient  for  this  size 
aircraft  [Ref.  18:p.  37].  In  the  takeoff/approach  flap  configuration  (18°), 
with  an  aft  C.G.  of  29*  M.A.C.,  a  stick  force  of  an  11-pound  pull  is 
needed  with  a  30-knot  decrease  from  trim  airspeed.  Stick  force 
increases  to  a  16  pound  pull  with  the  most  forward  C.G.  of  16?  M.A.C. 
Flight  test  data  [Ref.  18]  also  show  that  the  stick  force  to  trim 
airspeed  gradient  is  not  appreciably  effected  by  landing  gear  position 
or  gross  weight.  The  same  flap  setting  (takeoff /approach)  is  used  for 
landing  approach  and  takeoff.  Thus,  the  elevator  "feel"  is  expected  to 
be  the  same  for  either  approach  to  landing  or  takeoff  microburst 
encounters. 
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The  P-3  Flight  Simulator  (Device  2F87F)  was  used  to  confirm 
the  flight  test  data  relating  to  off-speed  to  stick  force  and  landing 
gear  drag.  The  simulator  allows  one  or  more  flight  parameters  to  be 
frozen.  This  feature  is  used  by  freezing  the  altitude  to  500  feet  and 
trimming  the  aircraft  with  maximum  power  at  reference  airspeed. 
The  airspeed  is  then  reset  and  frozen  to  the  low  airspeed  value 
expected  during  a  microburst  encounter.  The  control  force  required 
to  maintain  a  pitch  attitude  was  then  directly  measured  and  printed 
out.  The  P-3  flight  simulator  data  (Appendix  H)  showed  that  the 
elevator  pull  needed  at  Vtrim  minus  30kts  decreased  by  5  pounds 

when  full  power  was  applied.  Level  acceleration  maneuvers  were 
performed  with  the  landing  gear  extended,  retracted,  and  during 
retraction.  Simulator  data  showed  that  the  landing  gear  retraction 
cycle  does  not  increase  drag. 

5.  Weight,  Wing  Loading,  and  Thrust  to  Weight  Effects 
It  is  obvious  from  the  analysis  that  a  heavy  P-3  is  effected  by 
microburst  windshear  to  a  higher  degree  than  a  light  P-3.  The  T-44 
was  relatively  unaffected  by  the  windshear.  The  question  arises 
whether  weight,  wing  loading,  or  thrust  to  weight  has  the  largest 
effect  on  microburst  survival.  Table  12  compares  the  weight  (Wt), 
wing  loading  (W/S),  thrust  to  weight  (T/W),  and  specific  energy  values 
for  three  different  aircraft  during  an  approach  to  landing  microburst 
encounter.    The  specific  energy  loss  (Es  loss)  column  shows  the  energy 
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(in  feet)  lost  during  the  encounter.     The  specific  energy  spread  (E6 

spread)  column  shows  the  spread   (in  feet)  of   the  different  escape 

techniques  measured     after  the  encounter.     Apparently,  the  higher 

the  thrust  to  weight  ratio,  the  less  energy  is  lost.    Light  wing  loading 

equates  to  an  individual  escape  procedure  having  less  effect  upon  the 

result.      Weight,  as  a  separate  parameter,  is  irrelevant.      In  other 

words,    increased    thrust    to    weight    decreases    the    effect    of    the 

windshear  upon  the  aircraft  and  light  wing  loading  allows  increased 

efficiency  of  energy  transformation. 

TABLE  12.     WEIGHT,  WING  LOADING,  AND  THRUST  TO  WEIGHT 
EFFECTS  ON  Ps  VALUES. 
aircraft         Wt  (lbs)  W/£  WtL 

P-3  89,500  68.85  0.37 

P-3  114,000  88.08  0.29 

T-44  8,280  39.43  0.37 

3-eng.  362,000  79.07  0.35 

hvy. 

6.  Early  Liftoff  Speed  For  The  P-3 

The  analyses  did  not  look  at  the  effect  of  a  microburst 
windshear  encountered  during  a  takeoff  roll.  The  FA  A  Windshear 
Training  Aid  [Ref.  14]  points  out  that  If  windshear  Is  encountered  past 
V~i,  the  aircraft  should  be  committed  to  flight  no  later  than  with 

2000  feet  of  runway  remaining.    Although  the  airspeed  may  fall  and 
stay  below  Vr,     most  aircraft  can  go  airborne.     In  this  regard,  a 

theoretical  liftoff  speed  was  calculated  for  the  P-3.     Equations  38,  40, 
and  42  were  used  In  conjunction  to  solve  for  V.    The  C\a  data  were 
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s  loss 

Ec  spread 

380 

145 

455 

155 

285 

50 

345 

100 

extrapolated  from  flight  test  data  [Ref.  18].     Table  13  compares  the 
recommended  Vr  speed  [Ref.  20]  to  the  theoretical  liftoff  speed  for  a 

rotate  pitch  attitude  of  10o.     Ground  effect  was  not  considered  in  the 

calculations. 

TABLE  13.    THEORETICAL  LIFTOFF  SPEEDS 

weight  (lbs)  Vr  (Kts)  theoretical  liftoff  (kts) 

80,000  115  90 

90,000  115  95 

100,000  115  101 

110,000  117  105 

120,000  123  HO 

130,000  129  115 

140,000  134  119 

8.  Summary 

The  original  three  questions  are  answered  for  the  P-3: 

1)  The  optimum  approach  to  landing  microburst  escape 
procedure  is  to  set  and  maintain  maximum  power  while 
simultaneously  setting  a  10o  pitch  attitude.  For  the  takeoff 
microburst  encounter,  increase  rotate  speed  to  140  knots 
then  pitch  to  and  maintain  10°. 

2)  The  optimum  escape  procedures  remain  the  same  for  all 
gross  weights.  The  approach  to  landing  encounter  escape 
procedure  is  the  same  for  four  or  three  operating  engines. 

3)  The  airspeed  will  rapidly  decay  and  remain  abnormally  low 
during  a  microburst  penetration.  An  elevator  force  of  5  to 
10  pounds  can  be  expected  to  maintain  a  10°  pitch  attitude 
with  full  power.  Flight  tests  show  that  the  P-3  has  no 
unusual  short  period,  phugoid,  nor  cross-couple  dynamics  to 
contend  with. 
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The  answers  for  the  T-44  should  be  similar  to  the  P-3.  The  10°  pitch 
attitude  was  optimal  for  the  approach  to  landing  encounter.  A 
pitch  attitude  of  15°  is  optimal  for  takeoff,  the  difference  being  that 
the  flaps  are  up  in  the  takeoff  configuration.  Not  enough 
performance  data  were  available  to  predict  the  expected  stick  forces. 
However,  personal  flight  experience  has  shown  that  the  T-44  exhibits 
longitudinal  forces  and  responses  very  similar  to  the  P-3. 
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IV.    CONCLUSIONS 

The  aim  of  this  thesis  was  to  analytically  produce  microburst 
escape  flight  procedures.  Study  was  directed  at  two  turboprop 
aircraft  the  Navy  routinely  operates.  The  P-3  and  T-44  were  closely 
scrutinized  for  the  optimal  escape  procedure  if  confronted  with  a 
microburst  penetration.  The  microburst  modeled  was  patterned 
after  that  encountered  by  a  Delta  Airlines  L-1011  at  Dallas/Ft.  Worth 
on  August  2,  1985.  The  conclusions  drawn  are  presented  within  the 
format  of  the  FAA  Windshear  Training  Aid  [Ref.  14].  Viable  flight 
procedures  and  associated  precautions  are  presented  for  both  aircraft. 
The  precautions  should  be  adhered  to  if  an  encounter  with  a 
microburst  windshear  is  a  possible  expectation.  NATOPS  change 
recommendations  reflecting  these  conclusions  will  be  forwarded  to  the 
appropriate  aircraft  model  managers. 

A.  MICROBURST  ESCAPE  PROCEDURES  FOR  THE  P-3 

A  target  pitch  attitude  of  10°  is  optimal  for  all  weights  and  with 
four  or  three  engines  operating.  This  target  pitch  attitude  is  used  for 
approach  to  landing  and  takeoff  microburst  encounters  with  the  flaps 
in  the  Takeoff /Approach  position.  The  airspeed  at  which  the  P-3  lifts 
off  is  the  critical  factor  in  the  takeoff  scenario. 
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1.  Approach  To  Landing  Microburst  Encounter 

Precautions: 

1)  Avoid  thunderstorm  conditions.  Delay  the  approach  if 
possible. 

2)  Use  Takeoff/Approach  flap  setting.  Do  not  extend  flaps 
to  Land  position  until  the  runway  is  made. 

3)  Use  a  precision  approach  procedure. 

4)  The  approach  speed,  Vref,  should  be  equal  to  1.35VV5 
knots.  Additional  airspeed  is  not  normally  warranted. 
This  Vref  will  provide  a  comfortable  margin  above  stall 
speed.  Note  that  additional  airspeed  increases  landing 
distance  (5  knots  faster,  10$  farther). 

5)  Determine  expected  descent  rate. 

6)  Consider  using  the  autopilot  (ASW-31)  with  the  altitude 
hold  switch  in  the  off  position.  The  autopilot  will 
provide  wing  rock  damping  and  pitch  attitude 
augmentation  in  the  face  of  turbulence.  It  will  also 
provide  pitch  attitude  hold  if  a  microburst  escape 
maneuver  is  warranted.  Insure  autopilot  disconnect 
passing  200  feet  AGL. 

7)  Attain  a  stabilized  airspeed  approach  before  passing 
1000  feet  AGL.  Minimize  power  lever  movement 
beyond  this  point.  Maintain  glideslope  with  pitch.  The 
airspeed  indicator  will  serve  as  a  windshear  indicator. 

8)  Strong  consideration  should  be  given  to  executing  the 
microburst  escape  procedure  if  one  of  the  indications  is 
observed: 

a)  a  rapid  and  sustained  airspeed  loss  of  20  knots  below 
Vrcf; 

b)  a  descent  rate  500  feet  per  minute  greater  than  the 
predetermined  value; 

c)  greater  than  1  dot  low  on  the  ILS  glideslope  associated 
with  airspeed  10  knots  below  Vref;  or 

d)  a   "well   below  glideslope"   call  on   a   PAR   associated 
with  airspeed  10  knots  below  VrCf. 

100 


Escape  procedure: 

1)  Apply  maximum  power  (power  levers  to  the  stops). 

2)  Set  and  maintain  a  pitch  attitude  of  10°  on  the  AI.  Do 
not  attempt  to  recover  airspeed. 

3)  Once  a  positive  rate  of  climb  is  established,  select 
landing  gear  up. 

4)  Do  not  raise  the  flaps  until  the  airspeed  has  increased 
above  140  knots  indicating  exit  from  the  windshear. 

2.  Takeoff  Microburst  Encounter 

Precautions: 

1)  Delay  the  takeoff  if  able. 

2)  Select  the  longest  suitable  runway. 

3)  Perform  takeoff  planning  for  adverse  conditions  as 
prescribed  by  the  NATOPS  Flight  Manual. 

4)  Increase  rotate  airspeed  to  140  knots,  or  the  airspeed 
that  is  attained  with  2000  feet  of  runway  remaining. 
This  is  determined  by  using  the  Four  Engine 
Acceleration  Chart  in  NATOPS  and  following  this 
procedure: 

a)  Subtract  two  thousand  feet  from  the  available 
runway  length. 

b)  Enter  the  chart  with  the  adjusted  runway  length, 
pressure  altitude,  temperature  and  gross  weight.  Exit 
the  chart  with  the  corresponding  airspeed  value.  Use 
140  knots  if  the  chart  value  is  higher. 

c)  Corrections  for  runway  slope,  winds,  or  standing 
water  is  not  required  for  THIS  prediction. 

5)  Thoroughly  brief  the  takeoff  procedure,  voice  calls,  and 
windshear  indications  among  the  flight  crew. 
Penetration  of  a  windshear  during  takeoff  will  be 
indicated  by  a  loss  of  airspeed  or  no  airspeed 
acceleration. 
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Procedure: 

1)  Abort  the  takeoff  if  windshear  is  indicated  before 
reaching  refusal  airspeed.  Reaching  refusal  speed 
ground  roll  distance  with  airspeed  significantly  below 
OR  above  refusal  speed  indicates  windshear. 

2)  Continue  the  takeoff  if  windshear  is  experienced  after 
refusal  speed.    Rotate  the  nose  when: 

a)  predicted  increased  rotate  speed,  usually  140  knots,  is 
attained,  or 

b)  2000  feet  of  runway  remain. 

3)  Rotate  to  a  pitch  attitude  of  10°.  DO  NOT  delay  rotating 
the  nose  because  of  low  airspeed. 

4)  Increase  pitch  attitude  toward  15°  if  ground  impact  is 
imminent. 

5)  Raise  the  landing  gear  when  above  100  feet  AGL  with  a 
positive  rate  of  climb. 

6)  Do  not  raise  the  flaps  until  normal  climb  airspeed  is 
regained,  indicating  clear  of  the  microburst. 

3.  On-Station  Loiter  Microburst  Encounter 

No  immediate  effect  on  the  aircraft  performance  was 
observed.  However,  certain  precautions  should  be  taken  and 
procedures  followed  if  flying  through  convective  activity  at  low 
altitude: 

1)  Use  the  autopilot  altitude  hold  in  dual  axis  mode. 
Note  that  if  the  autopilot  is  not  used,  the  pilot  can 
expect  abnormal  elevator  control  force  changes. 

2)  Minimize  power  lever  movement.  Do  not  pull  off 
power  for  high  airspeed  indication. 
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3)  A  penetration  of  a  microburst  will  be  indicated  by  a 
significant  loss  of  airspeed.  The  first  and  foremost 
reaction  should  be  the  addition  of  full  power.  Monitor 
autopilot  input  and  altitude  hold. 

4)  If  a  climb  is  deemed  warranted,  set  10°  pitch  attitude 
using  the  CWS  function  of  the  autopilot. 


B.  MICROBURST  ESCAPE  PROCEDURES  FOR  THE  T-44 

Target  pitch  attitude  escape  procedures  are  just  as  effective  for  a 
light  turboprop  as  the  T-44.  In  its  case,  the  target  pitch  values  are 
different  between  takeoff  and  landing.  This  change  is  owed  to 
different  flap  settings.  Many  of  the  same  precautions  and 
procedures  are  the  same  between  the  T-44  and  P-3. 

1.  Approach  To  Landing  Microburst  Encounter 
Precautions: 

1)  Avoid  thunderstorm  conditions.  Delay  the  approach  if 
possible. 

2)  Use  Approach  flap  setting.  Do  not  extend  flaps  to  full 
down  position  unless  required  and  until  the  runway  is 
made. 

3)  Use  a  precision  approach  procedure. 

4)  The  approach  speed,  VrCf,  should  be  equal  to  120  knots. 
Additional  airspeed  is  not  warranted.  This  Vret  will 
provide  a  comfortable  margin  above  stall  speed.  Note 
that  additional  airspeed  increases  landing  distance  (5 
knots  faster,  10$  farther). 

5)  Determine  expected  descent  rate. 
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6)  Consider  using  the  autopilot  in  the  coupled  mode  if 
executing  an  ILS.  However,  close  monitoring  of  the 
autopilot  is  required  because  of  abnormal  elevator  pitch 
forces  expected  in  a  windshear.  Insure  autopilot 
disconnect  passing  200  feet  AGL. 

7)  Attain  a  stabilized  airspeed  approach  before  passing 
1000  feet  AGL.  Minimize  power  lever  movement 
beyond  this  point.  Maintain  glideslope  with  pitch.  The 
airspeed  indicator  will  serve  as  a  windshear  indicator. 

8)  Strong  consideration  should  be  given  to  executing  the 
microburst  escape  procedure  if  one  of  the  indications  is 
observed: 

a)  a  rapid  and  sustained  airspeed  loss  of  15  knots  below 

b)  a  descent  rate  500  feet  per  minute  greater  than  the 
predetermined  value; 

c)  greater  than  1  dot  low  on  the  ILS  glideslope  associated 
with  airspeed  10  knots  below  Vref;  or 

d)  a   "well  below  glideslope"  call  on   a  PAR  associated 
with  airspeed  10  knots  below  Vrcf. 

Escape  procedure: 

1)  Apply  maximum  power  (power  levers  to  the  stops). 

2)  Set  and  maintain  a  pitch  attitude  of  10°  on  the  AI.  DO 
NOT  attempt  to  recover  airspeed. 

3)  Once  a  positive  rate  of  climb  is  established,  select 
landing  gear  up. 

4)  Do  not  raise  the  flaps  until  the  airspeed  has  increased 
above  120  knots,  indicating  exit  from  the  windshear. 
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2.  Takeoff  Microburst  Encounter 

Precautions: 

1)  Delay  the  takeoff  if  able. 

2)  Select  the  longest  suitable  runway. 

3)  Perform  takeoff  planning  for  adverse  conditions  as 
prescribed  by  the  NATOPS  Flight  Manual. 

4)  Increase  rotate  airspeed  to  120  knots,  or  the  airspeed 
that  is  attained  with  2000  feet  of  runway  remaining. 

5)  Thoroughly  brief  the  takeoff  procedure,  voice  calls,  and 
windshear  indications  among  the  flight  crew. 
Penetration  of  a  windshear  during  takeoff  will  be 
indicated  by  a  loss  of  airspeed  or  no  airspeed 
acceleration. 

Procedure: 

1)  Abort  the  takeoff  if  windshear  is  indicated  before 
reaching  91  knots  airspeed. 

2)  Continue  the  takeoff  if  windshear  is  experienced  after  91 
knots  airspeed.    Rotate  the  nose  when: 

a)  120  knots;  or 

b)  2000  feet  of  runway  remain. 

3)  Rotate  to  a  pitch  attitude  of  15o.  DO  NOT  delay  rotating 
the  nose  because  of  low  airspeed. 

4)  Increase  pitch  attitude  toward  20°  if  ground  impact  is 
imminent. 

5)  Raise  the  landing  gear  when  above  100  feet  AGL  with  a 
positive  rate  of  climb. 
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V.    RECOMMENDATIONS 

Recommendations  are  broken  into  two  sections.  The  first  set  are 
for  analysis  refinement  and  further  research.  The  second  set  are 
directed  to  the  P-3  operators. 

A.  ANALYSIS  REFINEMENT  AND  CONTINUED  RESEARCH 

There  were  some  weak  areas  of  this  analysis  that  were  pointed 
out  earlier.  Although  valid  results  were  obtained,  refinement  of  the 
windshear  model  and  the  aircraft  equations  of  motion  would  further 
explore  the  edges  of  the  operating  envelopes.  The  following 
recommendations  are  provided  for  the  microburst  windshear  model: 

1)  Vortex  wind  field  models  seem  to  provide  the  closest 
representation  of  a  microburst  windshear,  short  of  applying 
Navier-Stokes  principles.  Schultz's  multiple  vortex  model 
combined  with  source  flows  should  be  combined  and  fitted  to 
other  previously  recorded  microbursts. 

2)  A  more  sophisticated  parameter  estimation  scheme  should  be 
applied  to  enhance  and  accelerate  parameter  fits  of  recorded 
data  to  the  vortex/source  flow  model. 

3)  Realistic  boundary  conditions  to  the  windshear  model  should  be 
defined.  This  would  help  explore  the  entry/exit  aircraft 
response. 

The   following   recommendations   are   for    the   aircraft   performance 
models: 
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1)  The  lift  and  drag  equations  should  be  improved  from  second  to 
fourth  order  equations.  This  change  will  enhance  the 
representation  of  the  aircraft  performance  at  low  airspeed/high 
angle  of  attack. 

2)  An  available  thrust  algorithm  should  be  integrated  into  the 
aircraft  performance  model.  Comparison  of  turbojet,  turbofan, 
turboprop,  and  reciprocating  engine  propulsion  may  lead  to 
significant  effects  on  performance  through  windshear. 

3)  A  state-space  matrix  of  aircraft  dynamics  should  be  integrated 
into  the  performance  algorithm.  Some  aircraft  may  exhibit 
unexpected  dynamics  associated  with  microburst  windshear 
and  escape  maneuvers. 


B.  P-3  OPERATIONS  AND  FLIGHT  CREW  TRAINING 

An  engineering  analysis  provides  the  performance  specifications 
for  penetrating  a  microburst  and  escaping.  However,  certain  issues 
must  be  addressed  by  the  P-3  community  so  as  to  successfully 
employ  the  results: 

1)  Flight  crew  coordination  training  is  essential  for  successful 
employment  of  microburst  escape  procedures.  All  members  of 
the  flight  deck  must  understand  the  teamwork  required  to 
execute  these  procedures.  Therefore,  comprehensive  windshear 
training  should  be  developed.  The  FA  A  Windshear  Training  Aid 
and  this  thesis  provides  an  initial  foundation. 

2)  A  windshear  algorithm  should  be  developed  for  the  P-3  Flight 
Simulators.  There  Is  no  effective  inflight  means  to  expose 
flight  crews  to  microburst  effects.  Nor  can  escape  procedure  be 
effectively  practiced.  The  only  proven  instructional  means  are 
with  a  flight  simulator. 

3)  A  brief  dissertation  on  microburst  and  gust  front  windshear 
should  be  provided  in  the  flight  station  NATOPS  manual. 
NATOPS  change  recommendations  reflecting  the  conclusions  of 
this  study  should  be  adopted. 
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APPENDIX  A 
WINDSHEAR  MODEL  ALGORITHMS 

The  following  are  MatLab®  functions  which  model  the  winds 
experienced  at  DFW  on  August  2,  1985  by  two  separate  aircraft. 
AA539  penetrated  the  microburst  at  2500  feet.  DAL191  penetrated 
the  same  microburst  while  attempting  a  landing. 

A.  AA539  WINDSHEAR  MODEL  ALGORITHM 


funct  ion  [vXjVzl-microburstUejx) 


X    INITIAL  UflLUES  AND  INPUTS 

X  Earth  posit  ion 

Xze  -  altitude  (ft) 

Xx  horizontal  position  (ft) 

X  Uortex  Ring  Dimensions 

AL-8503.3; 

riL-2004,1; 

kL-131968.8; 


ro»L-riL/0.371j 

yL*3350.1; 
LUoralt-3400.6; 
Rs-1701 .7; 
ris«323.9j 
ks-57204.9; 

ro»s»r is/0. 371; 

ys-830.9; 

sUoralt-2333.6; 


Xradius  of  ring  filament  (ft) 

Xcore  radius  (ft) 

X vortex  strength  of  the  ring 

Xfi latent* ft A2/sec) 

Xrelated  core  radius 

X large  ring  y  lateral  position  (ft) 

X large  vortex  altitude  AGL  (ft) 

Xradius  of  ring  fi latent  (ft) 

Xcore  radius  (ft) 

Xvortex  strength  of  the  ring 

Xfi latent  (ftA2/sec) 

Xrelated  core  radius 

Xstall  ring  y  lateral  position  (ft) 

Xstall  vortex  altitude  AGL  (ft) 
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X    PRinflRV  LARGE  RING  CALCULATIONS 
z«ze-LUoralt j 
if  x«=0, 

if  yL«0, 

vxLB0; 

vzL=-kL*RLA2/(2*(RLA2+zA2)A1.5); 

zerock=0; 
end 
else 

zerock=999; 
r«sqrt(xA2+yLA2)j 
r1-sqrt(zA2+(r-RL)A2); 


r2*sqrt(zA2+(r+RL)A2); 


Xradial  distance  from  the  z  axis 
Xclosest  distance  to  the  point  of 
X  interest  from  the  ring  filament. 
{farthest  distance  to  the  point  of 
X  interest  from  the  ring  filament. 
Xscal ing  term 


Iambda«(r2-r1)/(r2+r1); 
n1=(r+RL)/r2; 
n2=(r-RL)/H; 
sig1sz/r2+z/rl; 
sig2sz/r2-z/r1; 

tom»0.75*sqrt(1-lambdaA2)+0.25; 
delta1«(kL*lambda)/(pi*to»); 

de I ta2-(0 . 2955*kL*  i  ambdaA3)/(p  i  *tomA2*sqrt (1-1 ambdaA2) ) ; 
dSidr*0.786*delta1*(n1-n2)-0.394*lambda*delta1*(n1+n2) 
+delta2*((n1-n2)-lambda*(n1+n2)); 

Xderi votive  of  the  stream  function  trt  r 
dSidz=deita1*(0.?68*3ig2-0.394*lambda*sig1) 
+de I ta2*(s  i g2- 1 ambda*s  i  gt ) ; 

^derivative  of  the  stream  function  mrt  z 
dmpL=(1-exp(-r1A2/(0.1*romLA2))); 

X velocity  damping  factor 
vxL*(x/r)*(1/r)*dSidz;     Xhorizontal  velocity  in  x  direction 
vzL«-(1/r)*dSidr;        Ivertical  velocity 
end  Xif 
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X  PRinflRV  SHALL  RING  CRLCULRTIONS 

zBze-sUoralt; 
i  f  xBB0, 
if  ysBB0, 
vxs=0; 

vzsB-ks*RsA2/(2*(RsA2+zA2)A1.5); 
end 
else 

rssqrt((x-50)A2+ysA2);       Iradial  distance  from  the  z  axis 
HBsqrt(zA2+(r-Rs)A2);       Iclosest  distance  to  the  point  of 

I  interest  from  the  ring  filament, 
r2»sqrt(zA2+(r+Rs)A2);       Ifarthest  distance  to  the  point  of 

X interest  from  the  ring  f i lament . 
lambdaB(r2-r1)/(r2+H);       Iscaling  term 
n1B(r+Rs)/r2; 
n2-(r-Rs)/r1; 
sig1Bz/r2+z/r1; 
sig2Bz/r2-z/H; 

t  o««0 . 75*sqrt (1-1 ambdaA2)+0 . 25; 
delta1«(ks*lanbda)/(pi*tom); 

de I t a2-(0 , 2955*ks* I ambdaA3)/(p I *tomA2*sqrt (1-1 ambdaA2) ) ; 
dSidrB0.788*delta1*(n1-n2)-0.394*lambda*delta1 
*(n1+n2)+delta2*((n1-n2)-lambda*(n1+n2)); 

Iderivative  of  the  stream  function  vrt  r 
dSidzBdelta1*(0.788*sig2-0.394*lambda*sig1) 
+de I ta2*(s  i g2- 1 ambda*s  i  gl ) ; 

Iderivative  of  the  stream  function  mrt  z 
if  (rP2/(0.1*roMA2))>441, 

dmpss1; 
else  Ivelocity  damping  factor 

dmps-1-exp(-r1A2/(0. 1*romsA2)); 
end 

vxs*(x/r)*(1/r)*dSidz;       ^horizontal  velocity  x  direction 
vzs"-(1/r)*dSidr;  Xvertical  velocity 

end  Xif 
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X  MAGE  LARGE  RING  CALCULATIONS 

x«x; 

z=-ze-LUoralt; 

if  x«0, 

if  yL==0, 

vxLi=0; 

vzLi«-kL*ALA2/(2*(RLA2+zA2)A1.5); 
end 
else 

r»sqrt(xA2+yLA2); 
M*sqrt(zA2+(r-RL)A2); 


r2«sqrt(zA2+(r+RL)A2); 


Xradial  distance  from  the  z  axis 
Xclosest  distance  to  the  point  of 
X interest  from  the  ring  fi latent. 
Xfarthest  distance  to  the  point  of 
X  interest  from  the  ring  filament. 
Iscal ing  tern 


Iambdaa(r2-r1)/(r2+r1); 

n1*(r+RL)/r2j 

n2«(r-RL)/r1; 

sig1«z/r2+z/r1; 

sig2"z/r2-z/M; 

t o«-0 . 75*sqrt (1-1 ambdaA2)+0 . 25; 

delta1=(kL*lambda)/(pi*tom)j 

de I ta2«(0 . 2955*kL* I ambdaA3)/(p  i  *t omA2*sqrt (1-1 ambdaA2) ) ; 

dSidr«0.?88*delta1*(n1-n2)-0.394*lambda*delta1 

*(n1+n2)+delta2*((n1-n2)-lambda*(n1+n2)); 

Xderi votive  of  the  stream  function  «rt  r 
dS  i  dz*de 1 1  a 1  * ( 0 . 768*s  i  g2-0 . 394* I ambda*s  i  g 1 ) 

+de 1 1  a2* (s  i  g2- 1 ambda*s  i  g 1 ) ; 

Xderi vat ive  of  the  stream  function  mrt  z 
dmpLi-(1-exp(-r1A2/(0.1*romLA2))); 

Xvelocity  damping  factor 
vxLi«(x/r)*(1/r)*dSidz;      Xhorizontal  velocity  in  x  direction 
vzLi«-(1/r)*dSidr;  Xuertical  velocity 

end  Xif 
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X     IflRGE  SMALL  RING  CRLCULRTIONS 
zs-ze-sUoralt; 
if  x«0, 
if  y3=sO, 
vxsisO; 

vzsi*-ks*RsA2/(2*(RsA2+zA2n.5); 
end 
else 

rBsqrt((x+50)A2+ysA2);       Xradial  distance  from  the  z  axis 
r1=sqrt(zA2+(r-Rs)A2);       Xclosest  distance  to  the  point  of 

X  interest  from  the  ring  filament. 
r2=sqrt(zA2+(r+Rs)A2);       Xfarthest  distance  to  the  point  of 

Xinterest  from  the  ring  filament. 
Iambda=(r2-r1)/(r2+r1 );       Xscaling  term 
n1=(r+Rs)/r2; 
n2*(r-Rs)/r1; 
sig1-z/r2+z/r1; 
sig2»z/r2-z/r1; 

t  oipO . 75*sqrt (1-1 ambdaA2)+0 . 25; 
deltaHks*lambda)/(pi*tom); 

de I ta2«(0 . 2955*ks* I ambdaA3)/(p  i  *tomA2*sqrt (1-1 ambdaA2) ) ; 
dSidr-0.788*delta1*(n1-n2)-Q.394*lambda*delta1 
*(n1+n2)+delta2*((n1-n2)-lambda*(n1+n2))j 

Xderivative  of  the  stream  function  mrt  r 
dS  i  dz=de 1 1  a 1  * ( 0 . 788*s  i  g2-0 . 394* I ambda*s  i  g 1 ) 
♦delta2*(sig2-lambda*sig1); 

Xderivative  of  the  stream  function  trt  z 
if  (rr2/(Q.1*romsA2))>411, 

dmpsi-1; 
else  Xvelocity  damping  factor 

dmps I -1 -exp(-r 1A2/(0 . 1 *romsA2) ) ; 
end 

vxsis(x/r)*(1/r)*dSidz;       Xhorizontal  velocity  x  direction 
vzsi«-(1/r)*dSidr;  Xvertical   velocity 

end  Xif 
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X  SUMPTION 

if  zerock==999, 

dmp=dmpL*dmps*dff!pL i  *dnps  i ; 

uxs-(dmp*(uxL+uxs+uxLi+uxsi )); 

vzBdmp*(uzL+uzs-uzL  i -uzs  i ) ; 
else 

ux«-( vxL+uxs+uxL  i  *vxs  i ) ; 

uzsuzL+vzs-uzL  i -uzs  i ; 
end 


B.  DAL191   WINDSHEAR  MODEL  ALGORITHM 

funct  ion  [ux,uz]=microburst(ze,x) 

X    INITIAL  UflLUES  AND  INPUTS 

X  Earth  posit  ion 

Xze  =  altitude  (ft) 

Xx  horizontal  position  (ft) 


X  Uortex  Ring  Dimensions 

RL«7000; 

riL*2004.1; 

kL«43 1968. 8*1; 

rott-riL/0.371; 

yL--300; 

LUoralt-3400; 

Lxdisp=2500; 

Ra-1300; 

ris-323.9; 
ks=57204.9*2.3; 

ro»ssr is/0. 371 ; 
ys»1; 

slloralt-800; 
sxdispB300; 


Xradius  of  ring  f  i  I  anient   (ft) 

Xcore  radius  (ft) 

Xuortex  strength  of  the  ring 

Xfi lament  (ftA2/sec) 

Xrelated  core  radius 

X large  ring  y  lateral  position  (ft) 

X large  uortex  altitude  RGL  (ft) 

X large  ring  displacement  in  x  (ft) 

Xradius  of  ring  filament  (ft) 

Xcore  radius  (ft) 

Xuortex  strength  of  the  ring 

Xfi lament  (ftA2/sec) 

Xrelated  core  radius 

Xsmall  ring  y  lateral  position  (ft) 

Xsmall  uortex  altitude  RGL  (ft) 

Xsmall  ring  displacement  in  x  (ft) 


115 


X    PRinflRV  LARGE  RING  CALCULATIONS 
zBze-LUoralt; 
if  x«0, 
if  yL«0, 

vxL=0; 

uzL=-kL*RL"2/(2*(RL"2+z"2)"1.5); 

zerock*0; 
end 
else 

zerock=999; 
r*sqrt(x*2+yL*2); 
M«sqrt(zA2+(r-RL)A2); 


r2-sqrt(z"2+(r+RL)A2); 


Xradial  distance  from  the  z  axis 
Xclosest  distance  to  the  point  of 
I  interest  froi  the  ring  fi latent. 
Ifarthest  distance  to  the  point  of 
X interest  from  the  ring  filament. 
Xscal ing  tern 


Ianbda«(r2-r1)/(r2+r1); 
n1*(r*RL)/r2; 
n2-(r-AL)/r1; 
sig1*z/r2+z/r1; 
sig2sz/r2-z/r1; 

tof-0.75*sqrt(1-laibda"2)+0.25; 
delta1=(kL*laibda)/(pi*to»); 

de I ta2«(0 . 2955*kL* I aibdaA3)/(p  i  *toiA2*sqrt (1-1 aibda*2) ) ; 
dSidr«0.788*delta1*(n1-n2)-0.394*laibda*delta1*(n1+n2) 
+delta2*((n1-n2)-laibda*(n1*n2))j 

Xderiuative  of  the  strean  function  irt  r 
dSidz«delta1*(0.788*sig2-0.394*laibda*sig1) 
+delta2*(sig2-la«bda*sig1); 

Xderi votive  of  the  streai  function  »rt  z 
dtpL-Cl-expC-rl^/CO.I'roiL^))); 

X velocity  daiping  factor 
vxL*(x/r)*(1/r)*dSidz;     Xhorizontal  velocity  in  x  direction 
vzL*-0/r)*dSidr;        Xvertical  uelocity 
end  Xif 
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X    PRinflRV  SHALL  RING  CALCULATIONS 
zBze-sUoralt; 
if  x"0, 
If  ys«0, 
vxs=0; 

vz8*-ks*RsA2/(2*(RsA2+zA2)A1.5); 
end 
else 

r=sqrt((x-50)A2+ysA2);       Iradial  distance  from  the  z  axis 
r1*sqrt(zA2+(r-Rs)A2);       Xclosest  distance  to  the  point  of 

^interest  from  the  ring  filament. 
r2*sqrt(zA2+(r+Rs)A2)j       Xfarthest  distance  to  the  point  of 

linterest  from  the  ring  filament, 
Iambdas(r2-r1)/(r2+r1);       Iscaling  ten 
n1*(r+Rs)/r2; 
n2-(r-Rs)/r1; 
sig1=z/r2+z/r1; 
sig2*z/r2-z/r1; 

t  om»0 . 75*sqrt (1-1 ambdaA2)+0 . 25; 
delta1B(ks*lambda)/(pi*tom); 

de I ta2»(0 . 2955*ks* I ambdaA3)/(p  i  *tomA2*sqrt (1-1 ambdaA2) ) ; 
dSidr-0.788*delta1*(n1-n2)-0.394*lambda*delta1 
*(n1+n2)*delta2*((n1-n2)-lambda*(n1*n2)); 

Ideriuative  of  the  stream  function  «rt  r 
dSidz«delta1*(0.788*sig2-0.394*lambda*sig1) 
*delta2*(sig2-lambda*sig1); 

Iderivative  of  the  stream  function  art  z 
if  (rr2/(0.1*roM*2))>111, 

dmpsB1; 
else  ^velocity  damping  factor 

dmps-1 -exp(-r 1 A2/(0 , 1 *rotsA2) ) ; 
end 

vxs«(x/r)*( 1/r )*dSidz;       Ihorizontal  velocity  x  direction 
vzs«-(1/r)*dSidr;  ^vertical  velocity 

end  Xif 
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X     IttflGE  LRRGE  RING  CRLCULTIONS 

X8Xj 

z=-ze-LUoralt; 
if  x«0, 
if  yL«0, 
vxLi*0; 

vzLi=-kl_*RLA2/(2*(RLA2+zA2)A1.5); 
end 
else 

rssqrt(xA2+yLA2);  Xradial  distance  from  the  z  axis 

r1Bsqrt(zA2+(r-RL)A2);       Xclosest  distance  to  the  point  of 

Xinterest  from  the  ring  filament. 
r2»sqrt(zA2+(r+RL)A2);       Xfarthest  distance  to  the  point  of 

X interest  from  the  ring  filament. 
Iambdas(r2-r1)/(r2+r1);       Xscallng  term 
n1-(r+RL)/r2; 
n2«(r-RL)/M; 
sig1*z/r2+z/r1; 
sig2=z/r2-z/M j 

tom-0 . 75*sqrt (1-1 ambdaA2)+0 . 25; 
delta1*(kL*lambda)/(pi*toi); 

de I ta2«(0 . 2955*kL* I ambdaA3)/(p  i  *tomA2*sqrt (1-1 ambdaA2) ) ; 
dSidr«0.788*delta1*(n1-n2)-0.394*lambda*delta1 
*(n1+n2)+delta2*((n1-n2)-lambda*(n1+n2)); 

Xderivative  of  the  stream  function  trt  r 
dS  i  dz-de 1 1  a 1  * ( 0 . 788*s  i  g2-0 .394*1 ambda*s  i  g 1 ) 
+delta2*(sig2-lambda*sig1); 

Xderivatiue  of  the  stream  function  mrt  z 
dmpLi-(1-exp(-r1A2/(0.1*romLA2))); 

Xvelocity  damping  factor 
uxLi«(x/r)*(1/r)*dSidzj      Xhorizontal  velocity  in  x  direction 
uzl_i«-(1/r)*dSidr;  Xuertlcal  velocity 

end  Xif 
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%  IHRGE  SMALL  RING  CRLCULRTIONS 

z=-ze-sUoralt; 
if  x--Q, 
i  f  ysBB0, 
vxsi=0; 

vz3i=-ks*RsA2/(2*(RsA2+zA2)A1.5); 
end 
else 

r»sqrt((x+50)A2+ysA2);       Iradial  distance  from  the  z  axis 
r1asqrt(zA2*(r-Rs)A2);       Iclosest  distance  to  the  point  of 

I  interest  from  the  ring  filament. 
r2ssqrt(zA2+(r+Rs)A2);       Ifarthest  distance  to  the  point  of 

I  interest  from  the  ring  filament. 
Iambda»(r2-r1)/(r2+r1);       Iscaling  term 
nl«(r+Rs)/r2; 
n2«(r-Rs)/H; 
sig1»z/r2+z/r1; 
sig2sz/r2-z/r1 ; 

tom-0 . ?5*sqrt (1-1 ambdaA2)+0 . 25; 
deltaWks*lambda)/(pi*toi); 

delta2-(0.2955*ks*lambdaA3)/(pi*toi"2*sqrt(1-lambdaA2)); 
dSidr*0.?88*delta1*(n1-n2)-0.394*lambda*delta1 
*(n1+n2)+delta2*((n1-n2)-lambda*(n1+n2)); 

Iderivative  of  the  stream  function  vrt  r 
dSidz«delta1*(0.788*sig2-0.394*lambda*sig1) 
+de I ta2*(s  i  g2- 1 ambda*s  i  gl ) ; 

Xderi votive  of  the  stream  function  mrt  z 
if  (r 1*2/(0. 1*romsA2))>444, 

dmpsis1; 
else  ^velocity  damping  factor 

dmps  i -1 -exp(-r P2/(0 . 1 *romsA2) ) ; 
end 

vxsi«(x/r)*(1/r)*dSidz;      ^horizontal  velocity  x  direction 
vzsi«-(1/r)*dSidr;  ^vertical   velocity 

end  lif 
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X    SOURCE  PROFILE  FOR  UZ  CORRECTION 

Salt =10000;  Xsource  altitude 

gagazs2772000;  Xsource  strength/2i 

gagaxe 1355396;  Xsource  strength/2*  far 

uzsourceB0.0; 

uxsources0.0; 

if  ze>800, 

vzsource*-15*(LUoralt-zc)/((LUoralt-ze)A2+xA2)A0.5; 

else 

Radius1»(xA2+(Salt-ze)A2)A0.5; 

Xradius  fron  priuary  source  to  the  flight  pt 
Radius2»(xA2+(Salt+ze)A2)A0.5; 

Xradius  from  inage  source  to  the  flight  pt 
uzsourceagagaz*( (Sa I t+ze)/Rad  i  us2A2-(Sa I t-ze)/Rad  i  us1 A2) ; 
end 
vxsource=gagax*x*( 1 /Rad  i  usl A2+1 /Rad  i  us2A2) ; 

X  SUrWflTION 

if  zerock—999, 

dip»dnpL*dtps*dipL i  *dips  i ; 

if  x<«2000, 
vx« ( dip* ( vxL+vxs+vxL  i  +vxs  i ) ) ; 

else 
ux«vxsource; 

end 

uz«-(dip*( vzL+wzs-vzL i -vzs i )-uzsource) ; 
else 

ux«-( vxL+uxs+vxL I +VX3 I ) ; 

vz«vzL+vzs-vzL I -uzs I ; 
end 
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APPENDIX  B 
PARAMETER  SENSITIVITY  ALGORITHM 

The  following  are  MatLab®  m-file  and  function  algorithms  used  to 
determine  the  parameters  for  the  windshear  model  that 
approximates  the  DAL191  recorded  winds. 

A.  M-FILE  DRIVER 

IParameter  sensitivity  program  for  the  DRL1 91  wind  model. 

clear 

clc 

IRecorded  flight  data  from  DRL191 . 

X[ref  nuaber,  alt  itude(ft  AGO,  diet  fra  MB  center(ft),  x  wind 

(ft/s),  y  find  (ft/a)]. 


DflL«[1 

987  -8642 

2  1 

0; 

2  975 

-8361 

-1 

Hi 

3  962 

-8080 

-3 

11, 

4  948 

-7799 

-4 

11, 

5  935 

-7518 

-6 

11 

6  922 

-7237 

-9 

15, 

7  909 

-6956 

-12 

12, 

8  894 

-6676 

-12 

3, 

9  875 

-6396 

-9 

6, 

10  856 

-6117 

-11 

9, 

11  842 

-5838 

-16 

1i 

12  823 

-5561 

-25 

-18, 

13  807 

-5287 

-36 

-11, 

14  790 

-5016 

-39 

-13, 

15  777 

-4749 

-29 

-20, 

16  764 

-4486 

-26 

-17, 

17  755 

-4229 

-22 

-17, 
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18  747 

-3977 

-29 

-17 

19  742 

-3730 

-20 

-11 

20  738 

-3487 

-17 

-23 

21  730 

-3249 

-12 

-30 

22  714 

-3013 

-8 

-36 

23  698 

-2778 

0 

-34 

24  680 

-2542 

3 

-32 

25  662 

-2306 

0 

-31 

26  647 

-2069 

-8 

-31 

27  635 

-1831 

-18 

-24 

28  627 

-1593 

-20 

-26 

29  620 

-1355 

-10 

-36 

30  608 

-1117 

-8 

-40 

31  592 

-880 

18 

-4 

32  585 

-639 

32 

19 

33  582 

-391 

19 

16 

34  580 

-139 

18 

21, 

35  573 
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19 

-3, 

36  548 

385 

38 

-21, 

37  511 

656 

54 

-19, 

38  462 

936 

65 

9, 

39  411 

1228 

66 

9, 

40  351 

1531 

72 

-13, 

41  277 

1845 

71 

-42, 

42  193 

2169 

72 

-33, 

43  119 

2504 

76 

-2, 

44  70 

2851 

78 

15, 

45  42 

3211 

64 

17. 

46  20 

3579 

81 

4, 

47  0 

3954 

85 

2]; 

XPar dieter  aatri 

x: 

Ps«[2300 

2500 

1250 

1300 

57204.9*2  57204.9*2.3 

-400 

-300 

300 

500 

-150 

1 

2800;  XLarge  ring  x  displacement 

1 350 ;  XSaa II  r  i  ng  core  r ad  i  us 

57204.9*2.5;  XSaa 1 1  ring  vortex  strength 

-200;  XLarge  ring  y  posit 

700;  XSaa 1 1  ring  x  displacement 

150;  XSaall  ring  y  displacement 
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1386000*1.8     1386000*2     1386000*2.2;  XSource  z  strength 

193628*6.5        193628*7       193628*7.5];         XSource  x  strength 


yxerrs0.0;vzerr,s0.0;»nB0;bestrisa999;       linit  ial  izat  ion 
iB1 ; 

3=1; 

for  k«1 :3, 
hoie 
i 

z 
k 

bestris 

for  l«1;3,   for  ■■IjS,   for  n«1:3, 
for  o«1 :3,   for  p«1 :3, 
vxer r -0 . 0 ; vzerr "0 . 0 ; 
for  c«1 :47 , 
[Ux,Uh]-piicroburst(DRL(c,2),DRL(c,3),Ps(1, i),Ps(2,z), 

Ps(3,k),P3(4J),P3(5,n),Ps(6,n),P3(7,o),Ps(8,p)); 
vxerr-uxerr+CDflUcjD-Ux)^; 
vzerr-vzerr+(DflL(c,  5)-Uh)*2; 
end 

risvx-(vxerr/40) A0 .5; 
risvz«(vzerr/40)*0.5; 
nstot-( (vxerr+vzerr )/80)*0 .5; 
■n-in+1; 

Pn(in,:)-[Ps(1J)lPs(2iz)JPs(3,k)iP3(^J), 

Ps(5li),Ps(6>n),Ps(7,o)lPs(8lp)lrisvx,risvz,ristot]; 
if  bestris>ristotj 

BESTP-Pridn,:); 

bestris-nstot; 

save  BEST 
end 

end!  end,  end,   end,   end,   end, 
save  paraietersll 
clear  ix  ih  nsvx  nstot  in  Pfl  BESTP 

vxerr«0 . 0; uzerr«0 . 0; in"0;       I i  n  i  t  i  a  I i  zat  i  on 
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1-1; 
z*2; 

for  k«1:3, 
hone 
i 

z 
k 

bestrns 

for  1—1 : 3 ,  for  ms1:3,  for  n=1:3, 
for  o=1 :3,  for  ps1 :3, 
uxerrxO . 0; uzerr=0 .0; 
for  c=1  'Al , 
[Ux,Uh]«pi>icroburst(DAL(c, 2), DRL(eJ3)lPs(1i l),Ps(2,z), 

Ps(3,k),Ps(4,l),Ps(5,8),Ps(6,n),Ps(7,o),Ps(8,p)); 
vxerr-uxerr+(DflL(c,4)-Ux)A2; 
vzerr*uzerr+(DflL(c,5)-Uh)A2; 
end 

risvx«(vxerr/4D) AQ .5; 
risvz«(vzerr/40)A0.5; 
ristots((vxerr+yzerr)/80)A0.5; 
■n-m+1; 
Pn(8n,:)«[Ps(1,i),Ps(2,z),Ps(3,k),P8(4,l), 

Ps(5,i),Ps(6,n),Ps(7,o),P8(8,p),risvx,rtsuz,r«stot]; 
if  beatri8>rietot , 

BESTP-PfKin,:); 

bestns-rietot ; 

save  BEST 
end 

end ,  end,  end,  end,  end,  end, 
save  parameters 12 
clear  ix  »h  rtsvx  rnstot  in  Pfl  BE5TP 

uxerr"0.0;vzerr*0.0;in«0;   Xinitiol ization 

1-1; 

z-3; 

for  k-1:3, 
hone 
i 
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z 

k 

best p»s 

for  1*1:3,  for  8=1:3,  for  n=1:3, 

for  os1  :3,  for  p«1 :3, 

uxerra0 . 0; vzerrB0 .0; 

for  cs1  :47, 

[Ux,Uh]Bpincrobur3t(DRL(c, 2), DRL(c,3),P8(1, i),Ps(2,z), 

P3(3,k),P3(4,l),Ps(5,«),Ps(6,n),P3(7,o),P3(8,p)); 
vxerr=vxerr+(DRL(c,4)-Uxr2; 
vzerrBvzerr+(DRL(c,5)-Uh)A2; 
end 

r»8vx"( vxerr/40) A0 ,5; 
ri3vz«(uzerr/40) "0 .5; 
rwtot  ■(  (uxerr+uzerr  )/80)A0 .5; 
■n«»n+1; 
Pn(8n,:)*[P3(1,i),P3(2,z),Ps(3,k),Ps(4,l), 

P3(5,i),Ps(6,n),Ps(7,o),P3(8,p),r83vx,r»svz,rB3tot]; 
if  bestris>ristot, 

BESTP-PnCin,:); 

be3tri3"nstot; 

save  BEST 
end 

end;  end,  end,  end,  end,  end, 
saue  parameters  13 
clear  8X  in  nsux  ristot  in  Pfl  BESTP 

vxerr-0 . 0; vzerr«0 . 0; in-0;   I i  n  i  t  i  a  I i  zat  i  on 
i-2; 

1-1  i 

for  kB1 :3, 
hoie 
i 

z 
k 

best r is 

for  1*1:3,  for  ■■1:3,  for  n«1 :3, 
for  o«1 :3,  for  p«1 :3, 
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vxerrB0 . 0; gzerrB0 .0; 

for  c=1 :47, 

[Ux,Uh]«p»icrobur3t(DHL(c, 2) ,DflL(c, 3),Ps(1 , i),Ps(2,z), 

P3(3,k),P3(4,l),P3(5,»),Ps(6,n),Ps(7,o),Ps(8,p)); 
vxerr*vxerr+(Dfl!_(c,4)-Ux)A2; 
uzerr=vzerr+(DRL(c,5)-Uh)A2; 
end 

rnsvx«(vxerr/40) A0 .5; 
r«3vz=(uzerr/40)A0.5; 
r»stots((vxerr+vzerr)/80)A0.5; 
«nBin+1; 
Pn(»n,:HP3(1,i),P3(2,z),Ps(3,k),p3(4,l), 

Ps(5,B),Ps(6,n),P3(7,o),Ps(8,p),ri3ux,rm3uz,ristot]; 
if  bestri3>ri3tot, 

BESTP-Pndn,:); 

bestris*ristot; 

save  BEST 
end 

end,  end,  end,  end,  end,  end, 
save  parameter 821 
clear  bx  bH  ridux  mstot  in  Pf1  BESTP 

vxerr«0.0;vzerr"0.0;Bn"0;   Unit ial izat ion 

i-2; 

z-2; 

for  k«1 :3, 
hose 
i 

z 
k 

best r is 

for  1-1:3,  for  ■■1:3,  for  n«1:3, 
for  o«1 :3,  for  p»1 :3, 
vxerr«0 .0; vzerr«0 .0; 
for  c«1 :47, 
[Ux,Uh]-piicrobur3t(DflL(c,2),DRL(c,3),Ps(1,i),P3(2,z), 

P3(3,k),P3(4,l),P8(5,i),P3(6,n),P3(7,o),P3(8,p)); 
vxerr-vxerr+(DflL(c,4)-Ux)A2; 
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uzepp*vzepp+(DRL(c,5)-Uh)A2; 

end 

p«3vxB(vxepp/40)A0 .5; 

p«3vz«(uzepp/40)A0.5; 

p»stot*((vxepp+vzepp)/80)A0.5; 

«n=in+1 ; 

Pn(in,:HP3(1,i),P3(2,z),Ps(3,k),P3(4,l), 

P3(5,n»),P3(6,n),P3(7,o),Ps(8,p),Pisux,PBi3UZ,p«3tot]; 
if  be3tPK3>n»3totl 

BESTP-Pn(«n,:); 

bestrns^ptstot; 

save  BEST 
end 

end,  end,  end,  end,  end,  end, 
save  papaieteps22 
cleap  vx   ih  pisvx  puetot  in  PF1  BESTP 

vxepp«0.0;vzepp«0.0;in«0;   Xinit ial izat  ion 

i«2; 

z-3; 

fop  k«1:3, 
hone 
i 

z 
k 

bestpis 

fop  1*1:3,  fop  ■■1:3,  fop  n«1:3, 
fop  o«1 :3,  fop  p«1 :3, 
wxepp-0 . 0 ; vzepp«0 . 0 ; 
fop  c-1 :47, 
[Ux,Uh]=piicpobupst(DflL(c,2),DflL(c,3),P3(1,i),Ps(2,z), 

Pe(3,k),Pe(4,l),P3(5,i),Pe(6,n),P3(7lo),P3(8,p)); 
vxepp-vxepp+(DRL(c,4)-Ux)A2; 
vzepp«vzepp+(DRL(c,5)-Uh)A2; 
end 

pisvxs(vxepp/40)A0.5; 
PisuzB ( vzepp/40 ) AQ . 5 ; 
Pi3tot«( ( vxepp+vzepp)/80)A0 .5; 
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»ns8n+1 ; 
Pn(in,:MPs(1#l),P8(2,z),P8(3,k),Ps(4,l), 

P3(5,»),Ps(6,n),Ps(7,o),P8(8,p),r«3vx,rin3vz,r»3tot]; 
i  f  bestr»s>ra3tot, 

BESTP«Pn(«n,:); 

bestrBssrB8tot; 

saue  BEST 
end 

end,  end,  end,  end,  end,  end, 
saue  panaroeters23 
clear  «x  »h  rasvx  r»stot  in  PR  BESTP 

vxerrsO.O;vzerrsO.O;omsO;   linit ial izat  ion 
i»3; 

z=1j 

for  ke1 :3, 
hone 
i 

z 
k 

best r is 

for  Ms3,  for  ■■1:3,  for  n«1:3, 
for  o«1 :3,  for  p«1 :3, 
vxerr«0 .0; vzerr«0 .0; 
for  c*1 :47, 
[Ux,Uh]«piicrobur8t(DflL(c,2),DflL(c,3),Ps(1,i),P8(2,z)i 

p8(3,k),P8(4,l),P8(5,i)lP8(6,n),P8(7,o),P8(8,p)); 
vxerr«vxerr+(DRL(c, 4)-Ux) A2; 
vzerr-vzerr+(DflL(c,5)-Uh)A2; 
end 

r 18 vx« ( vxer r /10 ) A0 . 5 ; 
ri3uz«( vzerr/40) A0 .5; 
ri8tot-( ( vxerr+wzerr )/80) A0 .5; 
■n*in+1; 
Pn(in,:)«[p8(1,i),Ps(2,z),P8(3,k),P8(4,l), 

P8(5,i),P8(6,n),P8(7,o),P8(8,p),rt8VX,r«8vz,ristot]; 
if  be8tr88>ri8tot, 
BESTP-Pttdn,!); 
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bestrwristot; 

save  BEST 
end 

end,  end,  end,  end,  end,  end, 
save  parameters31 
clear  wx  wh  r»svx  rastot  »n  Pn  BESTP 

uxerna0.0;vzerps0.0;nnat0;   Unit  ial  izat  ion 

1-3; 

z«2; 

for  ks1 :3, 
hote 
i 

z 
k 

bestms 

for  l«1:3,  for  ■■1:3,  for  n«1:3, 
for  o*1 :3,  for  pB1 :3, 
vxerr«0 .0; vzerrB0 .0; 
for  c«1  :47, 
[Ux,Uh]-piicroburst(DflL(c,2),DflL(c,3),Ps(1,i),Ps(2,z), 

Ps(3,k),Ps(^,l),Ps(5,i),Ps(6,n),Ps(7,o),Ps(8,p)); 
vxerr»vxerr+(Dfll_(c,4)-Ux)A2; 
vzerr«vzerr+(DflL(c,5)-Uh)A2; 
end 

risvx"(vxerr/40)A0.5; 
rtsvz- ( vzerr/40 ) A0, 5 ; 
ristot«((vxerr+vzerr)/80)A0.5; 
■n*tn+1; 
Pf1(in,!)«[Ps(1,i),Ps(2,z),Ps(3,k),Ps(4,l), 

Ps(5,i),Ps(6,n),Ps(7,o),Ps(8,p),risvx,rBsyz,r«stot]; 
if  best ns>ns tot, 

BESTP-Pndn,:); 

bestris«ristot; 

save  BEST 
end 

end,  end,  end,  end,  end,  end, 
save  paraieters32 
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clear  wx  uh  rusux  rastot  an  PR  BESTP 

vxerr=O.Ojuzerr=0.0;im*0;   Xinit  ial izat  ion 
i-3; 

z*3; 

for  ks1 :3, 
hoie 
i 

z 
k 

oestrus 

for  ls1:3,  for  11*1:3,  for  n=1:3, 
for  os1 :3,  for  ps1 :3, 
vxerr«0 . 0 ; vzerr-0 . 0 ; 
for  c«1 :17 , 
[Ux,Uh]-p«icroburst(DflL(c>2),DflL(c>3),Ps(1>  1)^8(2^), 

Ps(3,k),Ps(4,l),Ps(5,ii),Ps(6,n),Ps(7,o),Ps(8,p)); 
vxerr*uxerr+(DRL(c,4)-Ux)A2; 
vzerr- vzer  r  ■►  ( DflL  ( c ,  5 ) -Uh  )*2 ; 
end 

risvx=(vxerr/40)"0.5; 
rts vzB ( uzerr/40 ) A0 . 5 ; 
ristot«((vxerr+vzerr)/80)A0.5; 
■n*in+1; 
Pri(inl:)-[P8(1,l),Ps(2,2)lP8(3,k)lPs(4ll), 

Ps(5,i),P8(6,n),Ps(7,o),P3(8,p),rtevx,rt3vz,r«3tot]; 
if  bestris>ristot, 

BESTP-Pndn,:); 

bestns^rtstot; 

eaue  BEST 
end 

end,  end,  end,  end,  end,  end, 
saue  paraieters33 
clear 

load  BEST 

vxerr«0 .0; vzerr«0 .0; 

for  i«1:47, 
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hone 

i 

X=DflL(i,3); 

[Ux,Uh]=p»icroburst(DflL(i,2),DflL(i,3),BESTP(1)>BESTP(2), 

BESTP(3),BESTP(4),BESTP(5),BESTP(6),BESTP(7),BESTP(8))j 

model  (i, :Hi, Ux,Uh]; 

vxerr=vxerr+(DFIL( i  ^)-Ux)A2; 

uzerr*vzerr+(DRL( i ,5)-Uh)A2; 
end 

p«svx=(vxerr/40)/'0.5; 
rmsvz=(vzerr/40)A0.5; 
nB3tots((uxerr+uzerr)/80)A0.5; 
ansmer«0; 
•hi  le  ansmer~*999, 

dispCThe  horizontal  RflS  is:     ') 

disp(rmsvx) 

dispCThe  vertical  RHS  Is:') 

disp(rmsvz) 

dispCThe  total  system  RMS  is:') 

disp(rmstot) 

dispCThe  following  plots  are  available:') 

dispC  1 "horizontal  minds') 

dispC  2«uertical  minds') 

dispC  999-end') 
ansmer« input ('plot  number  '); 
if  ansmer«*1, 

plot(DRL(:,3),DHL(:,4)) 

hold 

plot(DflL(:,3),iodel(:,2),,~') 

pause 

hold 
else  if  ansmerM2, 

plot(DRL(:,3),DRL(:,5)) 

hold 

plot(DRL(:,3),model(:,3),'--') 

pause 

hold 
end, end 
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B.  WINDSHEAR  FUNCTION 

function  [uxJuz]"piicroburst(zeixiLxdi3piR8^ksJyLJ3xdi3p>RL;kL) 

X     INITIAL  UflLUES  AND  INPUTS 

X  Earth  posit  ion 

Xze  *  altitude  (ft) 

Xx  horizontal  position  (ft) 


X  Uortex  Ring  Di tens  ions 

XRL 

riL«2004.1; 

XkL 

roiL-riL/0.371; 

*yL 

LUoralt-3400; 
XLxdisp 

XRs 

ris-323.9; 

Xks 

roms»r is/0. 371 j 
ys«830. 9/830.9; 
sUoralt-800; 
Xsxdisp 


Xradius  of  ring  filament  (ft) 

Xcore  radius  (ft) 

Xvortex  strength  of  the  ring 

Xfi latent  (ftA2/sec) 

Xrelated  core  radius 

X large  ring  y  laterial  position  (ft) 

X large  uortex  altitude  R6L  (ft) 

X large  ring  displacement  in  x  (ft) 

Xradius  of  ring  filament  (ft) 

Xcore  radius  (ft) 

Xvortex  strength  of  the  ring 

Xfi lament  (ftA2/sec) 

Xrelated  core  radius 

Xsmall  ring  y  laterial  position  (ft) 

Xsmall  uortex  altitude  RGL  (ft) 

Xsmall  ring  displacement  in  x  (ft) 


X    PRIttflRV  LRR6E  RING  CRLCULRTI0N5 

z*ze-LUoralt; 

zerock-999; 

r-sqrt ( (x-Lxd  i  sp) A2+yLA2) ; 

H-sqrt(zA2*(r-RL)A2); 

r2«sqrt(zA2*(r+RL)A2); 

I ambda-(r2-r 1 )/(r2*r 1 ) ; 

n1-(r+RL)/r2; 

n2«(r-RL)/M; 
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sig1Bz/n2+z/H; 

sig2=z/r2-z/r1 ; 

to««0.75*sqrt(1-la»bda"2)+0.25; 

delta1=(kL*la»bda)/(pi*toi); 

de 1 1 a2=(0 . 2955*kL* I a»bdaA3)/(p  i  »to«T2*3qpt (1-1 ambdaA2) ) ; 

dSidp=0.788*delta1*(n1-n2)-0.39^*la»bda*delta1 

»(n1+n2)+delta2*((n1-n2)-la»bda*(n1+n2)); 
dSidz=delta1*(0.788*3ig2-0.394*lainbda*3ig1) 

+de I  ta2*(s  i  g2- 1 ambda*s  i  gl ) ; 
d»pL*(1-exp(-n1A2/(0.1*ro«ir2))); 
vxL*(x/p)*(1/p)*dSidz; 
vzL«-(1/r)*dSidr; 

X    PRIflARY  SHALL  RING  CALCULATIONS 
zeze-sUoralt; 

p=sqpt ( (x-sxd  i  sp) A2+ysA2) ; 
p1=sqpt(zA2+(p-As)A2); 
p2«sqpt(zA2+(p+Rs)A2); 
Iatbda*(r2-p1)/(r2+p1)j 
nl»(r+Rs)/r2; 
n2-(p-Re)/p1; 
3ig1«z/p2+z/p1; 
eig2Bz/p2-z/p1; 

tot-0 . 75*eqpt (1-1 atbdcT2)+Q . 25; 
delta1-(ks«laabda)/(pi*toi); 

de I ta2=(0 . 2955*ks* I a«bda"3)/(p  i  *toi"2*sqpt (1-1 aibdaA2) ) ; 
dSidp*0.788*delta1*(n1-n2)-0.394*la^bda*delta1 
*(n1+n2)+delta2*((n1-n2)-  Iaibda*(n1+n2)); 
dSidz»delta1*(0.788*sig2-0.391*laibda*sig1) 

+delta2*(sig2-laibda*sig1); 
if  (pr2/(0.1*poisA2))>W4, 

dip3"1j 
else 

dipe-1 -exp(-p1A2/(0 . 1 5*POieA2) ) ; 
end 

yxs-(x/p)*(1/p)*dSldz; 
vz3s-(1/p)*dSidp; 
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X  ItlAGE  LARGE  RING  CRLCULRTIOHS 

ze-ze-LUorolt ; 

r=sqrt ( (x-Lxd  i  sp) A2+yL"2) j 

H*3qrt(z'42+(p-RLr2)j 

n2=3qrt(z"2+(r+RL)"2); 

Ia»bdcp(p2-r1)/(n2+r1); 

n1«(r+RL)/r2; 

n2«(r-RL)/r1; 

3igt«z/r2+z/r1; 

sig2sz/r2-z/r1 ; 

to»»0.75*sqrt(1-la»bda"2)+0.25; 

delta1»(kL*latbda)/(pi*toi)j 

delta2«(0.2955*kL*la»bda"3)/(pi*tor2*3qrt(1-laibda"2)); 

dSidp=0.788*delta1*(n1-n2)-0.394*laibda*delta1*(n1+n2)+delta2*((n1-n 

2)-laibda*(n1+n2)); 

dSidz«delta1*(0.788*sig2-0.394*l«bda*sig1) 

+de I ta2*(s  i  g2- 1 aibda*3  i  g1 ) ; 
d»pLi-(1-exp(-pr2/(0.1*roiL'42))); 
vxLi«(x/r)*(1/p)*dSidzj 
uzLi«-(1/r)*dSidr; 

X     IHRGE  SHALL  RING  CRLCULRTIONS 
z«-ze-sUoraltj 
r«sqrt ( (x-sxd I 8p)*2+y3*2) ; 
r1-3qrt(zA2+(r-R3)"2); 
r2-3qnt(z"2+(r+R3)"2); 
Iatbda«(r2-r1)/(r2+r1); 
n1«(r+R3)/r2; 
n2-(r-R3)/r1; 
sig1sz/r2+z/M; 
3lg2«z/r2-z/r1; 

toi-0.75*3qrt(1-laibda"2)+0.25; 
delta1«(ks»latbda)/(pl*toi); 

de  I  ta2=(0 .  2955*ks*  I  aibdaA3)/(p  I  *tor2*3qrt  ( 1  - 1  anbdaA2) ) ; 
dSidr«0.788*delta1*(n1-n2)-0.391*laibda*delta1 
*(n1+n2)+delta2*((n1-n2)-  Iaibda*(n1+n2)); 
dSidz«delta1*(0.788*3ig2-0.394*laibda*3lg1) 
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+de I ta2*(s  i g2- 1 a«bda*s  i  gl ) ; 
if  (rr2/(0.1*roisA2))>444, 
dups  i  =  1 ; 

else 

d»p3i  =  1-exp(-rr2/(0.15*ro»s"2)); 
end 

uxsi=(x/r)*(1/r)*dSidz; 
vzsi»-(1/r)*dSidr; 

X    SOURCE  PROFILE  FOR  UZ  CORRECTION 

Salt-10000;  Xsource  altitude 

Xgagaz  Xsource  strength/2i 

Xgagax  Xsource  strength/2f  far 

vzsource=0.0; 

vxsource=0.0; 

if  ze>600, 

vzsource--15*(LUoralt-ze)/((LUoralt-ze)/7*x*2r0.5; 

else 

Radius1-(x"2+(Salt-ze)"2)"0.5; 

Radius2«(xA2+(Salt+ze)"2)A0.5; 

uzsource-gagaz*( (Sa 1 1 +ze)/Rad  i  us2A2-(Sa I t-ze)/Rad  i  usT2) ; 

end 

vxsourcesgagax*x*(  1  /Rad  i  usP2+1  /Rad  i  us2A2) ; 

x   sunnflTiOH 

dipBdapL*dips*dipL  i  *dips i ; 
If  x<«2000, 

ux-(dip*(uxL+vxs+vxL I +vxs  i ) ) ; 
else 

ux«uxsource; 
end 
uz«-(dip*( vzL+vzs-vzL  i -vzs  i )-uzsource) ; 
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APPENDIX  C 
LANDING  APPROACH  ALGORITHMS 

Contained  in  this  appendix  are  the  MatLab®  programs  utilized  in 
calculating  the  aircraft  response  given  a  windshear  and  a  particular 
escape  maneuver.  They  are  listed  in  three  sections.  The  first  section 
contains  the  driver  routine  and  the  functions  required  to  solve  the 
differential  equations.  The  second  and  third  sections  contain  the 
specific  aircraft  parameter  functions  and  escape  maneuver  functions 
respectively.    The  windshear  function  is  listed  in  Appendix  A. 

A.  DRIVER  AND  DIFFERENTIAL  EQUATION  SOLVER 

X        Descent  Profile  tlodel   * 

X        This  program  ill!  calculate  the  Inert  I  a  I  flight  path  angle  for 
X   a  given  aircraft.  It  utilizes  the  function  '■icroburst'  ihich 
X  lodels  the  find  shear  encountered  during  the  DFU  Delta  Accident. 
X  This  prograt  uses  energy  height  theory  to  calculate  inert ial  space 
position. 

clear 

clc 

foriat  bank 

fonat  cotpact 

% 

XIHITIRLI2RTI0H  RHD  END  I RONflENTRL  CONSTANTS 

X   The  fo Noting  initialize  constants  and  variables  for  the  lodel. 
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XfllRCRflFT  CONSTANTS 

X  The  fol  loiiing  constants  are  for  the  particular  aircraft  being 

studied. 

dispC     LANDING  APPAOACH  MODEL') 

dispC  ') 

ans»ers0.0; 

dispCType  aircraft: ' ) 

dispC       1.   P3  at  89,5001b') 

dispC   2.  P3  at  111,5001b1) 

dispC   3.  P3  at  89,5001b  3  engine') 

dispC   4.  L1011 ') 

dispC   5.  T4V) 

ansier^inputC  Enter  desired  aircraft  by  number:  '); 

If  ansnerM1, 

[S,K,CDo,T»ax,Ut,CLo,CLalpha,alpha»ax,q,U1,T,theta]=P3LtApp; 

XP3  @  89, 500 lbs,  app  config. 
else  if  anstera=2, 

[S,K,CDo,Tiax,Ut,CLo,CLalpha,alpha»ax,q,U1,T,theta]*P3HvApp; 

XP3  e  1M,  500  lbs,  app  config. 
else  if  ansierss3, 

[S,K,CDo,Tiax,Ut,CLo,CLalpha,alphaiiax,q,U1,T,theta]sP3Eout; 

XP3  with  one  engine  out. 
el  self  ansiers*4, 

[S,K,CDo,Taax,Ut,CLo,CLalpha,alphaiax,q,U1,T,theta]=L10nApp;IL1011 
el  self  ansier"«5, 

[S^^DoJtax^t^CLo.CLalpha^alphaBax^q^UIJ^thetal-HIAppiXHI 
end 

XS«plane  fori  surface  area  in  ftA2 
XK=coeff icient  for  CD  calculation 
XCDo»coeff icient  of  drag 
XTiax=«ax  thrust  available  in  ft-lbs 
XUt "aircraft  teight  in  lbs 

f1«Ut/32.1?4;  X  mass  in  slugs 
XCLo«zero  flOA  lift  coefficient 
XCLalphaH ift  curve  slope 
Xalphaiax-stal I  buffet  alpha  ■/  approach  flaps 
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Xqsnominal  pitch  rate  of  5deg/sec 
XU1=approach  equivalent  airspeed  in  ft/s 
XT«initial  thrust  in  Ibf 


X***npoQPQined  i  n  i  t  i  a  I  va  I  ugs********1*********'''***'''*1''**************1''* 
iters*30;  Xnumber  of  iterations 

deltatB1.0;  Itime  increment  in  sec 

HBcenters10300;  ^distance  fIB  center  is 

from  runway 
xstart=-500;  Xst art ing  distance  from 

MB  center 

g************************************^ 

hstartB(f1Bcenter-xstart)*tan(0.052);  Xstarts  on  a  3deg 

gl ides  I  ope 

dispC  Initial  Input  Section') 

dispC  ') 

d  i  sp ( ' The  f o 1 1 00  i  ng  va I ues  are   preprogrammed : ' ) 

dispC  a.  The  approach  is  based  upon  a  3deg  gl ides  lope. ') 

dispC  b.  The  starting  distance  (distance  from  MB  center)  is:') 

disp(xstart) 

dispC  c.  The  flB  center  is  this  from  the  end  of  the  runmay:') 

disp(flBcenter) 

dispC  d.  The  starting  altitude  is:') 

dlsp(hstart) 

dispC  e.  The  time  step  in  sec:') 

disp(deltat) 

dispC  f.  The  number  of  iterations  is:') 

disp(iters) 

dispC  ') 

TUBTmax/Ut;  Xthrust/meight  ratio 

UL-Ut/S;  liing  loading 

d  i  sp ( ' The  f 0 1 1 om I ng  ua I ues  mere  supp I i  ed/ca I cu I  at  ed  from  a  i  rcra  f t 

constants') 

dispC  a.  Uref  in  ft/s  (1 .35Uso+5kts  for  P3  and  1.3Uso  for  other 

aircraft: ') 

disp(UI) 

dispC  b.  Thrust  to  meight  is:') 
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disp(TU) 

dispC  c.  Uing  loading  is: ') 

disp(UL) 

dispC  ') 

■aneuuer80.0; 

disp( 'Select  the  escape  naneuuer:') 

dispC  1.  Stay  on  glide  slope,  go  tissed  at  200ft  RGL;  NO 

■indshear') 

dispC  2,  Stay  on  glide  slope,  go  missed  at  200ft  RGL;UITH 

•indshear') 

dispC  3.  Escape  with  constant  airspeed') 

dispC  4.  Escape  iith  constant  altitude') 

dispC  5.  Escape  tith  constant  theta') 

dispC  6.  Escape  with  constant  alpha') 

dispC  7.  Escape  lith  tax  alpha') 

dispC  8.  Eiulate  DflL 191  final  few  seconds  using  recorded  theta') 

dispC  9.  Eiulate  DflL 191  final  fei  seconds  using  recorded  theta  and 

winds') 

■aneuversinput('  Enter  selection  nuiber:  '); 

if  ianeuver"5, 

dispCUhat  value  for  escape  theta?') 

t net aa i ■■ i nput ( '  Enter  in  radians:  '); 
el  self  ■aneuuer*"6l 

dispCUhat  value  for  escape  alpha?') 

alphaaii*input('  Enter  in  radians:  '); 
else  if  ianeuver"?, 

alphaaii*alphaiax; 
el  self  ianeuver>B8, 

hstartshstart+50; 
end 

xl-xstart;  linit ial izat ion 

ht-hstart;  linit ial izat ion 

Ug«U1;  Xground  speed  initialization  in  ft/s 

gaiai=-,052;  X  in  it ial  inert ial  flight  path  angle  in 

rad  (-3deg) 
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r o«fO. 002377; 
triggers0.0; 
check1a0.0; 
check2*0.0; 
X 


%a\r   density  in  s lugs/ft A3  @  S.L. 
Xut  i I ity  check  number 
Xut  i I ity  check  number 
Xut  i  I  ity  check  number 


XflLGORITHH 

X  The  actual  algorithm  is  broken  into  several  smaller  divisions  and 

funct  ions: 


starting  values  -  as  it  implies 

loop  -  beginning  of  iterating  process 

■indshear  inputs  -  uses  the  function  'microburst'  to  give 

point  x  and  z  axis  mind  speeds 
predictor  -  first  guess  at  the  solution  of  the  governing  DEs 
corrector  -  corrected  solutions  to  the  governing  DEs  and 

final  point  values, 
output  file  management  -  as  it  implies 
aircraft  control  -  uses  multiple  control  functions  for 

control  inputs  to  theta,  thrust,  and  flOfl 
final  output  -  as  it  implies 


Xcounter 

Xdynamic  pressure 

Xcoeff icient  of  I ift 

Xcoefficient  of  drag 

XROR  required  for  1g  flight 

Xairmass  flight  path  angle 

Xspecific  energy  in  ft 

Xtime  rate  change  in  specific  energy 


X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 

% 

Xstart ing  values: 
Pt-1; 

Q-0.5*rom*Ur2; 
CL-Ut/(Q*S); 
CD«CDo+K*Cr2; 
alpha«(CL-CLo)/CLalpha; 
gamaa-t net a-a I pha; 
Es*Ur2/(2*32.174)+h1; 
Esdot«-12.32; 


y(1)*x1;  Xy(1)«x  position  in  ft 

y(2)-h1;  Xy(2)«aircraft  altitude  in  ft 

y(5)«U1;  Xy(5)*aircraft  airspeed  in  ft/s 

y(6)-gamaa;  Xy(6)*aircraft  airmass  fit,  path  (rad) 

Aircraft (1,:)-[x1, hi, U1 ,-12. 36,theta, alpha, 1,T]j 

Inert  ial (1 , : )>[Ug,0,gamaa,gamal ,0,0, 1 ,Es,Esdot]; 

X  starting  output  files 
% 
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-999, 


X  loop : 

clc 

•hi  le  theta' 

pt°pt+1 ; 

% 

Xmind  shear   inputs: 
x2*Ug*deltat+x1; 
h2-tan(gamai)*(x2-x1)+h1; 

if  maneuver"1 , 
Ux11«0;Ux21«0;Ux12=0;Ux22 

else  if  maneuuer=ss9, 

[UxH ,Uh11 ]=recwinds(pt); 

[Ux22,Uh22>recminds(pt+1 

Ux12-(Ux11+Ux22)/2; 

Ux21=Ux12; 

Uh12=(Uh11+Uh22)/2; 

Uh21«Uh12; 

else 
[Ux11,Uhn]-microburst(h1 
[Ux21,Uh21]-microbur3t(h2 
[Ux12,Uh12]-microburst(h1 
[Ux22,Uh22]«microburst(h2 


XI oop  control 
Xpresent  pt  in  space 


^approximate  horizontal  displacement 
Xapproximate  vertical  displacement 


'0;Uh11«0;Uh21=0;Uh12=0;Uh22=0; 


Xrecorded  minds  from  DRL191 


end 

UindX(pt)-Ux11; 

UindH(pt)-Uhl1; 


,x1);  XThese  four  calls  to  the  function 
,x1)j  Xmicroburst  finds  the  point  mind 
,x2);  Xshear  for  determination  of  the 
,x2);  Xaverage  block  mind  shear. 
XThe  above  is  used  in  conjunction 
Xother  calls  mhen  no  mind  shear  is 
Xdesired. 

XUind  shear  matrices 


dUxdx=((Ux12-Ux11)+(Ux22-Ux21))/(2*(x2-x1));XThese  formulas 
dUxdh-((Ux22-Ux12)+(Ux21-Ux11))/(2*(h2-h1));Xcalculate  the 
dUhdx-((Uht24lh11)+(Uh22-Uh21))/(2*(x2-x1)) ;Xdi fferent iai  change 
dUhdh-((Uh22-Uh12)+(Uh21-Uh11))/(2*(h2-h1));Xin  mind  shear. 


y(3)=Ux11; 

y(4)=Uh11; 
I 


Xy(3)*mind  shear  in  the  x  dir.  ft/s 
Xy(4)smind  shear  in  the  z  dir.  ft/s 
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Xpredictor  and  interned! ate  values 
ydot=ypr  i  ies(y ,  dUxdx,  dlixdh,  dUhdx,  dllhdh,  theta  J,  S,  fl,  0,  CL,  CD, 
alpha^alphamax) ; 

Xca 1 1  to  the  OE  function 

for  ns1 :6, 
yl (n)«y(n)+deltat*ydot(n);     XEuler  first  step  fid  predictor 
yhalf(n)By(n)+deltat/2*ydot(n);  XEuler  half  step  fid  predictor 

end 

alpha«theta-y(6);         XROR 

CL*CLo+alpha*CLalpha;      Xcoefficient  of  lift 

CD*CDo+K*CLA2;  ^coefficient  of  drag 

QB0.5*roi*yhalf(5)A2;      Xdynaiic  pressure 

I 

(corrector  and  nei  value  formulation 
ydot-ypriiesCyhalf.dUxdx.dUxdh.dUhdx^Uhdh^thetaJ.S.rt^^CL.CD, 
alpha,alpha«ax); 

for  n«1 :6, 
y2(n)«yha I f (n)+de I tat/2*ydot (n) ; XR  i  chardson  extrapo I  at  i  on 
y(n)«2*y2(n)-y1(n);  Xcorrector  Scheie. 

end 

x1»y(1);  Xnei  x  posit  in  ft 

h1-y(2);  Xnei  altitude  in  ft  AGL 
Ug«y(5)*cos(y(6))+y(3);        Xground  speed  in  ft/sec 
R0Oy(5)*sin(y(6))+y(4);        XRate  of  Cliib  in  ft/sec 
gaiai«atan(ROC/Ug);  Xnei  a i mass  FP  angle  in  rod 

Es-y (5)*2/64. 348+y (2);         Xnei  specific  energy 
Esdot«y(5)*ydot(5)/32.174+ydot(2)jXnei  trc  of  specific  energy 

alpha=theta-y(6);         Xnei  ROR 

CL*CLo+alpha*CLalpha;      Xnei  coefficient  of  lift 

CD«CDo+K*CLA2;  Xnei  coefficient  of  drag 

Q-0.5*roi*y(5)A2;         Xnei  dynaiic  pressure 

% 
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Xoutput  file  sanagenent : 
R i rcraf t(pt, :)=[y(1),y(2),y(5),ydot(5),theta,alpha,pt,T]; 
Inertia  I  (pt,  :)=[Ug,R0C,y(6),ga»ai ,y(3),y(4),pt,Es,Esdot]; 

hoie 

dispC  pt  x  h 

airspeed  accel  ' ) 

opl -[pt , fl  i  rcraf t (pt , 1 ) , fl  i rcraf t (pt , 2) , A  i rcraf t (pt , 3) , 

flircraft(pt,4)]; 
disp(opl) 

dispC  pt      theta      alpha     ganaa      ganai') 
op2=[pt , fl  i  rcraf t (pt , 5) , R  i rcraf t (pt ,6),  I nert  i  a  I (pt ,  3) , 

InertiaKptj  4)]; 
disp(op2) 

dispC         pt         Ug        ROC        Thrust') 

op3«[pt,lnertial(pt,1),lnertial(pt,2),flircraft(pt,8)]; 

disp(op3) 

dispC         pt         Es         Esdot') 

op5»[pt, lnertial(pt,8),  lnertial(pt,9)]; 

disp(op5) 

dispC        pt         Ux         Uh') 

op4-[pt, Inert ial(pt,5), Inert  ia I (pt ,6)3; 

disp(opl) 

I 

laircraft  control : 

^functions  called  for  theta,  thrust 
land  alpha  control . 
if  taneuwer  <■  2, 
if  y(2)<200, 
trigger*1; 

[theta,T]«nflcontrol(q,deltat,Tiax,theta,T,alpha,alphaiax,y); 
else  if  trigger  ~*   1, 
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[theta^Tl-GPcontnol (q.del tat ^ROC^Taiax^ theta^T.alpha.alphanax.UI , 

hstart,xstart ,y); 
elsei  f  trigger  ss   1 , 

[thetaiT]anRcontnol(q,deltat,Tiax>theta;Tialphaialphaiiiax>y)j 
end 
elsei  f  »aneuver<8, 
if  y(5)  <  (U1-33.33), 

trigger*! ; 
end 
else 

if  y(2)<«1004, 

[thetaJ^heckn-DflLthetaCTiax.thetaJ^UI^y^checkl); 
else 
[theta J]-GPcontrol(q,deltat,ROC Jiax,theta,T,alpha,alphaiax,U1 , 

hstart,xstart,y); 
end 
end 

if  uaneuuer  >«3, 
if  trigger  ■■  1, 
if  taneuuer  ■■  3, 
[theta,T,check1]=CRSesc(q,deltat,Tiax,theta,T,alpha,alphaiiax, 

U1,y,check1) 
el  self  ianeuyer  ■■  4, 

[theta^TfCheckll^RLTescCq^deltat^Tiax^theta^T^alpha^alphaBax, 

y, check 1,ydot) 
else  if  ianeuver  ■■  5, 

[ t beta, T ]-CTHesc(q, de 1 1 at , Tiax, t het a, T, a  I pha , a  I phaiax, thetaa i ■) 
else  if  ianeuver  ■■  6, 

[thetaJl-CflOflescfq^deltatJiax^hetaJ^alpha^lphaaax^alphaaii) 
else  if  ianeuver  ■■  ?, 

[ t het a , T]anflXesc(q, de I  tat , Tiax, thet a, T , a  I pha, a  I phaa  i ■) 
end 
end 
end 

if  pt" iters,  Xnuiber  of  iterations 

theta-999; 
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elseif  y(2)<0,  Xstop  at  ground  iipact 

theta*999; 
else 
end 

end  Xwhile  algorithw  loop 

% 

loutput  files: 
answer* input ('Do  you  want  a  flight  path  plot?  1 =VES  '); 
i  f  answer*=1 , 

plot (Aircraft (: ,1), Aircraft (:  ,2)) 

pause 
answer*0.0; 

end 

answer* input ('Do  you  want  x,  h,  U,  and  airspeed  accel?  1*YES  '); 

if  answer"  1 , 

dispC      pt      x       h       U      airspeed  accel  ') 

op1-[Aircraft(:,7),Aircraft(: , 1), Aircraft (:,2), Aircraft (: ,3), 

Aircraft(:,4)]; 
disp(opl) 
else 
end 

answer* input ('Do  you  want  theta,  alpha,  gawaa,  and  gawai?  1 -YES  '); 

if  answer** 1, 

dispC      pt      theta     alpha     gawaa      gawai') 

op2-[Aircraft(:,7),Aircraft(:l5),Aircraft(:,6),lnertiaI(:,3), 

Inert ial(: :,4)]j 
disp(op2) 
else 
end 

answer* input ('Do  you  want  Ug,  ROC,  and  Thrust?  1-VES  '); 

if  answer** 1, 

dispC        pt        Ug         ROC       Thrust') 

op3«[Aircraft(:,7), Inert ial(:,1), Inert ial(:,2), Aircraft (: ,8)]; 

disp(op3) 
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else 
end 


answer* input( 'Do  you  iant  Es  and  Esdot?  1sYES  "); 

i  f  answer3* 1 , 

dispC         pt         Ed         Eedot') 

op5=[fl  i  rcraf t ( : , 7) , I nert  i  a  I ( : , 8) , I nert  i  a  I ( : , 9) ] ; 

disp(op5) 

else 

end 

answer*  input  ('Do  you  want  Ux  and  I4h?  1-VES  '); 

i  f  answer*"!, 

dispC         pt         Ux         Uh') 

op4*[R  i rcraf t ( : , 7) J nert  i  a  I ( : , 5) , I nert  i  a  I ( : , 6) ] ; 

disp(opl) 

else 

end 
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funct ion  [ydot]Bypriies(y,dllxdx,dUxdh,dUhdx,dUhdh,theta, 

T,S,f1,Q,CL,CD,alpha,aipha»ax) 
This  function  calculates  the  y'  solutions  to  the 
coupled  system  equations: 

y(1 ) ,=xdot=U*cos(gaiaa)+Ux 

y(2)'=hdot=U*sin(ga»aa)+Uh 

y(3),=Uxdot=dUx/dx*xdot*dUx/dh*hdot 

y(4)'=Uhdot=dUh/dx*xdot*dUh/dh*hdot 

y(5) '  =Udot  =  (T*cos(a  I  pha)  )/f1-D/f1-g*s  i  n(gataa)-Uxdot*cos(ga»aa) 

-Uhdot*sin(ga»aa) 
y(6)'«T*sin(alpha)/(l1*U)+L/(l1*U)-g*co8(ga»aa)/U 
+Uxdot*sin(ga»aa)/U-Uhdot*cos(ga»aa)/U 


ydot(1)«y(5)*cos(y(6))+y(3); 

ydot(2)«y(5)*sin(y(6))+y(4); 

ydot (3)=dUxdx*ydot ( 1 )+dUxdh*ydot (2) ; 

ydot (4)-dUhdx*ydot ( 1 )+dUhdh*ydot (2) ; 

if  alpha<alphaaax, 

ydot(5)-T/n*cos(theta-y(6))-(J*S*CD/n-32,174*sin(y(6)) 

-ydot(3)*cos(y(6))-ydot(4)*sin(y(6)); 
ydot(6)-T/(n*y(5))*sin(theta-y(6))+Q*S*CL/(n*y(5)) 

-32.174/y(5)*cos(y(6))+ydot(3)*sin(y(6))/y(5) 
ydot(4)/y(5)*cos(y(6)); 
else 
ydot  (5)«T/f1*cos(a  I phatax)-Q*S*CD/n-32 . 1 74*s I n(y (6) ) 

-ydot(3)*cos(y(6))-ydot(4)*  sin(y(6)); 
ydot(6)-T/(n*y(5))*sin(alph«ax)+Q*S*CL/(n*y(5))-32.17Vy(5) 

*cos(y(6))+ydot(3)»  sin(y(6))/y(5)-ydot(4)/y(5)*cos(y(6)); 
end 
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B.  AIRCRAFT  PERFORMANCE  PARAMETER  FUNCTIONS 


functlonCS^^CDo^Tiax^Ut^CLo^CLalpha^alphaiax^q^UI jT^theta] 

=P3Ltfipp() 
XP3  at  89, 500 lbs,  app  flaps,  gear  down 

S=1 300 j  Xplane  fort  surface  area  in  ftA2 

K*0. 05041;  Xcoefficient  for  CD  calculation 

CDos0.0567;  Xcoefficient  of  drag 

Twax=33400;  Xwax  thrust  avai I  able  in  ft- 1  be 

Ut=89500;  Xaircraft  weight  in  lbs 

CLo=0.800;  Xzero  flQfi  lift  coefficient 

CLalpha*5.73;  Xlift  curve  slope 

alphawaxs0.244;  Xstall  buffet  alpha  •>/  approach  flaps 

q=0.0873;  Xnowinal  pitch  rate  of  5deg/sec 

U1 =236;  Xapproach  equivalent  airspeed  in  ft/s 

T=5000;  Xinitial  thrust  in  Ibf 

theta«-0.013;  Xinitial  water  line  deck  angle  in  rad 


f unct  i  on[S, K , CDo , Twax, Ut , CLo, CLa I pha, a  I phawax, q , U1 , T , theta] 

«P3Hvfipp() 
XP3  at  114,000 lbs ,  app  flaps,  gear  down 

S-1300;  Xplane  fori  surface  area  in  ftA2 

K-0. 05041;  Xcoefficient  for  CD  calculation 

CDo*0.0567;  Xcoefficient  of  drag 

Twax=33400;  Xwax  thrust  available  in  ft- lbs 

Ut-1 14500;  Xaircraft  weight  in  lbs 

CLo-0.800;  Xzero  ROfl  lift  coefficient 

CLa I pha»5 .73;  Xlift  curve  s I  ope 

alphawax=0.244;  Xstall  buffet  alpha  w/  approach  flaps 

q«0.08?3;  Xnowinal  pitch  rate  of  5deg/sec 

U1«262;  Xapproach  equivalent  airspeed  in  ft/s 

T-8000;  Xinitial  thrust  in  Ibf 

theta«-0.013;  Xinitial  water  line  deck  angle  in  rad 
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funct  ion[S,K,CDo,Tnax, 

=L101 1flpp() 
XL1011  at  362,000163, 
S-4578; 
K-0.059; 
CDo-0.108; 
Tiax-126000; 
Ut-362000; 
CLo«0.532; 
CLalpha-4.96; 
alpha»ax«0.3M; 

q=0.0873; 
U1-227; 

W4500; 
theta-0. 1 18; 


Ut,CLo,CLalpha,alphaiax,q,U1 , T,theta] 

fipp  flaps,  gear  dotn 

Xplane  for»  surface  area  in  ft"2 

Xcoefficient  for  CD  calculation 

Xcoefficient  of  drag 

Imax  thrust  ava liable  in  ft- lbs 

laircraft  weight  in  lbs 

Xzero  ROR  lift  coefficient 

XI ift  curve  slope 

Xstall  farning  alpha  •/  approach  flaps 

X(18deg) 
Xnoiinal  pitch  rate  of  5deg/sec 
Xapproach  equivalent  airspeed  in  ft/s 

%(136kts) 
Xinitial  thrust  in  Ibf 
X initial  tater  line  deck  angle  in  rad 


funct  ion[S,K,CDo 
«H4Rpp( 
Xm  at  8280 lbs, 
S=210; 
K-0.0503; 
CDo»0.120; 
Tiax«3023; 
Ut»8260; 
CLo«0.587; 
CLalpha«6.24; 
alpha»ax*0.244; 
q«0.Q873; 
U 1-203; 
T-1250; 
theta*-0.02; 


, Tiax, Ut , CLo, CLa I pha, a  I phaiax, q, U1 , T, theta] 

) 

app  flaps,  gear  dom 

Xplane  fori  surface  area  in  ft*2 

Xcoefficient  for  CD  calculation 

Xcoefficient  of  drag 

Xiax  thrust  available  in  ft-lbs 

Xaircraft  ■eight  in  lbs 

Xzero  R0R  lift  coefficient 

XI ift  curue  slope 

Xstall  buffet  alpha  ■/  approach  flaps 

Xnoiinal  pitch  rate  of  5deg/sec 

Xapproach  equivalent  airspeed  in  ft/s 

Xinitial  thrust  in  Ibf 

Xinitial  later  line  deck  angle  in  rad 
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funct  i on[S, K, CDo , Tnax, Ut , CLo, CLa Ipha.a lpha»ax, q^U1 ^T^theta] 

=P3Eout() 
XP3  at  89, 500 lbs,  app  flaps,  gear  dcrnn,  one  engine  out. 

S«1300;  Xplane  fori  surface  area  in  ft "2 

K-0.050^1 ;  Xcoefficient  for  CD  calculation 

CDo«0. 0630;  Xcoefficient  of  drag  (inc  for  eng,  out) 

T»ax«s33^00*(3/^);  X»ax  thrust  available  in  ft-lbs 

Ut-89500;  Xaircraft  weight  in  lbs 

CLo«0.800;  Xzero  ROfl  lift  coefficient 

CLa I pha=5 . 73;  Xlift  curve  slope 

alphanaxs0.2'H;  Xstall  buffet  alpha  •/  approach  flaps 

q=0 . 0873;  Xnoiinal  pitch  rate  of  5deg/sec 

01=236;  Xapproach  equivalent  airspeed  in  ft/s 

T-5000;  Xinitial  thrust  in  Ibf 

theta*-0,013;  Xinitial  later  line  deck  angle  in  rad 
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C.  FLIGHT  PATH  CONTROL  AND  ESCAPE  MANEUVER  FUNCTIONS 


funct  i  on[t hetaout , Tout ]=GPcontro I (q, de I  tat  j ROC, T»ax, theta, T, 

a  I pha, a  I phaiax , U1 , net art , xst art , y ) 
XThis  control  function  tries  to  naintain  a  3deg  glide  path. 
Tout*T; 

thetaout«theta; 
•antedh«hstart-(abs(xstart-y(1))*0.0524); 

Xdesired  altitude  for  3deg  GP 


hdiff»y(2)-iantedh; 
Udiff=y(5)-U1; 
X   theta  input: 
if  alpha<0, 

thetaout*theta-alphaj 
e lee  if  alpha<alphaiax, 
if  R0O0, 

thetaoutBtheta-q*de I  tat ; 
eleeif  R0C<-25.0, 

thetaout«theta+q*deltat; 
eleeif  hdiff>100, 

thetaout«theta-q/2*deltatj 
eleeif  hdiff<-100, 

thetaout«theta+q/2*deltat; 
end 

if  thetaout>0.524, 
t hetaout -0.524; 
eleeif  thetaout<-0.175, 

thetaout«-0.175; 
end 
elee 

t  he t aout ■ t net  a- ( a  I pha-a I phaiax ) ; 
end 

XThruet  input: 
if  Udlff>8.5, 

Tout-T-0.20*Tt,ax*deltat; 
if  Tout<0, 
Tout*0; 
elee 


Xa I t  difference  (neg  when  low) 
Xairepeed  diff  (neg  nhen  slow) 

KDcheck  for  negative  ROR 

X(2)check  for  aax  available  ROR 
X(3)check  for  ROC  >  Ofpt 

% i  f  so  pitch  dwn  q*t  ise 
%(4)check  for  ROC  <  1500fp» 

%\  f  eo  pitch  up  q*t iie 
X(5)check  for  alt  MOO  diff 

Xif  so  pitch  dun  1/2*q*tiie 
X(6)check  for  alt  < 1 00  diff 

Xif  eo  pitch  up  l/2*q*tise 

Know  check  for  theta  lie  it  of 
%   -lOdeg  to  +30  deg 


Xcheck  speed  diff 
Xif  diff  >  5kte  above  app  speed 
Xreduce  thrust  by  20%  of  T*ax 
Xbut  not  less  than  0. 
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end 
elseif  Udiff<-8.5, 

Tout «T+0.25*T»ax*de! tat; 

if  Tout>Tiax, 
Tout*T»ax; 

else 

end 
end 


XI  f  diff  <  5kte  be  lot  app  speed 
Xinc  thrust  by  25X  of  Tnax 
Xbut  not  tone  than  Tnax. 


KOcheck  for  lax  ava liable  flOfl 
X(2) check  for  lOdeg  deck  angle 


f unct  i  on[thetaout , Tout ]=nflcont  ro I (q, de I  tat , T»ax , t  het a , T ,  a  I pha , 

alpha»ax,y) 
XThis  control  function  initiates  a  noma  I  aissed  approach. 
X  Max  poter  is  added  and  theta  increased  to  lOdeg.  flax  flOfl 
Xis  applied  «hen  required. 
Tout«T; 

thetaout-theta; 
X  theta  input: 
if  alpha<alphaiax, 

if  thetaout>»0.175, 
thetaout -0.175; 

else 

thetaout-theta+q/2*deltat;    X(3)if  theta  not  lOdeg,  pitch  to 
lOdeg 

end 
else 

thetaout-theta-(alpha-alphaiax); 
end 

XThrust  input: 
if  Tout<Tiax, 

Tout-T+Tiax/2*deltat; 

if  Tout>Tiax, 
Tout-Tiax; 

end 
else 
end 
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functiontthetaout Jout,check1]"CflSe5c(q>deltat Jiax,theta,T, 

alpha,alpha»ax,U1 ,y,check1 ) 
Xflircraft  escape  maneuver  consisting  of  constant  airspeed. 
thetaout«theta; 
if  check  1—0, 
if  y(5)<U1, 

thetaout=0; 
else 
thetaout«q*deltat; 
checkl-1; 
end 
else 

if  y(5)>«(Ul+8.4)> 

thetaout"theta+q*deltat; 
elseif  y(5)<-(U1-8.4) 

t  het  aout  «t  het  a-q*de 1 1  at ; 
end 
end 
if  alpha>alphaiax, 

thetaoutstheta-(alpha-aiphaiax); 
end 
if  KTiax, 

Tout-T+Tiax*0.2*deltat; 
else 

TouWiax; 
end 
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funct  ion[thetaout, Tout, check  1  ]BCflLTesc(q,deltat >T«ax,theta,T , 

alpha,alpha»ax,y,check1 , ydot) 
Xflircraft  control  using  constant  altitude  escape  laneuuer. 
thetaout*theta; 
i  f  check  1"0; 

check1sy(2); 
end 
if  y(2)<check1-20, 

if  ydot(2)<»0, 
thetaout«theta+q*de I  tat ; 

end 
elseif  y(2)>check1+20, 

if  ydot(2)>-0, 
thetaout«theta-q*deltat; 

end 
end 
if  alpha>alphaiax, 

thetaout-theta-(alpha-alphaiax); 
end 
if  T<Ttax, 

Tout-T+Tiax*0.2*deltat; 
else 

Tout-Tiax; 
end 
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f unct  i  on[ thetaout ,  Tout ]«CTHesc(q, de 1 1  at , Tiax, theta, T , 

alpha,alpha»ax,thetaai») 
{Constant  theta  escape  ianeuver. 
if  theta<thetaaii, 

thetaout«theta+q*deltat; 
else 

thetaout*thetaain; 
end 
if  alpha>alphaiax, 

thetaoutstheta-(alpha-alpha»ax); 
end 
if  KTiax, 

Tout  «T+T»ax*0 . 2*de 1 1  at ; 
else 

TouWiax; 
end 


f unct  i  on[ thetaout , Tout >CflOflesc(q, de I  tat , Tiax, theta, 1, 

alpha,alpha8>ax,alphaai») 
{Constant  alpha  escape  ianeuver , 
if  alpha<alphaaii, 
thetaout "theta* (a  I phaa  i i-a I pha) ; 
if  thetaout>(theta+q*deltat), 

t  het  aout  «t  het  a+q*de 1 1  at  j 
end 
else 
thetaout-theta-(alpha-alphaali); 
if  thetaout<(theta-q*deltat), 
thetaout «t het a-q*de I  tat  j 
end,  end 
if  alpha>aiphaiaxj 

t  het aout "t het  a- ( a  I pha-a I phaiax ) ; 
end 
if  T<Tiax, 

Tout=T+Tiax*0.2*deltat; 
else ,   ToutBTtax;  end 
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functiontthetaout  Jout]snflXesc(q,deltat  J»ax,thetaJ,alpha,alphaaii) 
Xflaxitui  alpha  escape  laneuuer. 
if  alpha<alphaai», 
t  het  aout = t het  a+ ( a  I phaa  i i-a  i  pha ) ; 
if  thetaout>(theta+q*deltat), 

thetaout«theta+q*deltat; 
end 
else 
t het aout «t het a-(a I pha-a I phaa  i ») ; 
if  thetaout<(theta-q*deltat), 

thetaout*theta-q*deltat ; 
end 
end 
i  f  T<T«ax, 

Tout«T+Tiax*0.2*deltat; 
else 

TouWtax; 
end 
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funct  ion[thetaout,Tout,check1out]sDflLtheta(Ttax>theta>T,U1 ,y,check1 ) 

X  This  control  function  tries  to  emulate  DRL 191 

Hast  48  sec  of  fl  ight. 

if  check1"«0,  Istarting  point 

check1B32; 
end 

%   theta  input : (theta.Xponer) 
DRL191-C0. 069813170071; 
0.0628318530  71;  0.054105206871; 

72;  0.047123889873;  0.038397243573; 

71;  0.022689280271;  0.038397243569, 

69;  0.054105206869;  0.054105206870, 

68;  0.123918376867;  0.132645023167, 

67;  0.185004900767;  0.191986217767, 

67;  0.200712863969;  0.200712863971 

72;  0.253072741585;  0.267035375589, 

92;  0.274016692592;  0.267035375595, 

92;  0.237364778292;  0.237364778292, 

87;  0.253072741587;  0.223402144292, 

93;  0.062631853098;  0.031415926596] 

99;  -0.005235987100;  -0.068067840100; 

100;  -0.130899693100;  -0.020943951100; 

100;  0.0925024503100;  0.0314159265100; 

100;  0.0541052068100;  0.0314159265100]; 
thetaout-DAL191 (checkl+1 , 1 ); 
Tout-DRL191  (checkl+1, 2)/100*Twx; 
check lout -check 1+1; 


0.0471238898 
0.0314159265 
0.0541052068 
0.0785398163 
0.1623156204 
0.2007128639 
0.2234021442 
0.2740166925 
0.2600540585 
0.2530727415 
0.1553343034 
0.0226892802 
-0.144862327 
0.0698131700 
-0.005235987 
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APPENDIX  D 
TAKEOFF  ALGORITHMS 

Contained  in  this  appendix  are  the  MatLab®  programs  utilized  in 
calculating  the  aircraft  response  given  a  windshear  and  a  particular 
escape  maneuver.  They  are  listed  in  three  sections.  The  first  section 
contains  the  driver  routine.  The  second  and  third  sections  contain 
the  specific  aircraft  parameter  functions  and  escape  maneuver 
functions  respectively.  The  windshear  function  is  listed  In  Appendix 
A.    The  differential  equation  solving  routine  is  listed  in  Appendix  C. 

A.  DRIVER 

X   Takeoff  Profile  Rode  I   * 
g***************************** 

X   This  prograt  ill  I  calculate  the  Inert lal  flight  path  angle  for 
X  a  given  aircraft.  It  utilizes  the  function  '■icroburst*  ehich 
X  models  the  find  ehear  encountered  during  the  DFU  Delta  Accident. 
X  This  prograi  uees  energy  height  theory  to  calculate  inert lal  epace 
position. 

clear 

clc 

fonat  bank 

fonat  compact 

% 

I  INITIALIZATION  AND  ENUIAONflENTAL  CONSTANTS 

X  The  fol losing  initialize  constants  and  variables  for  the  iodel. 
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XAIACAAFT  CONSTANTS 

X  The  fo Noting  constants  are  for  the  particular  aircraft  being 

studied. 

dispC     TAKEOFF  HODEL ' ) 

dispC  ') 

ansiers0.0; 

dispCType  aircraft : ') 

dispC   1.  P3  at  90,0001b') 

dispC   2.  P3  at  120,0001b') 

dispC   3.  P3  at  135,0001b') 

dispC   4.  L1011') 

dispC   5.  H4') 

dispC   6.  P3  at  120,0001b,  Ur-HO') 

answer" inputC  Enter  desired  aircraft  by  number:  '); 

i  f  ans«er«"1, 

[S,K,CDo,Ut,CLo,CLalpha,alphaiax,q,U1,U2,T,thetacliib]-P3LtToff; 

theta*0.087;  Xinit ial izat ion 

else  if  ans«er==2, 

[S,K,CDo,Ut,CLo,CLalpha,alphaiax,q,U1,U2,T,thetacliib]*P3HvToff; 

theta*0.130;  Xinit ial izat ion 

el  self  anster—3, 

[S,K,CDo,Ut,CLo,CLalpha,alphaaax,q,U1,U2J,thetacliib]«P3HuHuToff; 

theta-0.133;  Xinitial ization 

else  if  ansier*«1, 

[S,IC,CDoJUt,CLo,CLalpha,alphaiax,q,U1,U2,T,thetacliib]-L1011Toff;XL1 
011 

theta*0.2;  Xinitial ization 

else  if  ansier«e5, 

[S,K,CDo,Ut,CLo,CLalpha,alphaiax,q,U1,U2,T,thetaclimb]-T44Toff; 
XT14 

theta-0.200;  Xinit ial izat ion 

else  if  ansier«6, 

[S,K,CDo,Ut,CLo,CLalpha,alphawax,q,U1,U2,T,thetacliib]-aodP3HuTo; 

theta-0.130;  Xinitial ization 

end 
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XSsplane  fop»  surface  area  in  ftA2 
XK«coeff icient  for  CD  calculation 
XCDo=coeff icient  of  drag 
XUt "aircraft  weight  in  lbs 

l1»Ut/32.174;  X  boss  in  slugs 
XCLo=zero  flOfl  lift  coefficient 
XCLa I pha= I i  ft  curve  slope 
Xa lpha»ax=sta 1 1  buffet  alpha  •/  approach  flaps 
Xq»no»inal  pitch  rate  of  5deg/sec 
XU1 "rotate  speed  in  ft/s 
XU2"takeoff  safety  airspeed  In  ft/s 
XT-thrust  in  Ibf 


|***pr oar aied  initial  values*************************************** 
iters-65;  Xnuiber  of  iterations 

deltat80.5;  Xtine  increment  in  sec 

xstarts-12Q0;  Xstarting  distance  froi 

riB  center 
hstarts2.0;  Xstarts  at  rotate 

dispC     Initial   Input  Section1) 

dispC    ') 

dispCThe  fo Noting  values  are  preprogramed: ' ) 

dispC  a.  The  flight  path  is  based  on  tax  poter.') 

dispC  b.  The  starting  distance  (distance  frot  riB  center)  is:') 

disp(xstart) 

dispC  c.  The  MB  center  is  this  froi  the  end  of  the  runiay:') 

disp(abs(xstart)) 

dispC  d.  The  tite  step  in  sec:') 

disp(deltat) 

dispC  e.  The  nuiber  of  iterations  is:') 

disp(iters) 

dispC  ') 

TU-T/Ut;  Xthrust/ieight  ratio 

UL-Ut/S;  Wing  loading 

dispCThe  following  values  tere  supplied/calculated  froi  aircraft 
constants') 
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spC  a.  Takeoff  Safety  Speed:') 

ep(U2) 

spC  b.   Thrust  to  weight    is:') 

sp(TU) 

sp( '  c.  Uing  loading  is: ') 

sp(UL) 

spC  ') 


aaneuver=0.0; 

di spC Select  the  escape  aaneuver:') 

dispC  1.  Takeoff  flight  path;  HO  aindshear') 

dispC  2.  Takeoff  flight  path;  UlTH  aindshear') 

dispC  3.  Escape  with  constant  airspeed') 

dispC  4.  Escape  with  constant  altitude') 

dispC  5.  Escape  with  constant  theta') 

dispC  6.  Escape  with  constant  alpha') 

ianeuuer-inputC  Enter  selection  nutber:   '); 

if  aaneuverM5, 
dispCUhat  value  for  escape  theta?') 
thetaaii«input( '  Enter  in  radians:  '); 

else  if  ianeuver»*6, 
dispCUhat  value  for  escape  alpha?') 
alphaaii-inputC  Enter  in  radians:  '); 

end 


x1«xstart; 
h1«hstart; 


Unit  ial  ization 
linitial ization 


Ug-UI; 

gaiai =0.001 ; 
rad  (Ideg) 
roi-0. 002377; 
trigger-0.0; 
check 1-0.0; 
check2«0.0; 
X 

IflLGORITHfl 


{ground  speed  initialization  in  ft/s 
X  initial  inertia  I  flight  path  angle  in 

Xair  density  in  slugs/ftA3  f  S.L. 
Xuti I ity  check  nuaber 
Xuti I ity  check  nuaber 
Xuti I ity  check  nuaber 
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X  The 

actual  a  Igor 

it  hi 

is  broken  into  several 

snal ler  divisions  and 

funct 

ons: 

starting  va 

ues  - 

as  it  implies 

loop  -  beginning 

of  iterat  ing  process 

■indshear  inputs 

-  uses  the  funct  ion  'ni 

croburst '  to  give 

point  x  and  z  axis  wi 

nd  speeds 

predictor  - 

first 

guess  at  the  solut  ion 

of  the  governing  DEs 

corrector  - 

corrected  solutions  to  the  governing  DEs  and 

final 

point  values. 

output  file 

Ranagement  -  as  it  implies 

aircraft  control 

-  uses  mult iple  control 

functions  for 

control  inputs  to  theta,  thrust,  and  flOfl 

final  outpul 

.  -  as 

it  itplies 

1st  art ing  values: 
pt-1; 

Q-0.5*roi*U1A2; 
CL«Ut/(Q*S); 
CD«CDo+K*CLA2; 
a  I pha*(CL-CLo)/CLa I pha; 
gaiaa-theta-alpha; 
Es-Ur2/(2*32.174)+h1; 
Esdot«-12.32j 


Xcounter 

Xdynaiic  pressure 

Xcoeff icient  of  I  i  ft. 

Xcoeff icient  of  drag 

IflOfl  required  for  1g  flight 

lair mass  flight  path  angle 

Xspeclfic  energy  in  ft 

Xtiie  rate  change  in  specific  energy 


Xy(1)*x  position  in  ft 
Xy(2)-aircraft  altitude  in  ft 
Xy(5)aaircraft  airspeed  in  ft/s 
Xy(6)«aircraft  airnass  flight  path  in 


y(1)-x1; 

y(2)-h1; 

y(5)-U1; 

y(6)«g«aa; 

rad 

Aircraft (1,  :)-[x1,h1,U1 ,-12. 36, theta, alpha, 1  J]; 

I nert I  a  I ( 1 , s )-[Ug , 0 , gaiaa , gana i , 0, 0 , 1 , Es , Esdot ] ; 

X  starting  output  files 
% 

Xloop: 

clc 

■hile  theta~«999, 

pt-pt+lj 

% 


Xloop  control 
Xpresent  pt  in  space 
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Xiind  shear  inputs: 
x2sUg*deltat+x1;  Xapproximate  horizontal  displacement 

h2Btan(gamai)*(x2-x1 )+h1 ;    Xapproximate  vert ical  displacement 

if  maneuver~«1 , 
[Uxll ,Uh1 1]*iicroburst(h1 ,x1);  XThese  four  calls  to  the  function 
[Ux21,Uh21]amicroburst(h2,x1);  Xmicroburst  finds  the  point  find 
[Uxl2,Uh12]«microburst(h1,x2);  Xshear  for  determination  of  the 
[Ux22,Uh22]smicroburst(h2,x2);  Xaverage  block  wind  shear, 
else 
Ux11-0;Ux21-0;Ux12-0;Ux22-0;Uh11-0;Uh21=0;Uh12«0;Uh22-0; 

XThe  above  is  used  in  conjunction 
Xother  calls  when  no  mind  shear  is 
Xdesired. 
end 

UindX(pt)-Ux11;  XUind  shear  matrices 

UlndH(pt)-Uh11; 

dUxdx-((Ux12-Ux11)+(Ux22-Ux21))/(2*(x2-x1));XThese  formulas 
dUxdh-((Ux22-Ux12)+(Ux21-Ux11))/(2*(h2-h1));Icalculate  the 
dUhdx-( (Uhl 2-Uhl 1 )+(Uh22-Uh21 ) )/(2*(x2-x1 ) ) ; Xd  i  f f erent  i  a  I  change 
dUhdh«((Uh22-Uh12)+(Uh21-Uh11))/(2*(h2-h1));Xin  mind  shear. 

y(3)*Ux11;  Xy(3)*mind  shear  in  the  x  direction 

ft/s 

y(4)-UM1;  Xy(4)-mind  shear  in  the  z  direction 

ft/s 

I 

Xpredictor  and  intermediate  values 
ydot-yprimesCy.dUxdx.dUxdh^dUhdx^Uhdh^thetaJ^S.n.Q.CL.CD.alpha.alp 
hamax); 

Xcal I  to  the  DE  function 

for  n*1 :6, 
y1(n)By(n)+deltat*ydot(n);     XEuler  first  step  fid  predictor 
yhalf(n)*y(n)+deltat/2*ydot(n);  XEuler  half  step  fid  predictor 

end 
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alpha=theta-y(6);         XROR 

CL=CLo+alpha*CLalpha;       Xcoefficient  of  lift 

CD=CDo+K*CLA2j  Xcoefficient  of  drag 

Q«0.5*roi*yhalf(5)A2;      Xdynaiic  pressure 

% 

Xcorrector  and  nee  value  formulation 
ydot«Bypriiee(yhalfJdUxdxJdUxdh,dUhdxJdUhdh(thetaJTJS,n,0;CL,CD,alpha 
,alphaiax); 

for  ns1 :6, 
y2(n)=yha I f (n)+de I tat/2*ydot (n) ; XR  i  chardson  extrapo I  at  i  on 
y(n)*2*y2(n)-y1(n);  Xcorrecter  Scheie. 

end 

x1»y(1);  Xnei  x  posit  in  ft 

M«y(2);  Xnei  altitude  in  ft  RGL 
Ug*y(5)*cos(y(6))+y(3);         Xground  speed  in  ft/sec 
R0C-y(5)*sin(y(6))+y(4)j        XRate  of  CI  lib  in  ft/sec 
gaiai-atan(ROC/Ug);  Xnei  alriass  FP  angle  in  rad 

Es*y(5K2/64.348+y(2)j         Xnei  specific  energy 
Esdot*y(5)*ydot(5)/32.174*ydot(2)j  Xnei  trc  of  specific  energy 

alpha-theta-y(6);         Xnei  ROR 

CL-CLo+alpha*CLalpha;      Xnei  coefficient  of  lift 

CD-CDo+K*CLA2;  Xnei  coefficient  of  drag 

Q*0.5*roi*y(5)*2;         Xnei  dynaiic  pressure 

I 

Xoutput  file  lanageient: 
flircraft(pt,:)«[y(1),y(2),y(5),ydot(5),theta,alpha,pt,T]; 
Inert ial(pt,:)«[Ug, ROC, y(6),gaiai,y(3),y(4),pt,Es,Esdot]; 

hoie 

dispC  pt  x  h  U 

airspeed  accel  ') 

op1-[pt aircraft (ptl1)ifllrcraft(ptl2)iRircraft(pt ,3) , Aircraft (pt ,4)] 
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disp(opl) 

dispC  pt  theta        alpha        ganaa 

ganai ' ) 

op2*[pt ,  fl  i  rcraf t (pt , 5) , fli rcraf t (pt , 6) ,\ nert  i  a  I (pt , 3) , I nert  i  a  I (pt , 

t)]; 

disp(op2) 

dispC  pt  Ug  ROC  Thrust') 

op3=[pt , I nert  i  a  I (pt , 1 ), I nert  i  a  I (pt , 2) , fl  i rcraf t (pt , 8) ] ; 
disp(op3) 

dispC         pt         Es         Esdot') 

op5«[pt,  Inert ial (pt, 8),  Inert  ial  (pt  ,9)]; 

disp(op5) 

dispC         pt         Ux  Uh1) 

op4s[pt, Inert ial (pt, 5), Inert ial (pt,6)]; 

dispCopD 

% 

Ial rcraf t  control : 

^functions  called  for  theta,  thrust 
If  y(2)>50,  land  alpha  control, 

if  ROC  <-0, 
trigger«1; 
elseif  y(5XU2, 

trigger«1; 
end 
elseif  y(5)<U1, 

trigger=1; 
end 

if  trigger«"0, 
t net a«theta+q*0 . 6*de I  tat ; 
if  thetaHhetacliib, 
theta*thetacl iib; 
end 
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else 

i  f  «aneuuer«3, 

[theta^checkll^ToffCRSCq^deltat^theta^alpha^alphaiax^UZ.y.checkl) 
else  if  taneuver-M, 

[theta^checkl ]sToffCflLT(q,deltat,thetQ,alpha,alphaB»axJy,check1 ,ydot) 
elsei  f  ianeuueras5, 

[ thet a]=Tof fCTH(q, de 1 1 at , t het a, a  I pha, a  I phaiax,t net aa it) 
else  if  »aneuveras6, 

[theta]eToffCflOfl(q,deltatJtheta,alpha/alphao>ax,alphQai») 
end 
end 

if  pt" Iters,  Inunber  of  iterations 

theta-999; 
elseif  pt>5, 

if  y(2)<0.0, 
theto-999; 

end 
end 

end  Xthlie  algorithm  loop 

I 

lout put  files: 
anseer ■ I nput ( ' Do  you  lant  a  flight  path  plot?  1 ■ VES  ' ) ; 
If  ansierM1, 

plot(flircraft(:  ,1),  Aircraft  (:,  2)) 

pause 
ansier-O.O; 

end 

anseer- input ('Do  you  iant  x,  h,  U,  and  airspeed  accel?  1 -VES  '); 

if  anster«1, 

dispC  pt  x  h  U 

airspeed  accel  ') 

op1-[fllrcraft(:l?)ifllrcraft(:J)lfllrcraft(:l2)lflircraft(:l3),fllrcraf 

id  Mi 
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disp(opl) 

else 

end 

answer* input  ( 'Do  you  want  theta,  alpha,  ga»aa,  and  ga&ai?  1  «=VES  '); 

if  ans»er==1 , 

dispC  pt  theta        alpha        ganaa 

gaiai ') 

op2*[fl  I  rcraf  t  ( : , 7) , fl  i  rcraf  t  ( : ,  5) ,  fl  i  rcraf t  ( : ,  6) ,  I  nert  i  a  I  ( : ,  3)  ,  I  nert  I  a 

disp(op2) 

else 

end 

ans*ers input ('Do  you  »ant  Ug,  ROC,  and  Thrust?  1=VES  '); 

if  ans»er==1 , 

dispC         pt         Ug         ROC       Thrust') 

op3s[fl  i  rcraf t (:,?),! nert  i  a  I  ( : ,  1 ) ,  I nert  i  a  I ( : , 2) , fl  i rcraf t ( : , 8) ] ; 

disp(op3) 

else 

end 

ansier" input ('Do  you  iant  Es  and  Esdot?  1"VES  '); 

if  ansier«"1, 

dispC         pt         Es        Esdot') 

op5«[flircraft(:,7), Inert ial(:,8), Inert ial(:, 9)]; 

disp(op5) 

else 

end 

answer* input ('Do  you  tant  Ux  and  Uh?  1«YES  '); 

if  anster*"1, 

dispC        pt         Ux         Uh') 

op4«[fl  i  rcraf t (:,?),! nert  i  a  I ( : ,5) , I nert  i  a  I ( : , 6) ] ; 

disp(op4) 

else 

end 
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B.  AIRCRAFT  PERFORMANCE  PARAMETER  FUNCTIONS 


function[S,K,CDo,Ut,CLo, 

«P3LtToff() 
XP3  at  90,000168,  app  fl 
S»1300; 
KO.  05041  ; 
CDo«0.0551; 
Ut -90000; 
CLo-1.000; 
CLalpha«5.73; 
alpha»ax«0,244; 
q-0.0873; 
U1-204; 
U2-220; 
T-33400; 
thetacl inb-0. 175; 


CLa I pha, a  I phatax , q, U1 , U2, T , t het ac I i  «b] 

aps,  gear  up 

Xplane  fori  surface  area  in  ftA2 

Xcoefficient  for  CD  calculation 

Xcoefficient  of  drag 

Xaircraft  weight  in  I  be 

Xzero  flOfl  lift  coefficient 

XI ift  curve  elope 

Xetall  buffet  alpha  ■/  approach  flape 

Xnonlnal  pitch  rate  of  5deg/eec 

XLiftoff  epeed  (121kte) 

Xtakeoff  eafety  airepeed  in  ft/e 

Xthruet  in  Ibf 

Xtater  line  deck  angle  in  rad 


functlon[S,K,CDo,Ut,CLo 

•P3HvToff() 
XP3  at  120 ,000 1  be,  app 
S-1300; 
K-0.05041; 
CDo-0.0551; 
Ut-120000; 
CLo*1.Q00; 
CLalpha«5.73; 
alphaiax«0.244; 
q-0.0673; 
U1-214; 
U2-227; 
T-33400; 
thetacl iib«0. 175; 


,CLalpha,alphaiax,q,U1 ,U2,T, thetacl lib] 

flape,  gear  up 

Xplane  fori  surface  area  in  ftA2 

Xcoefficient  for  CO  calculation 

Xcoefficient  of  drag 

Xaircraft  weight  in  lbs 

Xzero  flOfl  I i  ft  coe f f  i  c  i  ent 

X 1 1 f t  curve  e I  ope 

Xetall  buffet  alpha  •/  approach  flape 

Xnoiinal  pitch  rate  of  5deg/eec 

XLiftoff  epeed  (127kte) 

Xtakeoff  eafety  airepeed  in  ft/e 

Xthruet  in  Ibf 

Xvater  line  deck  angle  in  rad 
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f  unct  1  on  [  S ,  K ,  CDo ,  Ut ,  CLo ,  CLa  I  pha ,  a  I  phanax  ^,111,112,1,1  het  ac  I  i  «b  ] 

«P3HvHvToff() 
XP3  at  135, 000 lbs,  app 
S=1300; 


K-=0. 05041 

CDo«0.055' 

Ut-135000, 

CLo-1.000; 

CLalpha»5.73; 

alphaiax*0.244; 

qs0.0873; 

U1-229:, 

U2«=239; 

T«33400j 

thetacliwb=0,1745j 


flaps,  gear  up 

Xplane  fori  surface  area  in  ft*2 

Xcoefficient  for  CD  calculation 

Xcoefficient  of  drag 

Xaircraft  weight  in  lbs 

Xzero  ROfl  lift  coefficient 

XI ift  curue  slope 

Xstall  buffet  alpha  w/  approach  flaps 

Xnoninal  pitch  rate  of  5deg/sec 

XL  if  toff  speed  (J36kts). 

Xtakeoff  safety  airspeed  in  ft/s 

Xthrust  in  Ibf 

Xwater  line  deck  angle  in  rad 


function[S,K,CDo,Ut,CLo,CLalpha,alphawax,q,U1,U2,T,thetacl iwb] 

«iodP3HvTo() 
XP3  at  120, 000 lbs,  app 
S-1300; 
K-0. 05041 
CDo-0.0551; 


Ut- 120000] 

CLo-I.OOOj 

CLalpha-5.73; 

alphaiax-0,244; 

q-0.0673; 

U1-242; 

U2-254; 

T-33400; 

thetacliwb-0.175; 


flaps,  gear  up.  Increased  rotate  speed. 
Xplane  fori  surface  area   in  ftA2 
Xcoefficient  for  CD  calculation 
Xcoefficient  of  drag 
Xaircraft  weight  in  lbs 
Xzero  ROfl  lift  coefficient 
XI ift  curue  slope 

Xstall  buffet  alpha  w/  approach  flaps 
Xnowinal  pitch  rate  of  5deg/sec 
XL  if  toff  speed  (145kts) 
Xtakeoff  safety  airspeed  in  ft/s 
Xthrust  in  Ibf 
Xwater  line  deck  angle  in  rad 
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funct  ion[S,K,CDo,Ut,CLo,CLalpha,alphciiax,q,U1 ,U2,T,thetacl i 

-LIOIIToffO 
XL  1 0il  at  462, 000 lbs,  flpp  flaps,  gear  up 


ib] 


S-4578; 
K-0.059; 

CDo-0.098; 

Ut»462000; 

CLo-0.532; 

CLalpha-4,96; 

alphaiax-0.3H; 

qs0.0873; 
U1=238; 
U2-255; 
T« 126000; 

thetaclitb-0.209; 


Xplane  fori  surface  area  In  ftA2 

Xcoefficient  for  CD  calculation 

Xcoefficient  of  drag 

Xaircraft  weight  in  lbs 

Xzero  R0R  lift  coefficient 

XI ift  curve  slope 

Xstall  warning  alpha  •/  approach  flaps 

X(18deg) 
Xnowinal  pitch  rate  of  5deg/sec 
Xrotate  speed  in  ft/sec  (Mlkts) 
Xtakeoff  safety  airspeed  ft/s  (151 kts) 
X  thrust  in  Ibf 
Xinitial  deck  angle  rad(12deg) 


funct  i  on[S,  K,  CDo,  lit ,  CLo,  CLa 

■TMToffO 
XT41  at  781 7 lbs,  app  flaps, 
S=210; 
K-0.040; 
CDo*0.010; 
Ut-7817; 
CLo-0.0524; 
CLalpha-6.24; 
alphawax"0.227; 
q-0.0873; 
U1-152; 
U2-180; 
T-3023; 
thetacl iwb=0.262; 


pha,alphawax,q,U1,U2,T,thetacl iwb] 

gear  down 

Xplane  fori  surface  area  in  ft*2 
Xcoefficient  for  CD  calculation 
Xcoefficient  of  drag 
Xaircraft  weight  in  lbs 
Xzero  flOfl  lift  coefficient 
XI ift  curve  slope 
Xstall  buffet  alpha  w/flaps  up 
Xnowinal  pitch  rate  of  5deg/sec 
Xliftoff  airspeed  in  ft/s(90kts) 
Xtakeoff  safety  speed  ft/s  (107kts) 
Xinitial  thrust  in  Ibf 
Xinitial  deck  angle  in  rod  (15  deg) 
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C.  FLIGHT  PATH  CONTROL  AND  ESCAPE  MANEUVER  FUNCTIONS 

funct  i  on[ t hetaout , checkl ] 

=Tof  f CRS(q, de I  tat , t beta, a  I pha, a  I phaiax, U2, y , checkl ) 
XRIrcnaft  escape  ianeuver  consitlng  of  constant  airspeed. 
thetaout»theta; 
if  checkl —0, 
If  y(5XU2, 
thetaout*theta-q*deltat; 
if  thetaout<0.0, 
thetaouta0.0; 
end 
else 
thetaout-q*deltat; 
checkl "1; 
end 
else 

If  y(5)>«(U2+8.4), 

thetaoutBtheta*q*deftat; 
elself  y(5)<«(U2-8.4) 

thetaout*theta-q*deltat; 
end 
end 
If  a  I pha >a I phaiax , 

thetaout-theta-(alpha-alphatax); 
end 
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funct  iontt hetaout, check  1] 

■ToffCRLT(q,deltat,theta,alpha,alphaiiax,y,check1,ydot) 
Xfl  i  rcra  ft  cont  no  I  us  i  ng  const  ant  a  1 1  i  t  ude  escape  aaneuuer . 
thetaout-theta; 
if  check1*s0; 

check1«y(2)+20; 
end 
if  y(2)<check1-10, 

if  ydot(2)<=0, 
t hetaout «theta+q*de I  tat ; 

end 
elseif  y(2)>checkl+IO, 

if  ydot(2)>-0, 
thetaout«theta-q*deltat; 

end 
end 
if  thetaouKO, 

thetaout^O; 
end 
if  alpha>alphatax, 

thetaout*theta-(alpha-alphaiax); 
end 


funct ion[t hetaout >ToffCTH(q, del  tat, thet a, alpha, a  I phaiax,thetaai ■) 
^Constant  thet a  escape  ■aneuuer. 
if  theta<thetaaia, 

t hetaout "theta+q*de I  tat ; 
else 

thetaout8thetaaia; 
end 
if  alpha>alphaiax, 

thetaout-theta-(alpha-alphaiax); 
end 
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f unct  i  on[ thetaout  3-Tof f CflOfi(q, de I  tat  ^  theta ^  a  I pha , a  I pha»ax, a  I phaa  in) 

XConstant  alpha  escape  ■aneuuer. 
if  alphcKalphaaii, 
t net aout»t heta+ (a  I phaa it-alpha); 
if  thetaout>(theta+q*deltat), 

thetaoutBtheta+q*deltat; 
end 
else 

thetaout=theta-(a I pha-a Iphaa  i ») ; 
if  thetaout<(theta-q*deltat), 

thetaout stheta-q*deltatj 
end 
end 
if  alpha>alphaaax, 

thetaout-theta-(alpha-alphanax); 
end 
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APPENDIX  E 
ON-STATION  ALGORITHMS 

Contained  in  this  appendix  are  the  MatLab®  programs  utilized  in 
calculating  the  response  of  a  P-3  encountering  a  windshear  while 
operating  close  to  the  surface.  They  are  listed  in  three  sections.  The 
first  section  contains  the  driver  routine.  The  second  and  third 
sections  contain  the  aircraft  parameter  function  and  flight  path 
control  function  respectively.  The  windshear  function  is  listed  in 
Appendix  A.  The  differential  equation  solving  routine  is  listed  in 
Appendix  C. 

A.  DRIVER 

%     P3  Loiter  Profile  flodel       * 

g***************************** 

clear 

clc 

fonat  bank 

fornat  coepact 

% 

XINITIALIZATION  AND  END  I  AQNflENTAL  CONSTANTS 

X  The  fo Noting  initialize  conetante  and  variables  for  the  node  I. 

XAIRCRAFT  CONSTANTS 

%    The  following  conetante  are  for  the  particular  aircraft  being 

studied. 

dispC    P3  ON-STATION') 
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diapC  ') 

check1"0.0; 

U"0.0; 

answer  =1 .0; 

check13input('  Enter  on-station  loiter  alt itude(ft):   '); 

Ueinput('  Enter  loiter  airspeed  in  knots:   '); 

U*U*1 .687;  Xconvert  knots  to  ft/sec 

answer* input ('Enter  0  for  no  windshear:  '); 

[S,K,CdotUt  fClo>Cla\pbQ,a\phQ9Qx,qll  ttheta]B?3on9tQ(\))'t 

Xinitial izat ion 

XS"plane  fori  surface  area  In  ft "2 

XK«coeff icient  for  CD  calculation 

XCDoBcoeff icient  of  drag 

XUt "aircraft  weight  in  lbs 

M*Ut/32. 1 74;  X  wass  in  slugs 

XCLo-zero  ROfl  lift  coefficient 

XCLalpha"! ift  curve  slope 

Xalphawax"stal I  buffet  alpha  ■/  approach  flaps 

Xq»notinal  pitch  rate  of  5deg/sec 

XU1 "rotate  speed  in  ft/s 

XU2"takeoff  safety  airspeed  in  ft/s 

XT"thrust  in  Ibf 

% 

|***prograied  initial  ua I ues*************************************** 
iters*H5;  Inuwber  of  iterations 

del  tat "0.25;  Xtiae  increment  in  sec 

xstart--5000;  Istarting  distance  froi 

flB  center 

dispC  Initial  Input  Section') 

dispC  ') 

dispCThe  following  values  are  preprogramed: ') 

dispC  a.  The  flight  path  is  based  on  constant  power.') 

dispC  b.  The  starting  distance  (distance  froi  FIB  center)  is:') 

disp(xstart) 

dispC  c.  The  tiie  step  in  sec:') 
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disp(deltat) 

dispC  d.  The  number  of  iterations  is:1) 

disp(iters) 

dispC  ') 

TU-T/Ut; 

UL-Ut/S; 

dispCThe  following  values  were  supplied/calculated  from  aircraft 

constants') 

dispC  a.  Loiter  flir  Speed:') 

disp(U) 

dispC  b.  Thrust  to  weight  is:') 

disp(TU) 


Ithrust/weight  ratio 
Jwing  loading 


dispC  c. 

disp(UL) 

dispC  ') 

xl-xstart; 

h 1 -check 1; 

U1-U; 

Ug-U; 

gatai -0.001; 

rad  (Ideg) 

rom-0. 002377; 

trigger-0.0; 

check2-0.0; 


Uing  loading  is: ') 


Unit  ial  izat  ion 
linit ial izat ion 
I  in  it ial izat ion 

Xground  speed  initialization  in  ft/s 
Xi nit ial  inert ial  flight  path  angle  in 

lair  density  in  slugs/ftA3  •  S.L. 
Xut i I ity  check  number 
Xut 1 1 ity  check  number 


RLGORITHn 

The  actual  algorithm  is  broken  into  several  smaller  divisions  and 
unctions: 

starting  values  -  as  it  implies 

loop  -  beginning  of  iterating  process 

■indshear  inputs  -  uses  the  function  'microburst'  to  give 

point  x  and  z  axis  mind  speeds 
predictor  -  first  guess  at  the  solution  of  the  governing  OEs 
corrector  -  corrected  solutions  to  the  governing  DEs  and 

final  point  values, 
output  file  management  -  as  it  implies 
aircraft  control  -  uses  multiple  control  functions  for 

control  Inputs  to  theta,  thrust,  and  flOfl 
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final  output  -  as  it  implies 


Xstart  ing  values: 
pt-1; 

Q«0.5*ro»*Ur2; 
CL-Ut/(Q»S); 
CD-=CDo+K*CL"2; 
alpha*(CL-CI_o)/CLalpha; 
ga»aa«theta-alpha; 
Es«Ur2/(2*32.1?4)+h1; 
"hsdot=U; 


Xcounter 

Xdynamic  pressure 
^coefficient  of  lift 
^coefficient  of  drag 
XROR  required  for  1g  flight 
Xairiass  flight  path  angle 
Xspecific  energy  in  ft 
Xtiie  rate  change  in  spec* i "tic  energy 


y(1)-x1; 

y(2)*h1j 

y(5)-U1; 

y(6)«ganaa; 

rad 

Aircraft (1 , :)«[x1,h1,U1,0,theta, alpha,  1,T]j 

Inert  ial (1 , : )«[Ug,0,ga»aa,gaiai ,0,0, 1 ,Es,Esdot ]; 

X  starting  output  files 


Xy(1 )=x  posit  ion  in  ft 
Xy(2)Baircraft  altitude  in  ft 
Xy(5)«aircraft  a 
Xy(6)«aircraft 


rspeed  in  ft/s 

a i mass  f  I  ight  path  in 


Xloop: 
clc 

■hi  le  theta' 
pt-pt+1; 
I 


'999, 


Xpresent  pt  in  space 


Xiind  shear  inputs: 
x2-Ug*deltat+x1; 
h2-tan(ganai)*(x2-x1)+h1; 


Xapproxiiate  horizontal  displacement 
Xapproxiiate  vertical  displacement 


If  ansfer'^O, 
[Ux11,Uh11]-iicroburst(h1,xD 
[Ux21,Uh21]*iicroburst(h2,xD 
[Ux12,Uh12>ticroburst(h1,x2) 
[Ux22,Uh22]-iicroburst(h2,x2) 
if  pt<-8, 


XThese  four  calls  to  the  function 
Xticroburst  finds  the  point  find 
Xshear  for  determination  of  the 
Xaverage  block  find  shear. 


XThese  next  if  logic  conands 
Ux11«Ux11*pt/8;Uh11-Uh11*pt/8;  Xcontrol  the  entry  and  exit 
Ux 1 2-Ux 1 2*pt /8 ; Uh 1 2-Uh 1 2*p t /8 ;  U\ ndshear  effects. 
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elseif  y(1)>sabs(xstart), 
if  check2S!S0, 
check2«pt; 
end 

if  (8-pt+check2)>«0, 
Ux21«Ux21*(8-(pt-check2))/8;Uh21»Uh21*(8-(pt-check2))/8; 
Ux22»Ux22*(8-(pt-check2))/8;Uh22«Uh22*(8-(pt-check2))/8; 
else 

Uxl1«0;Ux21»0;Ux12»0;Ux22«0;Uh11-0;Uh21=0;Uh12«0;Uh22=0; 
end 
end 
else 
Ux11»0;Ux21-0;Ux12-0;Ux22«0;Uh11=0;Uh21«0;Uh12-0;Uh22-0; 

XThe  above  is  used  in  conjunction 
Xother  calls  then  no  lind  shear  is 
Xdesired. 
end   Xansier 

UindXCpt)-Ux11;  XUind  shear  latrices 

UindH(pt)-Uh11; 

dUxdx«((Ux12-Ux11)+(Ux22-Ux21))/(2*(x2-x1));XThese  fonulas 
dUxdh-((Ux22-Ux12)^(Ux21-Ux11))/(2*(h2-h1));Xcalculate  the 
dUhdx-((Uh12-Uh11)+(Uh22-Uh21))/(2*(x2-x1)); differential  change 
dUhdh-((Uh22-Uh12)+(Uh21-Uh1l))/(2*(h2-h1));«in  lind  shear. 

y(3)«Ux11;  Xy(3)«iind  shear  in  the  x  direction 

ft/s 

y(4)-Uh11;  Xy(1)-iind  shear  in  the  z  direction 

ft/s 

% 

Xpredictor  and  intermediate  values 
ydot-ypriiesCy.dUxdx.dUxdh^Uhdx.dUhdh.thetaJ.S^.Q^L.CD, 
alpha,alphaiax);    Xcall  to  the  DE  function 
for  n«1 :6, 
y1(n)«y(n)+deltat*ydot(n);     XEuler  first  step  fid  predictor 
yhalf(n)-y(n)+deltat/2*ydot(n);  XEuler  half  step  fid  predictor 
end 
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alpha-theta-y(6)j         XROR 

CL=CLo+alpha*CLalpha;       Xcoefficient  of  lift 

CD-CDo+K*CLA2;  Xcoefficient  of  drag 

Q«0.5*roe*yhalf(5)A2;       Xdynatic  pressure 

% 

Xcorrector  and  new  ualue  formulation 

ydot»yprites(yhalf,dUxdx,dUxdh,dUhdx,dUhdh,theta,T,S,f1,Q> 
CL , CD, a  I pha , a  I phaiax) ; 

for  ns1 :6, 
y2(n),!yha  I  f  (n)+de  I  tat/2*ydot  (n);%Ri  chardson  extrapo  I  at  i  on 
y(n)»2*y2(n)-y1(n);  Xcorrector  scheme. 

end 

x1«y(1);  Xnei  x  posit  in  ft 

h1-y(2);  Xnei  altitude  in  ft  RGL 

Ug«y(5)*cos(y(6))+y(3)j  Xground  speed  in  ft/sec 

R0C-y(5)*sin(y(6))+y(4);  XRate  of  Cliib  in  ft/sec 

gaaaiaatan(ROC/Ug);  Xnei  a i mass  FP  angle  in  rad 

Es«y(5)A2/64.348+y(2);  Xnei  specific  energy 
Esdot»y(5)*ydot(5)/32.1?4+ydot(2);  Xnei  trc  of  specific  energy 

alpha-theta-y(6);         Xnei  ROR 

CL-CLo+alpha*CLalpha;       Xnei  coefficient  of  lift 

CD-CDo+K*CLA2j  Xnei  coefficient  of  drag 

Q«0.5*roi*y(5)A2;         Xnei  dynanic  pressure 

% 

Xoutput  file  lanageient: 
Aircraft (pt, :)-[y(1),y(2),y(5),ydot (5), theta, alpha, pt ,T]; 
I nertial(pt,:)«[Ug, ROC, y(6),gaiai,y(3), y(1), pt,Es,Esdot]; 
hoie 

dispC  pt  x  h  U 

airspeed  accel     ') 
opWpt,fllrcraft(pt ,1), Aircraft (pt ,2), Aircraft (pt, 3), 

flircraft(pt,4)]; 
disp(opl) 
dispC  pt  theta  alpha  gaaaa 

gaiai ') 
op2«[pt,flircraft(pt,5),Rircraft(pt,6), Inert ial(pt,3), 
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InertiaKpt,  4)]; 
disp(op2) 

dispC         pt         Ug        ROC         Thrust') 
op3-[pt , I nert  i  a  I (pt , 1 ) ,  I nert  i  a  I (pt , 2) , fl I rcraf t (pt , 8) ] ; 
disp(op3) 

dispC         pt         Es         Esdot') 
op5e[pt, Inert  ial (pt ,8), Inert  ial (pt ,9)]; 
disp(op5) 

dispC         pt         Ux  Uh') 

op4«[pt, Inert  ial (pt ,5)^ Inert ial (pt, 6)]; 

dlsp(op4) 

% 

Xa i rcraf t  control : 

[theta,T]»CflLT(q;deltat,theta,alpha,alpha»ax>yJcheck1, 

ydot,T,U) 

if  pt" iters,  Xnu«ber  of  iterations 

theta*999; 

elseif  y(1)>5000, 

theta«999; 

end 

end  Xihlle  a  Igor  it hi  loop 

X 

Xoutput  files: 
anster- input ('Do  you  iant  a  flight  path  plot?  1-YES  '); 
If  ansierM1, 

plot(fllrcpoft(:  J), Aircraft  (:  ,2)) 

pause 
ansierB0.0; 
end 

answer- input ('Do  you  iant  x,  h,  U,  and  airspeed  accel?  1-VES  '); 
if  ans»er=*1, 

dispC  pt  x  h  U 

airspeed  accel  ') 
op1-[fllrcraft(: ,7), Aircraft (:,1)lfllrcraft(:,2),fllrcraft(:, 3), 

Alrcraft(!,4)]j 
disp(opl) 
else 
end 
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answer" input ('Do  you  want  theta,  alpha,  ganaa,  and  gamai?  1-VES  '); 

if  answer"  1, 

dispC  pt  theta        alpha        ganaa 

gaial ' ) 

op2B[R  i  rcraf t ( : , 7) , fl  i rcraf t ( : , 5) , fl  i rcraf t ( : , 6) , I nert  i  a  I ( : , 3) , 

Inert ial(:,4)3j 
disp(op2) 
else 
end 

answer- input ('Do  you  want  Ug,  ROC,  and  Thruet?  1-VES  '); 
if  answer— 1, 

dispC         pt         Ug         ROC       Thrust') 
op3-[flircraft ( :,7), Inert ial(:,1), Inertia  I (:, 2) ,Rircraft (: ,8)]; 
diep(op3) 
else 
end 

answer = input ('Do  you  want  Es  and  Esdot?  1-VES  '); 
if  answer-- 1, 

dispC        pt        Es        Esdot') 
op5-[fllrcraft(:,7), Inert lal(: ,8), Inert ial(:,9)l; 
disp(op5) 
else 
end 

answer- input ('Do  you  want  Ux  and  Uh?  1-VES  '); 
if  answer "1, 

dispC        pt         Ux         Uh') 
op4-[R  i rcraf t ( : , 7) ,  I nert  i  a  I ( : , 5) , I nert  ial ( : ,6)] ; 
disp(opl) 
else 
end 
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B.  AIRCRAFT  PERFORMANCE  PARAMETER  FUNCTION 

functiontS^^CDcUt^Lo^Lalpha^alphaiaXiqJ^thetal'Paonstadl) 
XP3  at  120,000163, flaps  and  gear  up. 

S=1 300;  Xplane  for*  surface  area  in  ft*2 

K«0. 05041 j  ^coefficient  for  CD  calculat ion 

CDo=0.0213;  Xcoeff Iclent  of  drag 

Ut-120000;  Xaircraft  weight  in  lbs 

CLo-0,350;  Xzero  flOfl  lift  coefficient 

CLalpha»6.25;  IN  ft  curve  slope 

alpha»ax-0.209;  Xstall  buffet  alpha  */  flaps  up 

q-0.0873;  Inoninal  pitch  rate  of  5deg/sec 

T- . 5*0 . 00237?*ir2*S*(CDo+K*(Ut /( . 5* . 002377*(T2*S) ) A2) ; 

Xthrust  in  Ibf 
theta«(Ut/(.5*.002377*U"2*S)-CLo)/CLalphaj 

Xiater  I ine  deck 
angle  in  rad 
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C.  FLIGHT  PATH  CONTROL  FUNCTION 

f unct I on[thetaout , Tout ]=CflLT(q , de I  tat , theta, a  I pha ,  a  I phanax, y , 

check1,ydot,T,U) 
XRircraft  control  using  constant  altitude  escape  maneuver, 
thetaoutBtheta; 
ToutBT; 
if  y(2)<check1-5, 

If  ydot(2)<=0J 
thetaout=theta+q/2.5*deltat; 

end 
elseif  y(2)>check1+5, 

if  ydot(2)>»0, 
thetaout«theta-q/2.5*deltat; 

end 
end 
If  abe(ydot(2)*deltat)>abs(y(2)-checkl), 

thetaoufthetaout-sign(ydot(2))*q/1*deltat; 
end 
if  alpha>alphaiax, 

thetaout"theta-(alpha-alphaaax); 
end 
If  y(5)<»(U-67.5), 

Tout*T+deltat/2*33400; 
end 
if  Tout>33400, 

Tout-33400; 
end 
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APPENDIX  F 
APPROACH  TO  LANDING  ENCOUNTER  GRAPHIC  DATA 

This  appendix  contains  the  calculated  performance  of  different 
aircraft  and  weight  combinations  upon  encountering  a  microburst 
windshear  during  an  approach  to  landing.  Each  figure  is  for  a 
particular  aircraft  performing  a  specified  escape  maneuver.  Each 
figure  contains  four  graphs  depicting  altitude,  theta,  alpha,  airspeed, 
specific  energy,  and  time.  The  abscissa  In  all  graphs  is  x  (distance 
from  microburst  center). 
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Figure  Fl.    P-3  at  89,5001bs  gross  weight  performing  a 
constant  airspeed  escape. 
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Figure  F2.    P-3  at  89,5001bs  gross  weight  performing  a 
constant  altitude  escape. 
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Figure  F3.    P-3  at  89,5001bs  gross  weight  performing  a 
constant  5°  theta  escape. 
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Figure  F4.    P-3  at  89,5001bs  gross  weight  performing  a 
constant  8«  theta  escape. 
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Figure  F5.    P-3  at  89,5001bs  gross  weight  performing  a 
constant  10°  theta  escape. 
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Figure  F6.    P-3  at  89,5001bs  gross  weight  performing  a 
constant  15°  theta  escape. 
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Figure  F7.    P-3  at  89,500 lbs  gross  weight  performing  a 
constant  12  unit  AOA  escape. 
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Figure  F8.    P-3  at  89,5001bs  gross  weight  performing  a 
constant  15  unit  AOA  escape. 
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Figure  F9.    P-3  at  89,5001bs  gross  weight  performing  a 
constant  20  unit  AOA  escape. 
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Figure  F10.    P-3  at  114,0001bs  gross  weight  performing  a 
constant  airspeed  escape. 
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Figure  Fll.    P-3  at  114,0001bs  gross  weight  performing  a 
constant  altitude  escape. 
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Figure  F12.    P-3  at  114,0001bs  gross  weight  performing  a 
constant  5°  theta  escape. 
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Figure  F13.    P-3  at  114,0001bs  gross  weight  performing  a 
constant  8°  theta  escape. 
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Figure  F14.    P-3  at  114,0001bs  gross  weight  performing  a 
constant  10°  theta  escape. 
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Figure  F15.    P-3  at  114,0001bs  gross  weight  performing  a 
constant  15"  theta  escape. 
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Figure  F16.    P-3  at  11 4,000 lbs  gross  weight  performing  a 
constant  12  unit  AOA  escape. 
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Figure  F17.    P-3  at  114,0001bs  gross  weight  performing  a 
constant  15  unit  AOA  escape. 
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Figure  F18.    P-3  at  114,0001bs  gross  weight  performing  a 
constant  20  unit  AOA  escape. 


202 


0  -\ — i — i — i — i — i— i — i — j — i — i — " — i — " — i — ■ — i — i — i — ■ — j — ■ — r 
-500  0   500  1000  1500  2000  2500  3000  3500  4000  4500  5000  5500 


20 
15 
10 


,|_        | —»■""■?     alpjia [ 


C7> 

*  S  -I 


...■»■■ 


0--== 


■?•  - V ?-v.-*—r 

-♦4iwta* »       i \- 


,rt 


i/:. 


..» — - ♦•. 


-o-i— i — i — i — i — i — i — i — i — i — i — i — i — i — i  »  i — < — i — i — i — i — i — i- 

-500  0   500  1000  1500  2000  2500  3000  3500  4000  4500  5000  5500 


ioo  H — ■ — i — ■ — j — ■ — j — i — j — i — j — « — j — « — j — ■ — j — ■ — i — i — r 

-500  0   500  1000  1500  2000  2500  3000  3500  4000  4500  5000  5500 


-500  0   500  1000  1500  2000  2500  3000  3500  4000  4500  5000  5500 

x(ft) 


Figure  F19.    Three-engine  heavy  airline  transport 
performing  a  constant  airspeed  escape. 
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Figure  F20.    Three-engine  heavy  airline  transport 
performing  a  constant  altitude  escape. 
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Figure  F21.    Three-engine  heavy  airline  transport 
performing  a  constant  5°  theta  escape. 
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Figure  F22.    Three-engine  heavy  airline  transport 
performing  a  constant  10°  theta  escape. 
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Figure  F23.    Three-engine  heavy  airline  transport 
performing  a  constant  15°  theta  escape. 
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Figure  F24.    T-44  performing  a  constant  airspeed  escape. 


208 


600 
500  - 

^  400" 

£  300 
*  200 

100  H 


■ ♦ *-• 


— * -0— 


♦ 


-1 — I — 1 — I — 1 — I — 1 — I — I — I — i — I — I — I — I — f— I — I — I — j- 


-I — "- 


-500  0   500  1 000  1 500  2000  2500  3000  3500  4000  4500  5000  5500 


20 
15 
10 

04 


j    / 


th«ta 


».SJL«  •■.!". 


••■■■• 
~ ;--*••■' ♦ r--» -g- 

•r"  j       '.  i      .■/ 

"I *  ;  fil     '  : 


r.j, ^ i : -,»■»■■*',    ,i 

...! ?""  I 

mi.  m.^K.'.'.'.    '.   ' 

•••?••■••■?     alpjva 


-5  1 — " — j — i — i — i — i — i — i — ■ — j — i — j — * — j — ' — i — i — i — " — i — i — r 
-500     0      500    1 000  1 500  2000  2500  3000  3500  4000  4500  5000  5500 


120 


"»— I — ' — l — ' — l — ' — I — ' — I — ' — I — ' — l — ' — l — "—I — ■ — I — ■ — r 
-500  0   500  1 000  1 500  2000  2500  3000  3500  4000  4500  5000  5500 


1100 


1000 


$     900 


800  -i — i — j — i — j — i — |— « — I — ' — I — « — I — ' — I — « — I — « — P 

-500  0   500  1 000  1 500  2000  2500  3000  3500  4000  4500  5000  5500 

x(ft) 


Figure  F25.    T-44  performing  a  constant  altitude  escape. 
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Figure  F26.    T-44  performing  a  constant  5°  theta  escape. 
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Figure  F27.    T-44  performing  a  constant  10"  theta  escape. 
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Figure  F28.    T-44  performing  a  constant  15"  theta  escape. 
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Figure  F29.    T-44  performing  a  constant  10  unit  AOA 
escape. 
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Figure  F30.    T-44  performing  a  constant  16  unit  AOA 
escape. 
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Figure  F31.    T-44  performing  a  constant  25  unit  AOA 
escape. 
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Figure  F32.    P-3  at  89,5001bs  gross  weight  with  three 
engines  operating  performing  a  constant  5°  theta  escape. 
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Figure  F33.    P-3  at  89,5001bs  gross  weight  with  three 
engines  operating  performing  a  constant  10"  theta  escape. 
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Figure  F34.    P-3  at  89,5001bs  gross  weight  with  three 
engines  operating  performing  a  constant  15°  theta  escape. 
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APPENDIX  G 
TAKEOFF  ENCOUNTER  GRAPHIC  DATA 

This  appendix  contains  the  calculated  performance  of  different 
aircraft  and  weight  combinations  upon  encountering  a  microburst 
windshear  at  or  immediately  after  lift  off.  Each  figure  is  for  a 
particular  aircraft  performing  a  specified  escape  maneuver.  Each 
figure  contains  four  graphs  depicting  altitude,  theta,  alpha,  airspeed, 
specific  energy,  and  time.  The  abscissa  in  all  graphs  Is  x  (distance 
from  microburst  center). 
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Figure  Gl.    P-3  at  90,0001bs  gross  weight  performing  a 
constant  airspeed  maneuver  during  takeoff. 
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Figure  G2.    P-3  at  90,0001bs  gross  weight  performing  a 
constant  10°  theta  maneuver  during  takeoff. 
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Figure  G3.    P-3  at  90,0001bs  gross  weight  performing  a 
constant  12  unit  AOA  maneuver  during  takeoff. 
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Figure  G4.    P-3  at  90,0001bs  gross  weight  performing  a 
constant  15  unit  AOA  maneuver  during  takeoff. 
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Figure  G5.    P-3  at  120,0001bs  gross  weight  performing  a 
constant  10o  theta  maneuver  during  takeoff. 
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Figure  G6.    P-3  at  120,0001bs  gross  weight  performing  a 
constant  15°  theta  maneuver  during  takeoff. 
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Figure  G7.    P-3  at  120,0001bs  gross  weight  with  an  increased 
rotate  speed  followed  by  a  constant  10°  theta  climb. 
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APPENDIX  H 
P-3  FLIGHT  SIMULATOR  DATA 
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FUEL 
TANK 

1 

2 

3 

4 

5 

FUEL 
LOAD 

3937 

3937 

3937 

3935 

996 

ENGINE        1      2     3 
SHP           483   "484   499 
TIT           549   549   551 
FUEL  FLOU     849   849   855 

4 

597 
564 
893 

BRAKES 
PILOT 
COPILOT 

FLIGHT 

IAS    115 

RUNUAY-  100  MOFFETT   NA5  RUY  32R 

ALT   490 
HDG   265.9 

FLIGHT  TIMER 


FLIGHT  CONDITIONS  PAGE 

00:00:00  MET  TIMER 

CONF I GURAT I ON/COND I T I ONS 


GROSS  UEIGHT 

89960 

C.G. 

20.00 

FLAP  POSITION 

19.  1 

GEAR  POSITION 

1.0 

FLIGHT/AERO 

PITCH  ANGLE 

7.3 

ANGLE  OF  ATTACK 

6.4 

HEADING  ANGLE 

283.2 

PITCH  UELOCITY  ( D/S  ) 

-0. 156 

ROLL  UELOCITY  (D/S) 

-0. 164 

YAU  UELOCITY  (D/S) 

-0.273 

NORTH-SOUTH  UELOCITY 

111.71 

EAST-UEST  UELOCITY 

-160.32 

UERTICAL  UELOCITY 

-2.99 

LONGITUDINAL  ACCEL 

8. 1216 

LATERAL  ACCEL 

0.0117 

UERTICAL  ACCEL  (G'S) 

-0.9739 

CONTROL  LOADING 

ELEUATOR  POSITION 

-5.43 

COLUMN  FORCE 

16.06 

RUDDER  POSITION 

-4.34 

PEDAL  FORCE 

-11.06 

AILERON  POSITION 

-0.30 

UHEEL  FORCE 

-2.20 

ENGINES 

TOTAL  THRUST 

5024 

THROTTLE  ANGLE 

46.4 

ENGINE  S.H.P. 

460 

UEIGHT  AND  BALANC 

I XX  INERTIA  (/  1024) 

725 

122  INERTIA  (/  1024) 

1549 

PRESSURE  ALTITUDE 
CALIBRATED  AIRSPD 
EQUIUALENT  AIRSPD 
TRUE  AIRSPD  ( F/S  ) 
MACH  NUMBER 

BANK  ANGLE 
SIDESLIP 

RATE  OF  CLIMB  (FPM) 
PITCH  ACCELERATION 
ROLL  ACCELERATION 
YAU  ACCELERATION 
NORTH-SOUTH  ACCEL 
EAST-UEST  ACCELERATION 
UERTICAL  ACCELERATION 
TOTAL  PITCHING  MOMENT 
TOTAL  ROLLING  MOMENT 
TOTAL  YAUING  MOMENT 

ELEUATOR  TRIM  TAB 
COLUMN  POSITION 
RUDDER  TRIM  TAB 
PEDAL  POSITION 
AILERON  TRIM  TAB 
UHEEL  POSITION 

THRUST  COEFFICIENT 

LATERAL  T.C.  . 

ENGINE  T. I.T. 

E 

IYY  INERTIA  (/  1024) 

CROSS  PRODUCT/ INERT  I A 


00:03:58 

524.0 
114.7 
114.72 
195.41 
0.  18 

-1.4 
-1.4 
181 
0.0004 
-0.0024 
-0.0023 
-0.373 
-0. 138 
0.579 
336 
-1664 
-3575 

12.69 

3.24 

3.79 

-0.47 
0.42 

-1.  16 

0.09 
0.08 
549 

862 
43240 


NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRI 
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FUEL   FUEL 

ENGINE 

1 

2 

3     4 

BRAKES 

TANK   LOAD 

SHP 

485 

486 

496   600 

PILOT 

1    3929 

TIT 

550 

550 

551    564 

COPIL 

2   3929 

3   3929 

FUEL  FLOU 

848 

848 

853   893 

FLIGH1 

" 

IAS 

108 

4    3927 

ALT 

490 

5    996 

RUNU 

NAS  RU*i 

HDG 

256.8 

FLIGHT 

CONDITIONS 

PAGE 

FLIGHT  TIMER 

00: 00: 00 

MET  TIMER              00:04:27 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

89928 

PRESSURE  ALTITUDE 

524.0 

C.G. 

20.00 

CALIBRATED  AIRSPD 

108.6 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

108.59 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 
MACH  NUMBER 

185.00 
0.  17 

FLIGHT/A 

PITCH  ANGLE 

11.9 

BANK  ANGLE 

-0.3 

ANGLE  OF  ATTACK 

7.8 

SIDESLIP 

-1.0 

HEADING  ANGLE 

273.8 

RATE  OF  CLIMB  ( FPM  ) 

804 

PITCH  UELOCITY  (  D/S ) 

0.  195 

PITCH  ACCELERATION 

-0. 0002 

ROLL  UELOCITY  (D/S) 

0.086 

ROLL  ACCELERATION 

0. 0003 

YAU  UELOCITY  (D/S) 

-0.023 

YAU  ACCELERATION 

-0. 0007 

NORTH-SOUTH  UELOCITY 

80.  19 

NORTH-SOUTH  ACCEL 

-0.056 

EAST-UEST  UELOCITY 

-166.20 

EAST-UEST  ACCELERATION 

0.  108 

UERTICAL  UELOCITY 

-13.38 

UERTICAL  ACCELERATION 

0.455 

LONGITUDINAL  ACCEL 

0. 1999 

TOTAL  PITCHING  MOMENT 

-176 

LATERAL  ACCEL 

0.0846 

TOTAL  ROLLING  MOMENT 

256 

UERTICAL  ACCEL  (G'S) 

-0.9636 

TOTAL  YAUING  MOMENT 

-1092 

CONTROL  LOADING 

ELEUATOR  POSITION 

-6.40 

ELEUATOR  TRIM  TAB 

12.69 

COLUMN  FORCE 

16.87 

COLUMN  POSITION 

2.76 

RUDDER  POSITION 

-4.98 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-14.28 

PEDAL  POSITION 

-0.55 

AILERON  POSITION 

0.44 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

0.03 

UHEEL  POSITION 

1.46 

ENGINE 

TOTAL  THRUST 

5232 

THRUST  COEFFICIENT 

0.  10 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C.  . 

0.  10 

ENGINE  S.H.P. 

461 

ENGINE  T. I.T. 

550 

UEIGHT  AND 

BALANCE 

IXX  INERTIA  (/  1024) 

725 

IYY  INERTIA  (/  1024) 

862 

IZZ  INERTIA  (/  10 

24) 

1549 

CROSS  PRODUCT/ INERT  I A 

43230 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL   FUEL 

ENGINE 

1 

2 

3     4 

BRAKES 

TANK   LOAD 

SHP 

486 

487 

505   604 

PILOT 

1    3923 

2   3923 

TIT 

549 

549 

551    564 

COPILOT 

3   3923 

FUEL  FLOU 

850 

850 

856   894 

FLIGHT 

IAS 

117 

4   3921 

ALT 

490 

5    996 

RUNU 

NAS  RU\ 

HDG 

256.2 

FLIGHT 

CONDITIONS 

PAGE 

FLIGHT  TIMER 

00:00: 00 

MET  TIMER              00:04:49 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

89904 

PRESSURE  ALTITUDE 

524.0 

C.G. 

20.00 

CALIBRATED  AIRSPD 

117.4 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

117.37 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  ( F/S  ) 
MACH  NUMBER 

200.03 
0.  18 

FLIGHT/A 

PITCH  ANGLE 

2.4 

BANK  ANGLE 

-0.2 

ANGLE  OF  ATTACK 

6.4 

SIDESLIP 

-0.6 

HEADING  ANGLE 

273.3 

RATE  OF  CLIMB  ( FPM  ) 

-821 

PITCH  UELOCITY  (D/S) 

-0.016 

PITCH  ACCELERATION 

-0. 0000 

ROLL  UELOCITY  (D/S) 

0.000 

ROLL  ACCELERATION 

-0. 0008 

YAU  UELOCITY  (D/S) 

0.023 

YAU  ACCELERATION 

-0. 0000 

NORTH-SOUTH  UELOCITY 

86.  11 

NORTH-SOUTH  ACCEL 

-0.025 

EAST-UEST  UELOCITY 

-180.05 

EAST-UEST  ACCELERATION 

0.210 

UERTICAL  UELOCITY 

13.63 

UERTICAL  ACCELERATION 

-0.482 

LONGITUDINAL  ACCEL 

0.0376 

TOTAL  PITCHING  MOMENT 

-176 

LATERAL  ACCEL 

0.0046 

TOTAL  ROLLING  MOMENT 

-640 

UERTICAL  ACCEL  (G'S) 

-1.0149 

TOTAL  YAUING  MOMENT 

-16 

CONTROL  LOADING 

ELEUATOR  POSITION 

-5.43 

ELEUATOR  TRIM  TAB 

12.69 

COLUMN  FORCE 

15.94 

COLUMN  POSITION 

3.25 

RUDDER  POSITION 

-3.51 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-8.97 

PEDAL  POSITION 

-0.38 

AILERON  POSITION 

0.  19 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.52 

UHEEL  POSITION 

0.52 

ENGINE 

TOTAL  THRUST 

4848 

THRUST  COEFFICIENT 

0.08 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C.  . 

0.08 

ENGINE  S.H.P. 

463 

ENGINE  T. I.T. 

549 

UEIGHT  AND 

BALANCE 

I XX  INERTIA  (/  1024) 

725 

IYY  INERTIA  (/  1024) 

861 

IZZ  INERTIA  (/  1024) 

1549 

CROSS  PRODUCT/ INERT  I A 

43222 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRU 
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FUEL 

FUEL 

ENGINE        1     2     3 

4 

BRAKES 

TANK 

LOAD 

SHP           483   1484   499 

597 

PILOT 

1 

3983 

2 

3983 

TIT          550   550   552 

564 

COPILOT 

3 

3983 

FUEL  FLOU     849   849   855 

893 

FLIGHT 

IAS    114 

^ 

3983 

ALT   490 

5 

997 

RUNUAY-  100  MOFFETT  NAS  RUY  32R 

HDG   321.6 

FLIGHT  TIMER 


FLIGHT  CONDITIONS  PAGE 

00:00:00  MET  TIMER 

CONF I GURAT I ON/COND I T I ONS 


GROSS  UEIGHT 

90144 

C.G. 

24.65 

FLAP  POSITION 

19.  1 

GEAR  POSITION 

1.0 

FLIGHT/AERO 

PITCH  ANGLE 

7.7 

ANGLE  OF  ATTACK 

6.6 

HEADING  ANGLE 

338.8 

PITCH  UELOCITY  (D/S) 

0.250 

ROLL  UELOCITY  (D/S) 

-0. 141 

YAU  UELOCITY  (D/S) 

-0.359 

NORTH-SOUTH  UELOCITY 

195.01 

EAST-UEST  UELOCITY 

5.56 

UERTICAL  UELOCITY 

-4.21 

LONGITUDINAL  ACCEL 

0. 1284 

LATERAL  ACCEL 

-0.0010 

UERTICAL  ACCEL  (G'S) 

-0.9914 

CONTROL  LOADING 

ELEUATOR  POSITION 

-3.96 

COLUMN  FORCE 

8.94 

RUDDER  POSITION 

-3.  16 

PEDAL  FORCE 

-5.06 

AILERON  POSITION 

0.68 

UHEEL  FORCE 

-0.72 

ENGINES 

TOTAL  THRUST 

5056 

THROTTLE  ANGLE 

46.4 

ENGINE  S.H.P. 

460 

UEIGHT  AND  BALANC 

IXX  INERTIA  (/  1024) 

728 

IZZ  INERTIA  (/  1024) 

1551 

PRESSURE  ALTITUDE 
CALIBRATED  AIRSPD 
EQUIUALENT  AIRSPD 
TRUE  AIRSPD  ( F/S  ) 
MACH  NUMBER 

BANK  ANGLE 
SIDESLIP 

RATE  OF  CLIMB  ( FPM  ) 
PITCH  ACCELERATION 
ROLL  ACCELERATION 
YAU  ACCELERATION 
NORTH-SOUTH  ACCEL 
EAST-UEST  ACCELERATION 
UERTICAL  ACCELERATION 
TOTAL  PITCHING  MOMENT 
TOTAL  ROLLING  MOMENT 
TOTAL  YAUING  MOMENT 

ELEUATOR  TRIM  TAB 
COLUMN  POSITION 
RUDDER  TRIM  TAB 
PEDAL  POSITION 
AILERON  TRIM  TAB 
UHEEL  POSITION 

THRUST  COEFFICIENT 

LATERAL  T.C. 

ENGINE  T.  I.T." 

E 

IYY  INERTIA  (/  1024) 

CROSS  PRODUCT/ INERT  I A 


00:01: 13 

524.0 
114.5 
114.53 
195. 12 
0.  18 

-1.8 
-0.2 
252 

0.0025 
-0. 0007 

0.0011 
-0.238 
-1.073 
-0.013 
2128 
-576 
1837 

12.69 
4.05 
3.79 
-0.34 
0.42 
2.  15 

0.09 
0.08 
550 

862 
43270 


NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL   FUEL 

ENGINE 

1 

2 

3     4 

BRAKES 

TANK    LOAD 

SHP 

484 

483 

500   597 

PILOT 

1    3917 

2   391? 

TIT 

549 

549 

551    564 

COPILOT 

3   3917 

FUEL  FLOU 

849 

849 

855   893 

FLIGH- 

r 

IAS 

114 

4   3915 

ALT 

"490 

5   996 

RUNU 

NAS  RU*i 

HDG 

255.8 

FLIGHT 

CONDITIONS 

PAGE 

FLIGHT  TIMER 

00: 00: 00 

MET  TIMER              00:05: 10 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

89880 

PRESSURE  ALTITUDE 

524.0 

C.G. 

30.00 

CALIBRATED  AIRSPD 

114.7 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

114.72 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 
MACH  NUMBER 

195.28 
0.  18 

FLIGHT/A 

PITCH  ANGLE 

8.  1 

BANK  ANGLE 

0.4 

ANGLE  OF  ATTACK 

6.2 

SIDESLIP 

-0.3 

HEADING  ANGLE 

272.9 

RATE  OF  CLIMB  ( FPM  ) 

407 

PITCH  UELOCITY  (D/S) 

-0.039 

PITCH  ACCELERATION 

0.0028 

ROLL  UELOCITY  (D/S) 

-0.086 

ROLL  ACCELERATION 

-0.0011 

YAU  UELOCITY  < D/S ) 

0.062 

YAU  ACCELERATION 

-0.0014 

NORTH-SOUTH  UELOCITY 

84.02 

NORTH-SOUTH  ACCEL 

0.438 

EAST-UEST  UELOCITY 

-176. 15 

EAST-UEST  ACCELERATION 

0.376 

UERTICAL  UELOCITY 

-6.78 

UERTICAL  ACCELERATION 

0.  123 

LONGITUDINAL  ACCEL 

0. 1567 

TOTAL  PITCHING  MOMENT 

2384 

LATERAL  ACCEL 

0.0002 

TOTAL  ROLLING  MOMENT 

-768 

UERTICAL  ACCEL  (G'S) 

-0.9844 

TOTAL  YAUING  MOMENT 

-2103 

CONTROL  LOADING 

ELEUATOR  POSITION 

-1.54 

ELEUATOR  TRIM  TAB 

12.55 

COLUMN  FORCE 

1.  19 

COLUMN  POSITION 

5.33 

RUDDER  POSITION 

-3.  10 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-5.31 

PEDAL  POSITION 

-0.34 

AILERON  POSITION 

0.  19 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.  14 

UHEEL  POSITION 

0.53 

ENGINE 

TOTAL  THRUST 

5056 

THRUST  COEFFICIENT 

0.09 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C. 

0.08 

ENGINE  S.H.P. 

460 

ENGINE  T.  I.T.' 

549 

UEIGHT  AND 

BALANCE 

I XX  INERTIA  (/  1024) 

725 

IYY  INERTIA  (/  1024) 

861 

I2Z  INERTIA  (/  1024) 

1548 

CROSS  PRODUCT/ INERT  I A 

43214 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRI 
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FUEL 

FUEL 

ENGINE 

1 

2 

3 

4 

BRAKES 

TANK 

LOAD 

SHP 

486 

485 

496 

601 

PILOT 

1 
2 

3911 
3911 
3910 

TIT 

550 

550 

551 

564 

COPILOT 

FUEL  FLOU 

848 

848 

853 

893 

FLIGHT 

3 

IAS    108 

4 

3908 

ALT   490 

5 

996 

RUNUAY-  100 

MOFFETT 

NAS  RUY 

32R 

HDG   252.4 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00: 00: 00 

MET  TIMER 

00:05:32 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

89856 

PRESSURE  ALTITUDE 

524.0 

C.G. 

30.00 

CALIBRATED  AIRSPD 

108.6 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

108.59 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (  F/S ) 
MACH  NUMBER 

185.09 
0.  17 

FLIGHT/l 

PITCH  ANGLE 

12.8 

BANK  ANGLE 

1.  1 

ANGLE  OF  ATTACK 

7.3 

SIDESLIP 

0.7 

HEADING  ANGLE 

269.5 

RATE  OF  CLIMB  (  FPM  ) 

1085 

PITCH  UELOCITY  (D/S) 

-0. 109 

PITCH  ACCELERATION 

0. 0002 

ROLL  UELOCITY  (D/S) 

0.070 

ROLL  ACCELERATION 

-0. 0003 

YAU  UELOCITY  (D/S) 

0.062 

YAU  ACCELERATION 

-0.0002 

NORTH-SOUTH  UELOCITY 

71.96 

NORTH-SOUTH  ACCEL 

0.255 

EAST-UEST  UELOCITY 

-169.59 

EAST-UEST  ACCELERATION 

-0.033 

UERTICAL  UELOCITY 

-18.07 

UERTICAL  ACCELERATION 

0.483 

LONGITUDINAL  ACCEL 

0.2229 

TOTAL  PITCHING  MOMENT 

80 

LATERAL  ACCEL 

-0.0112 

TOTAL  ROLLING  MOMENT 

-192 

UERTICAL  ACCEL  (G'S) 

-0.9587 

TOTAL  YAUING  MOMENT 

-276 

CONTROL 

LOADING 

ELEUATOR  POSITION 

-1.72 

ELEUATOR  TRIM  TAB 

12.42 

COLUMN  FORCE 

1.41 

COLUMN  POSITION 

5.23 

RUDDER  POSITION 

-2.91 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-4.37 

PEDAL  POSITION 

-0.32 

AILERON  POSITION 

1.38 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.  17 

UHEEL  POSITION 

4.28 

ENGIN 

TOTAL  THRUST 

5168 

THRUST  COEFFICIENT 

0.  10 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C. 

0.  10 

ENGINE  S.H.P. 

463 

ENGINE  T.  I.T.' 

550 

UEIGHT  AND 

BALANCI 

I XX  INERTIA  (/  1024) 

724 

IYY  INERTIA  (/  1024) 

861 

122  INERTIA  (/  1024) 

1548 

CROSS  PRODUCT/ INERT  I A 

43206 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL 
TANK 

FUEL 
LOAD 

1 

10000 

2 

10000 

3 

10000 

4 

10000 

5 

995 

ENGINE  1  2  3  4 

SHP  4744  4744  4744  4744 

TIT  1075  1075  1075  1075 

FUEL  FLOU  2467  2466  2466  2467 


RUNUAY-  100  MOFFETT   NAS  RUY  32R 


BRAKES 
PILOT 
COPILOT 


FLIGHT 

IAS    125 
ALT   490 
HDG   228.6 


FLIGHT  TIMER 

GROSS  UEIGHT 
P  P 

FLAP  POSITION 
GEAR  POSITION 


PITCH  ANGLE 
ANGLE  OF  ATTACK 
HEADING  ANGLE 
PITCH  UELOCITY  (D/S) 
ROLL  UELOCITY  (D/S) 
YAU  UELOCITY  (D/S) 
NORTH-SOUTH  UELOCITY 
EAST-UEST  UELOCITY 
UERTICAL  UELOCITY 
LONGITUDINAL  ACCEL 
LATERAL  ACCEL 
UERTICAL  ACCEL  (G'S) 


FLIGHT  CONDITIONS  PAGE 

00: 00: 00  MET  TIMER 

CONF I GURAT I ON/COND I T I ONS 


114208 
26.01 

19.  1 

1.0 


FLIGHT/AERO 


8.0 

4.6 
246.  1 
0.273 
-0.430 
-0.500 
-5.82 
-213.54 
-12.72 
0. 1443 
-0.0097 
-1.0163 


PRESSURE  ALTITUDE 
CALIBRATED  AIRSPD 
EQUIUALENT  AIRSPD 
TRUE  AIRSPD  (  F/S ) 
MACH  NUMBER 

BANK  ANGLE 
SIDESLIP 

RATE  OF  CLIMB  ( FPM  ) 
PITCH  ACCELERATION 
ROLL  ACCELERATION 
YAU  ACCELERATION 
NORTH-SOUTH  ACCEL 
EAST-UEST  ACCELERATION 
UERTICAL  ACCELERATION 
TOTAL  PITCHING  MOMENT 
TOTAL  ROLLING  MOMENT 
TOTAL  YAUING  MOMENT 


00:07:52 

524.0 
125.6 
125.59 
213.97 
0.  19 

-2.7 

-0.8 
759 
0. 0000 
0.0010 

-0.0048 

-1.923 

0.015 

-0.801 

-176 

1344 

-9720 


CONTROL  LOADING 

ELEUATOR  POSITION 

-2.78              ELEUATOR  TRIM  TAB 

10.69 

COLUMN  FORCE 

6.28              COLUMN  POSITION 

4.66 

RUDDER  POSITION 

-10.06             RUDDER  TRIM  TAB 

5.08 

PEDAL  FORCE 

-46.25             PEDAL  POSITION 

-1.  12 

AILERON  POSITION 

1.29              AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.56              UHEEL  POSITION 
ENGINES 

4.04 

TOTAL  THRUST 

29280             THRUST  COEFFICIENT 

0.42 

THROTTLE  ANGLE 

90.0              LATERAL  T.C.  . 

0.42 

ENGINE  S.H.P. 

4699             ENGINE  T. I.T. 
UEIGHT  AND  BALANCE 

1075 

I XX  INERTIA  (/  1024) 

1154              IYY  INERTIA  (/  1024) 

872 

IZ2  INERTIA  (/  1024) 

1984             CROSS  PRODUCT/ INERT  I A 

47766 

NOTE:  UALUES  INUALID 

DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET 

COLSNP  TRU 

234 

FUEL 

FUEL 

ENGINE 

1 

2 

3 

4 

BRAKES 

TANK 

LOAD 

SHP 

4720 

4722 

4722 

4720 

PILOT 

1 
2 

10000 
10000 
10000 

TIT 

1075 

1075 

1075 

1075 

COPILOT 

FUEL  FLOU 

2459 

2459 

2459 

2459 

FLIGHT 

3 

IAS    113 

M 

10000 

ALT   490 

5 

995 

RUNUAY-  100 

MOFFETT   NAS  RUY 

32R 

HDG   203. 1 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00: 00: 00 

MET  TIMER 

00:08:32 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

114208 

PRESSURE  ALTITUDE 

524.0 

C.G. 

26.01 

CALIBRATED  AIRSPD 

112.7 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

112.72 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  ( F/S ) 
MACH  NUMBER 

192. 12 
0.  17 

FLIGHT/i 

PITCH  ANGLE 

12.9 

BANK  ANGLE 

-0.5 

ANGLE  OF  ATTACK 

5.9 

SIDESLIP 

-2.3 

HEADING  ANGLE 

220.  1 

RATE  OF  CLIMB  ( FPM  ) 

1425 

PITCH  UELOCITY  ( D/S ) 

-0.297 

PITCH  ACCELERATION 

0.0022 

ROLL  UELOCITY  (D/S) 

0.  102 

ROLL  ACCELERATION 

-0.0014 

YAU  UELOCITY  (D/S) 

0.039 

YAU  ACCELERATION 

-0.0004 

NORTH-SOUTH  UELOCITY 

-92.79 

NORTH-SOUTH  ACCEL 

-0.024 

EA5T-UEST  UELOCITY 

-166.59 

EAST-UEST  ACCELERATION 

-0.234 

UERTICAL  UELOCITY 

-23.70 

UERTICAL  ACCELERATION 

0.929 

LONGITUDINAL  ACCEL 

0.2163 

TOTAL  PITCHING  MOMENT 

1872 

LATERAL  ACCEL 

0.0105 

TOTAL  ROLLING  MOMENT 

-1728 

UERTICAL  ACCEL  (G'S) 

-0.9460 

TOTAL  YAUING  MOMENT 

-710 

CONTROL 

LOADING 

ELEUATOR  POSITION 

-3.  15 

ELEUATOR  TRIM  TAB 

10.69 

COLUMN  FORCE 

8.00 

COLUMN  POSITION 

4.45 

RUDDER  POSITION 

-12.70 

RUDDER  TRIM  TAB 

5.08 

PEDAL  FORCE 

-48.78 

PEDAL  POSITION 

-1.41 

AILERON  POSITION 

0.75 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-0.83 

UHEEL  POSITION 

2.34 

ENGIN 

TOTAL  THRUST 

30176 

THRUST  COEFFICIENT 

0.54 

THROTTLE  ANGLE 

90.0 

LATERAL  T.C. 

0.54 

ENGINE  S.H.P. 

4671 

ENGINE  T. I.T. 

1075 

UEIGHT  AND 

BALANCI 

I XX  INERTIA  (/  1024) 

1154 

IYY  INERTIA  (/  1024) 

872 

I2Z  INERTIA  (/  1024) 

1984 

CROSS  PRODUCT/ INERT  I A 

47766 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL   FUEL 

ENGINE 

1 

2 

3     4 

BRAKES 

TANK    LOAD 

SHP 

4743 

4742 

PILOT 

1    10000 

TIT 

1075 

107E 

COPIL 

2  10000 

3  10000 

FUEL  FLOU 

2465 

2465   2465   2465 

FLIGHT 

IAS 

124 

4    10000 

ALT 

490 

5    995 

RUNU 

193.7 

HDG 

FLIGHT 

CONDITIONS 

PAGE 

FLIGHT  TIMER 

00: 00: 00 

MET  TIMER             00:09:06 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

114208 

PRESSURE  ALTITUDE 

524.0 

C.G. 

21.00 

CALIBRATED  AIRSPD 

124.2 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

124. 16 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  ( F/S  ) 
MACH  NUMBER 

211.66 
0.  19 

1 

rLIGHT/AE 

PITCH  ANGLE 

8.5 

BANK  ANGLE 

-1.2 

ANGLE  OF  ATTACK 

5.  1 

SIDESLIP 

-4.  1 

HEADING  ANGLE 

210.7 

RATE  OF  CLIMB  (FPM) 

757 

PITCH  UELOCITY  (  D/S ) 

0.086 

PITCH  ACCELERATION 

-0.0016 

ROLL  UELOCITY  (D/S) 

0.250 

ROLL  ACCELERATION 

-0.0037 

YAU  UELOCITY  (D/S) 

-0. 180 

YAU  ACCELERATION 

-0.0048 

NORTH-SOUTH  UELOCITY 

-136.83 

NORTH-SOUTH  ACCEL 

0.518 

EAST-UEST  UELOCITY 

-160.99 

EAST-UEST  ACCELERATION 

-0.250 

UERTICAL  UELOCITY 

-12.68 

UERTICAL  ACCELERATION 

-0.932 

LONGITUDINAL  ACCEL 

0. 1521 

TOTAL  PITCHING  MOMENT 

-1456 

LATERAL  ACCEL 

0.0393 

TOTAL  ROLLING  MOMENT 

-4160 

UERTICAL  ACCEL  (G'S) 

-1.0176 

TOTAL  YAUING  MOMENT 

-9577 

CONTROL  LOADING 

ELEUATOR  POSITION 

-4.99 

ELEUATOR  TRIM  TAB 

10.69 

COLUMN  FORCE 

16.50 

COLUMN  POSITION 

3.46 

RUDDER  POSITION 

-m.6i 

RUDDER  TRIM  TAB 

8.98 

PEDAL  FORCE 

-52. MM 

PEDAL  POSITION 

-1.62 

AILERON  POSITION 

0.  17 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.M1 

UHEEL  POSITION 

0.45 

ENGINES 

TOTAL  THRUST 

29376 

THRUST  COEFFICIENT 

0.43 

THROTTLE  ANGLE 

90.0 

LATERAL  T.C. 

0.43 

ENGINE  S.H.P. 

4693 

ENGINE  T. I.T. 

1075 

UEIGHT  AND  BALANCE 

I XX  INERTIA  (/  1024) 

1154 

IYY  INERTIA  (/  1024) 

872 

IZZ  INERTIA  (/  1024) 

1984 

CROSS  PRODUCT/ INERT  I A 

47766 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL 

FUEL 

ENGINE 

1 

2 

3 

4 

BRAKES 

TANK 

LOAD 

SHP 

4722 

4721 

4721 

4722 

PILOT 

1 
2 

10000 
10000 
10000 

TIT 

1075 

1075 

1075 

1075 

COPILOT 

FUEL  FLOU 

2459 

2459 

2459 

2459 

FLIGHT 

3 

IAS    113 

^ 

10000 

ALT   490 

5 

995 

RUNUAY-  100 

MOFFETT   NAS  RUY 

32R 

HDG    171.9 

FLIGHT  TIMER 

GROSS  UEIGHT 

C.G. 

FLAP  POSITION 

GEAR  POSITION 


PITCH  ANGLE 
ANGLE  OF  ATTACK 
HEADING  ANGLE 
PITCH  UELOCITY  ( D/S  ) 
ROLL  UELOCITY  (D/S) 
YAU  UELOCITY  (D/S) 
NORTH-SOUTH  UELOCITY 
EAST-UEST  UELOCITY 
UERTICAL  UELOCITY 
LONGITUDINAL  ACCEL 
LATERAL  ACCEL 
UERTICAL  ACCEL  (G'S) 

ELEUATOR  POSITION 
COLUMN  FORCE 
RUDDER  POSITION 
PEDAL  FORCE 
AILERON  POSITION 
UHEEL  FORCE 

TOTAL  THRUST 
THROTTLE  ANGLE 
ENGINE  S.H.P. 


I XX 
IZ2 


INERTIA 
INERTIA 


(/ 
(/ 


1024) 
1024) 


FLIGHT  CONDITIONS  PAGE 

00:00:00  MET  TIMER 

CONF I GURAT I ON/COND I T I ONS 


114208 
21.00 

19.  1 

1.0 


FLIGHT/AERO 


12.6 
6.8 
188.7 
-0.070 
0.273 
0.055 
-166.58 
-93.36 
-19.51 
0.2287 
0.0112 
-0.9812 

CONTROL 
-6.02 

15.87 

-14. 19 

-41.87 

0.90 

1.61 


LOADING 


PRESSURE  ALTITUDE 
CALIBRATED  AIRSPD 
EQUIUALENT  AIRSPD 
TRUE  AIRSPD  (F/S) 
MACH  NUMBER 

BANK  ANGLE 
SIDESLIP 

RATE  OF  CLIMB  (FPM) 
PITCH  ACCELERATION 
ROLL  ACCELERATION 
YAU  ACCELERATION 
NORTH-SOUTH  ACCEL 
EAST-UEST  ACCELERATION 
UERTICAL  ACCELERATION 
TOTAL  PITCHING  MOMENT 
TOTAL  ROLLING  MOMENT 
TOTAL  YAUING  MOMENT 

ELEUATOR  TRIM  TAB 
COLUMN  POSITION 
RUDDER  TRIM  TAB 
PEDAL  POSITION 
AILERON  TRIM  TAB 
UHEEL  POSITION 


ENGINES 


30208  THRUST  COEFFICIENT 

90.0  LATERAL  T.C.  . 

4674  ENGINE  T. I.T. 

UEIGHT  AND  BALANCE 

1154  IYY  INERTIA  (/  1024) 

1984  CROSS  PRODUCT/ INERT  I A 


00: 10:07 

524.0 
112.6 
112.56 
191.91 
0.  17 

-0.7 
-2.5 

1169 
-0. 0010 
-0.0012 
0.0032 
-0.232 
-0. 130 
-0.237 
-944 
-1600 
6515 

10.69 
2.96 
8.98 
-1.56 
0.42 
3.  13 

0.54 
0.54 
1075 

872 
47766 


NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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